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ABSTRACT 


A major  part  of  this  study  is  the  specification  of  a modem 
optimized  for  transmission  of  digitized  speech  over  a HF  channel. 

A tenth-order  linear  predictive  coder  (LPC-10) , yielding  a data 
rate  of  2400  b/s,  is  used  as  the  speech  processor  and  a HF  chan- 
nel with  an  available  bandwidth  of  2.4  kHz  is  assumed.  A modula- 
tion format  which  is  frame  synchronous  to  that  of  the  speech  pro- 
cessor allows  coding  to  be  applied  to  critical  bits  generated  by 
the  LPC-10  speech  analyzer.  In  addition  to  the  specification  of 
modems  for  digitizing  speech,  an  equally  important  part  of  this 
study  is  the  development  of  new  preambles  which  can  operate  at 
low  signal-to-noise  ratios  in  the  presence  of  time-varying  multi- 
path  profiles  typical  of  HF  channels.  This  task  includes  the 
specification  of  a waveform  for  signal  presence  detection,  Doppler 
estimation,  frame  estimation,  and  the  decoding  of  a KG  synchro- 
nizing sequence.  Complexity  estimates,  based  on  an  implementation 
by  a programmable  processor,  are  obtained  for  the  preamble  and 
digitized  speech  modem  specified  in  this  study. 
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SECTION  1 
INTRODUCTION 

A major  part  of  this  study  is  the  specification  of  a modem 
optimized  for  transmission  of  digitized  speech  over  a HF  channel. 

A tenth-order  linear  predictive  coder  (LPC-10) , yielding  a data 
rate  of  2400  b/s,  is  used  as  the  speech  processor  and  a HF  chan- 
nel with  an  available  bandwidth  of  2.4  kHz  is  assumed.  A modu- 
lation format  which  is  frame  synchronous  to  that  of  the  speech 
processor  allows  coding  to  be  applied  to  critical  bits  generated 
by  the  LPC-10  speech  analyzer.  In  addition  to  the  specification 
of  modems  for  digitizing  speech,  an  equally  important  part  of 
this  study  is  the  development  of  new  preambles  which  can 
operate  at  low  signal-to-noise  ratios  in  the  presence  of  time- 
varying  multipath  profiles  typical  of  HF  channels.  This  task 
includes  the  specification  of  a waveform  for  signal  presence 
detection,  Doppler  estimation,  frame  estimation,  and  the  decoding 
of  a KG  synchronizing  sequence.  Complexity  estimates,  based  on 
an  implementation  by  a programmable  processor,  are  obtained  for 
the  preamble  and  digitized  speech  modem  specified  in  this  study. 

An  effective  technique  for  signal  presence  detection  and 
Doppler  estimation,  based  on  using  four  unkeyed  parallel  tones,  is 
developed  in  Section  2.  Frame  detection  is  achieved  by  the  use 
of  Carley  sequences  composed  of  wideband  pulses  with  raised  cosine 
envelopes.  The  total  time  for  this  preamble  is  .449  second  which 
is  well  within  our  design  goal  of  0.72  second.  Also,  successful 
operation  is  possible  at  an  average  signal-to-noise  density  ratio 
of  40  dB-Hz.  This  corresponds  to  a signal-to-noise  ratio  measured 


in  a 2.4  kHz  bandwidth  of  only  6.2  dB.  Appendix  C discusses  var- 
ious definitions  of  SNR  and  develops  simple  relations  between  them. 

Once  Doppler  and  framing  estimates  are  obtained  a parallel 
tone  format  is  used  for  the  transmission  of  the  KG  sequences  and 
the  digitizing  speech  data.  A new  Doppler  tracking  algorithm,  based 
on  using  all  the  keyed  data  tones  is  also  developed  in  Appendix  E 
and  Section  5.  Frame  tracking  is  achieved  by  comparing  energy 
levels  in  an  early/late  empty  tone  slot,  as  is  done  in  conventional 
HF  modems.  It  should  be  noted  that  this  slot  sync  technique  for 
tracking  represents  the  only  conventional  technique  used  in  this 
ANDVT  modem.  The  inherent  in-band  diversity  obtained  from  this 
tracking  approach  justifies  its  use. 

In-band  time  and  frequency  diversity  combining  used  in  con- 
junction with  a strong  error  correcting  block  is  used  in  Section  3 
for  reception  of  the  KG  synchronization  sequence.  As  an  example, 
a combined  KG  sync  sequence  and  net  control  sequence  of  up  to  100 
bits  can  be  received  by  the  use  of  a (250,100;46)  block  code  with 
eighth-order  in-band  diversity.  Assuming  a 4-phase  DPSK  modem, 
with  39  tones  keyed  every  22.5  ms,  a transmission  time  of  26  frames 
yields  1014  tones,  of  which  1000  are  required  for  obtaining  eighth- 
order  diversity.  If  we  allow  one  reference  frame  and  an  additional 
three  frames  for  processing,  the  total  time  required  for  the 
reception  of  the  KG  sync  sequence  is  .675  second  (30  frames)  which 
is  well  within  the  design  goal  of  less  than  1.53  seconds  for  the 
reception  of  this  sequence.  Successful  reception  of  this  preamble 
is  obtained  at  a signal-to-noise  density  ratio  of  35  dB-Hz  even 
with  just  hard  decoding.  This  number  is  quite  compatible  with  the 
results  for  the  preamble  since  the  acceptable  peak-to-average  ratio 
for  this  39  tone  format  should  be  about  5 dB  greater  than  that 
required  for  the  preamble. 
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In  Section  4 specific  modem  designs  for  an  LPC-10  speech 
processor  are  evaluated  for  several  HF  channel  models.  In 
addition  to  obtaining  results  in  terms  of  bit  error  probabil- 
ities, speech  tapes  of  the  synthesizer  speech  after  transmission 
over  a HF  channel  are  obtained.  The  major  emphasis  has  been  on 
comparing  a conventional  16  tone  channel,  operating  at  a 75  f/s 
rate,  to  a 39  tone  modem,  operating  at  the  same  frame  rate  of 
the  speech  processor,  i.e.,  44.44  f/s.  The  39  tone  modem  uses 
two  (24,12;8)  Golay  codes  to  protect  24  critical  bits  in  each 
frame  of  54  bits.  The  performance  results  indicate  quite  clearly 
the  superiority  of  the  39  tone  modem  over  the  more  conventional 
16  tone  modem.  Acceptable  voice  communication  appears  possible 
at  an  E^/Nq  of  around  11  dB,  which  corresponds  to  a signal-to- 
noise  density  ratio  of  around  45  dB-Hz.  This  number  is  several 
dB  greater  than  required  for  the  successful  reception  of  preamble 
and  KG  sync  sequence.  This  desired  result  indicates  that  the 
performance  of  the  ANDVT  modem  will  not  be  limited  by  signal 
acquisition  requirements,  as  is  the  case  for  most  conventional 
HF  modems . 


Some  results  have  been  obtained  by  using  a 41  tone  and  a 
multiple  rate  coding  approach  which  operates  with  a decoding  delay 
of  10  frames  (.225  second).  The  preliminary  results  to  date 
indicate  that  the  simpler  Golay  coding  approach  is  superior  to 
this  multiple  rate  coding  technique.  However,  further  evaluation 
in  this  area  is  still  required  since  multiple  burst  correction 
codeshave  not  yet  been  evaluated,  and  the  HF  channel  models  have 
not  yet  been  modified  to  include  impulse  noise. 
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The  computational  complexity  of  the  modem  operations  is 
treated  in  Section  5 from  the  point  of  view  of  implementation 
on  a programmable  processor.  The  number  of  operations  for  each 
task  has  been  determined  for  a model  of  a quite  basic  processor 
which  certainly  does  not  exceed  the  capabilities  of  the  (undefined) 
Quintrell  processors.  Half-duplex  operation  appears  well  within 
the  capabilities  of  the  processor  model  assumed. 

Several  miscellaneous  topics  are  discussed  in  Section  6. 

These  include  signal  level  optimization,  a data  transmission 
mode  for  the  ANDVT,  and  a performance  monitoring  technique  for 
this  HF  modem. 

It  is  helpful  to  provide  an  example  of  a potential  half- 
duplex format  for  the  ANDVT  which  can  be  obtained  from  the 
detailed  results  given  in  the  remainder  of  this  report. 

A four-tone  preamble  has  been  chosen  for  signal  presence 
detection  and  Doppler  offset  estimation.  The  combined  time 
required  for  signal  presence  detection  and  Doppler  offset  esti- 
mation is  .359  second.  After  the  Doppler  offset  is  corrected, 
frame  estimation  is  achieved  by  the  use  of  wide  band  pulses  to 
form  an  aperiodic  sequence  with  good  autocorrelation  properties. 

A pair  of  Carley  sequences  with  appropriate  spacing  is  used  for 
framing.  A duration  of  .09  second  (four  frames  at  22.5  ms  per 
frame)  is  estimated  as  the  time  required  for  frame  estimation. 

Thus  the  duration  of  the  entire  acquisition  format  is  .449  second. 
Following  the  acquisition  format,  a parallel  tone  format  with 
39  tones,  keyed  at  a rate  of  44.44  f/s,  can  be  used.  The  modula- 
tor is  4-phase  differential  phase  shift  keying  (DPSK) . The  KG  sync 
sequence  and  any  net  control  information  can  be  obtained  in  a time 
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/ period,  including  the  processing  for  a (250,100;46)  code,  of 

.675  second.  After  the  reception  of  the  KG  sync  information, 

the  39  tone  format  may  be  used  with  two  Golay  codes  to  protect 

i 

the  critical  LPC-10  bits. 

Further  studies  and  detailed  simulations  are  expected  to 
follow  this  initial  effort.  During  this  next  phase  more  detailed 
results  will  be  obtained  and  other  formats  in  addition  to  that 
given  on  Fig.  1.1  will  be  evaluated. 
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Figure  1.1  Example  of  a Potential  Half-Duplex  Format 


SECTION  2 


PREAMBLE  DESIGN 


2 . 1 Introduction 

The  preamble,  provided  at  the  beginning  of  the  message 
just  prior  to  the  actual  transmission  of  data,  has  three  primary 
purposes . 

1)  To  assert  presence  of  the  signal.  This  is 
done  via  comparative  energy  measurements  made 
at  the  receiver.  The  operation  is  referred 
to  as  signal  presence  detection  (see  descrip- 
tion in  Section  2.2  and  detailed  analytical 
treatment  in  Section  2.4). 

2)  To  make  available  to  the  receiver  a signal 
that  can  be  processed  to  obtain  an  estimate 
of  Doppler  shift  and  hence  effect  an  initial 
Doppler  correction.  This  operation  is  referred 
to  as  Doppler  shift  estimation  (see  overview 
description  in  Section  2.2  and  detailed  analy- 
tical treatment  in  Section  2.5). 

3)  To  provide  the  receiver  with  a signal  that 
can  be  processed  in  such  a way  as  to  initiate 
frame  synchronization  (see  overview  in  Section 
2.2  and  detailed  discussion  in  Section  2.6). 

The  sections  to  follow  describe  preliminary  design  specifica- 
tions for  the  preamble  and  preamble  processor  at  the  receiver. 
Numerical  values  are  provided  for  the  preamble  parameters  that 
most  seriously  impact  operation.  Additionally,  detailed  descrip- 
tions of  all  processing  operations  are  provided.  The  numerical 
values  and  procedures  are  summarized  in  Section  2.2  which  is 
dedicated  to  a descriptive  overview  of  the  preamble;  the  design 
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guidelines,  general  structure,  signal  and  processor  timing,  and 
specification  of  numerical  values  for  relevant  parameters  are 
presented. 

Subsequent  sections  describe  the  performance  capabilities 
of  the  preamble  and  preamble  processor.  In  particular,  because 
part  of  the  preamble  consists  of  parallel  tones,  the  problem  of 
reducing  the  crest  factor  via  phasing  of  the  tones  is  discussed 
in  Section  2.3. 

The  signal  presence  detector  is  described  and  analyzed  in 
detail  in  Section  2.4.  Similarly,  performance  of  the  Doppler 
estimation  stage  is  analyzed  in  Section  2.5,  and  the  sync  pre- 
amble is  analyzed  in  Section  2.6. 
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2 . 2 Description  of  Preamble  and  Processor 

2.2.1  General  Structure 

The  total  time  allotted  to  the  preamble  is  divided  into 
three  sections  corresponding  to  the  three  major  tasks  of  the 
preamble.  Similarly,  the  tasks  are  performed  sequentially  at 
the  receiver.  In  the  order  performed,  these  tasks  are:  • 

1)  Detection  of  signal  presence 

2)  Doppler  estimation 

3)  Frame  sync  estimation. 


The  g T,Qral  timing  and  nomenclature  of  the  transmitted  preamble 
are  illustrated  in  Fig.  2.1.*  Note  that  the  same  signal  structure 


As  prespecified,  the  preamble  has  been  designed  so  that  its 
overall  length  is  no  more  than  .72  second  (32  frames  at  22.5 
ms/ frame) . 
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Figure  2.1  Timing  of  Transmitted  Preamble 


(four  tones  separated  one  from  the  other  by  approximately  600  Hz) 

is  used  for  the  detection  preamble  as  for  the  Doppler  section  of 

the  preamble.  The  tones  are  "on"  for  T . + seconds,  T . 

win  Dopp  win 

seconds  allotted  to  detection  of  signal  presence,  T^  seconds 
allotted  to  estimation  of  Doppler  shift. 

Because  of  dynamic  limitations  of  equipment  it  is  desirable 
to  limit  the  crest  factor  (peak-to-average  power  ratio)  of  the 
four-tone  composite  signal.  This  can  be  achieved  easily  by  + 1 
multipliers.  Specifically,  a phasing  sequence  (0,0,0, tt]  on  the 
tones  yields  a crest  factor  of  1.77.  If  no  consideration  is  given 
to  simplicity  of  implementation  the  "best"  tone  phases  that  were 
determined  via  numerical  procedures  described  in  [2.1]  can  be 
used.  The  sequence  given  there,  (0,1.81,1.81,0],  results  in  a 
crest  factor  of  1.525.  Specific  details  are  provided  in  Section 
2.3. 

Processor  timing  is  illustrated  in  Fig.  2.2.  Note  that  sync 
processing  is  not  allowed  to  commence  while  the  tones  are  "on"; 
hence,  there  is  a waiting  period  before  sync  processing  begins. 

The  specifics  are  clear  from  Fig.  2.2.  The  first  line  shows  the 
preamble  timing.  The  second  and  third  lines  show  the  timing  of 
the  preamble  processor  in  two  extreme  cases  of  interest  that  are 
respectively,  Case  I (detection  in  zero  time)*  and  Case  II  (detec- 
tion in  full  allotted  time).  The  most  probable  case,  of  course, 
is  that  detection  will  occur  with  some  intermediate  number  of 
samples . 


This  is  an  unattainable  extreme,  of  course;  the  closest  one 
can  get  to  this  extreme  is  through  a threshold  exceedance  being 
obtained  with  the  collection  of  a single  sample. 
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It  can  be  seen  from  the  figure  that  the  preamble  timing  is 

to  a certain  extent  dominated  by  the  requirement  that  the  sync 

processing  not  start  until  the  tone  transmission  has  ended.  From 

consideration  of  Case  I,  this  requires  a wait  of  T . seconds 

after  the  Doppler  processing  has  been  completed.  The  same  wait, 

in  the  event  the  individual  detection  event  requires  the  full 

allotted  time  (Case  II),  means  that  the  processor  will  have 

waited  T . seconds  after  the  tone  transmission  has  ended.  Thus, 
win 

transmission  of  the  signals  required  for  frame  synchronization 

must  continue  for  at  least  T . seconds  after  the  tone  trans- 

wm 

mission  has  ended. 

The  problems  relating  to  the  detection  of  signal  presence, 
Doppler  estimation,  and  frame  synch  estimation  are  analyzed  and 
described  in  detail  in  Sections  2.4,  2.5,  and  2.6,  respectively. 
The  design  implications  of  those  extended  developments  are  sum- 
marized respectively  in  Sections  2.2.2,  2.2.3,  and  2.2.4  below. 

2.2  2 Detection  of  Signal  Presence 

The  four-tone  structure  of  the  preamble  segment  that  is 
used  for  both  detection  and  Doppler  estimation  was  described 
above  in  Section  2.2.1. 

In  the  receiver,  a filter  is  centered  at  each  of  the  four 
tone  frequencies.  This  filtering  operation  can  be  done  by  fre- 
quency shifting  the  input  separately  so  that  each  tone  frequency 
appears  at  d.c.  and  uses  a complex  low  pass  filter.  If  recur- 
sive, this  low  pass  filter  must  be  second-order,  at  least,  and 
a double  pole  is  sufficient.  Since  the  a priori  range  of  Doppler 
shifts  is  ± 75  Hz,  we  specify  a double-tuned  filter 
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H(f)  = ~ (2.1) 

(1  + j2Trftc)Z 

where  t is  the  time  constant  of  each  stage.*  We  additionally 
c 

specify  a filter  half-power  frequency  of  75  Hz.  Generally,  when 
the  input  to  the  filter  in  (2.1)  is  white  noise,  it  can  be  deter- 
mined that  sampling  of  the  output  every  A = 3t  seconds  results 
in  adjacent  samples  that  have  a correlation  coefficient  of  .2 
(see  Section  2.5).  The  constraint  on  the  half-power  frequency 
results  in  a numerical  specification  of  the  value  of  A (see  Sec- 
tion 2.5).  This  value  is 

A = 4.1  ms  (2.2) 

In  order  to  detect  signal  presence,  the  filter  output  sam- 
ples are  squared  and  added  and  fed  to  a filter  which  averages 
over  several  successive  samples.  We  call  this  latter  filter  the 
signal  energy  filter.  Signal  presence  is  denoted  by  the  fact 
that  this  filter  output  exceeds  a threshold.  Because  of  the 
action  of  the  radio's  AGC,  which  is  not  under  control  of  the 
modem,  the  total  signal-plus-noise  power  appearing  at  the  modem 
input  will  not  vary  much.  Thus,  a fixed  threshold  cannot  be  used 
in  deciding  signal  presence.  It  is  necessary  to  have  other  fil- 
ters (we  specify  four)  centered  where  no  tones  are  transmitted  to 
estimate  noise  level.  The  same  computations  are  made  as  for  the 
signal  filters  with  the  summation  of  squares  and  a low  pass  aver- 
aging filter.  This  latter  filter  is  called  the  noise  energy  filter. 


Time  has  not  permitted  an  extensive  comparison  of  filter  shapes. 
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The  configuration  of  the  signal  presence  detector  is  shown 
in  Fig.  2.3.  If  we  denote  S^(4A)  as  the  4'th  time  sample  at  the 
output  of  the  k'th  tone  filter  and  N^(4A)  as  the  T'th  time  sample 
at  the  output  of  the  k'th  noise  filter  the  threshold  test  is 
given  by 


N 4 

L E \S  .(4A)12 

4=1  k=l  > 

__  - < a 

£ E \N  (4A)12 

4=1  k=l 


(2.3) 


with  threshold  exceedance  (>)  resulting  in  a detection  event. 

The  operation  is  implemented  as  a moving  window;  continuous 
sampling  of  the  arriving  signal  (or  noise)  energy  takes  place 
with  only  the  N most  recent  time  samples  being  used  in  the 
threshold  test.  If  no  signal  is  present  and  the  N most  recent 
time  sampling  instants  provide  samples  that  result  in  threshold 
exceedance  a false  alarm  event  takes  place.  If  the  same  result 
is  obtained  when  a signal  i_s  present  a detection  event  takes  place. 
This  latter  event  may  take  place  even  though  the  tones  are  "on" 
for  considerably  less  than  N x A seconds.  This  is  the  cause  of 
the  fundamental  uncertainty  of  signal  arrival  time  discussed  in 
the  previous  section. 

To  relate  to  the  discussion  of  the  previous  section  we  note 

that 

T . = N x A (2.4) 
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Presence  Detection 


An  extensive  analysis  of  the  signal  presence  detector, 
accompanied  by  performance  curves,  is  provided  in  Section  2.4. 
As  a tentative  design  goal  we  have  chosen  the  bounds 


PFA  £ 10 


PID  £ 10 


(2.5) 


where  is  the  probability  of  false  alarm  and  P_^  is  the  prob- 

FA  ID 

ability  of  incorrect  dismissal.  The  specified  performance  levels 
are  to  be  achieved  with  a signal-to-noise  density  ratio  of 


rf-  = 40  dB-Hz 

No 


(2.6) 


where  P is  the  total  power  in  the  four-tone  signal  and  Nq  is 
the  single-sided  noise  power  density.  The  results  are  based  on 
the  assumption  of  independently  fluctuating  tones.* 

As  can  be  determined  from  the  results  of  Section  2.5,  the 
design  specifications  above  are  met  with  the  numerical  parameter 
values  itemized  in  Table  2-1. 

2.2.3  Doppler  Shift  Estimator 

The  composite  four-tone  signal  used  for  detection  of  signal 
presence  is  used,  immediately  following  the  detection  event, 
for  estimation  of  Doppler  shift.  The  four  tone  filters  used  pre- 
viously are  U6ed  unchanged  as  the  initial  stage  of  a three-stage 
adaptive  Doppler  estimator  that  utilizes  the  "current"  estimate 


The  validity  of  this  assumption  is  questioned  and  discussed  in 
Section  2.4. 
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TABLE  2-1 


SIGNAL  DETECTOR  DESIGN  PARAMETERS 


Signal-to-Noise  Density  Ratio: 


77-  = 40  dB-Hz 

No 


Number  of  Time  Samples: 


N = 10 


Corresponding  Duration  of  Detection  Window: 


T . = .041  second 

win 


Achieved  Performance: 


PFA  ' 10 

Pj  - .25  x 10' 


Threshold  Setting: 


a = 2.65 
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of  Doppler  shift  to  effect  a Doppler  correction  and  a narrowing 
of  the  tone-filter  bandwidth;  this  allows  (via  admission  of  less 
noise)  an  improvement  in  the  subsequent  estimate. 

The  algorithm  for  one-stage  Doppler  estimation  is  illus- 
trated in  Fig.  2.4.*  Note  that  N in  this  context  represents  the 
number  of  times  that  the  output  of  each  tone  filter  is  sampled. 
This  number,  as  we  show  below,  typically  varies  from  stage  to 
stage. 

An  extensive  discussion  of  the  single-stage  Doppler  esti- 
mator illustrated  in  Fig.  2.4  is  provided  in  Section  2.5.1.  Per- 
formance results  are  also  provided  in  the  form  of  an  analytical 
expression  for  a,  the  rms  error  in  estimating  the  Doppler  shift. 

Multi-stage  estimation  of  Doppler  shift  is  discussed  in 
Section  2.5.2  and  performance  is  determined  via  a worst-case 
analysis  described  in  that  section. 

The  basic  design  goals  and  preliminary  specifications  of 
numerical  parameters  for  the  estimator  are  provided  in  Table  2-2. 

Summarizing  the  operation  of  the  three-stage  estimator  we 
note  that  16  samples  of  the  output  of  the  originally  configured 
tone  filters  (half-power  frequency  = 75  Hz)  are  used  to  make  an 
initial  estimate  of  Doppler  shift.  This  estimate  is  used  to 
effect  a Doppler  correction  and  the  tone  filter  bandwidths  ai.e 
narrowed  to  a half-power  frequency  of  16.2  Hz.  The  outputs  of 
the  adapted  filter  [same  basic  structure  given  in  (2.1)]  are  now 

* . 

Questions  relative  to  estimation  of  the  derivatives  x^,  y,, 

etc.  are  discussed  in  Section  2.5 
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IN-PHASE  AND 
QUADRATURE 
DIFFERENTIATORS 


Figure  2.4  Illustration  of  a Doppler  Estimation  Technique 


TABLE  2-2 


DOPPLER  ESTIMATOR  DESIGN  PARAMETERS 


Signal-to-Noise  Density  Ratio: 

- AO  dB-Hz 

N0 

Design  Goal  for  RMS  Doppler  Shift  Estimation  Error: 
a < .25  Hz 


Stage 

#1: 

N = 16 

A = 4.1  ms 

( f 1 / 2 = 75  Hz^ 

Stage 

#2: 

N = 3 

A - 19.0  ms 

(f^/2  = 16.2  Hz) 

Stage 

#3: 

N = 1 

A = 88.2  ms 

(f1/2  = 3.5  Hz) 

Performance  Achieved: 

ct  = .234  Hz 

Total  Processing  Time  (Including  Transient  Build-Up): 

T = .318  second 
tot 

N = Number  of  samples  taken  per  stage 
A s Time  between  samples  used  in  computation 
f 1 / 2 5 Half-power  frequency  of  adaptive  filter 


k 
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sampled  three  times  at  the  new  rate.  The  basic  adaptation  pro- 
cess described  above  is  repeated  resulting  in  tone  filters  nar- 
rowed to  a 3.5  Hz  half-power  frequency.  The  outputs  are  sampled 
once  in  the  final  stage. 

Note  that  the  design  goal  for  Doppler  shift  rms  estimation 
error  is  a < .25  Hz.  It  would  be  illusory  to  attempt  to  estimate 
Doppler  to  a resolution  much  finer  than  .25  Hz  because  of  the 
occurrence  of  Doppler  spread  on  the  channel.  Available  measure- 
ments of  HF  channel  Doppler  spread  [2.2]  indicate  that  207,  of 
the  measurements  (expressed  as  a sigma  width  of  the  Doppler 
power  spectrum)  exceed  this  value.* 

2.2.4  Frame  Synch  Estimation 

In  this  section  we  describe  the  detailed  structure  of  the 
frame  synchronization  preamble  occupying  the  second  seg- 

ment of  the  preamble  illustrated  in  Figs.  2.1  and  2.2. 

The  sync  preamble  consists  of  a pulse  train  designed  both 
for  high  crest  factor  and  low  autocorrelation  sidelobes.  The 
pulse  shape  used  in  this  pulse  train  is  assumed  to  be  of  the 
Nyquist  type,  having  normalized  value  unity  at  t=0  and  zero  value 
at  time  samples  equal  to  multiples  of  A.  Another  requirement  is 
that  the  pulse  be  bandlimited  so  that  its  spectrum  lie  within  the 
flat  portion  of  the  HF  terminal  equipment.  Further  details  and 
considerations  are  provided  in  Section  2.6.  Because  time  did  not 
permit  a detailed  comparison  of  various  pulse  shapes,  we  tenta- 
tively propose  a pulse  with  a raised  cosine  spectrum  of  band- 
width 2/A.  This  satisfies  the  Nyquist  requirement.  The  low 
pass  spectrum  is  illustrated  later  in  Fig.  2.17. 


Statistical  measurements  of  Doppler  spread  are  very  limited.  The 
measurements  referred  to  were  taken  on  a Fort  Monmouth,  New  Jersey 
to  Palo  Alto,  California  path. 
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The  remainder  of  this  section  relates  to  specifics  of 
the  pulse  train  timing  and  the  discrete  sequence  used  for 
synchronization . 


In  a general  context  the  sync  preamble  timing  is  illus- 
trated in  Fig.  2.5.  The  length-K^  segment  of  the  preamble  will 
be  referred  to  as  the  correlation  sequence.  This  discrete 
sequence  is  chosen  both  for  high  crest  factor  and  low  autocor- 
relation sidelobes.  The  specific  choice  of  a length-15  Carley 
sequence  (peak-to-sidelobe  ratio  = 15)  is  recommended  for  reasons 
described  below  in  greater  detail  in  Section  2.6.  As  shown  in 
Fig.  2.5  the  correlation  sequence  is  repeated  periodically  with 
period  K2  (K^  > 2K^  - 1) .* 

The  second  line  of  Fig.  2.5  illustrates  the  timing  of  the 
output  of  the  matched  filter  illustrated  in  Fig.  2.6.  The  cor- 
relation peaks  and  sidelobes  (non-overlapping)  are  seen  to  repeat 
every  K2  samples.  Note  that  an  observation  window  of  length  K2 
will  include  one  correlation  peak. 

For  our  initial  recommendation  we  simplify  the  processor 
illustrated  in  Fig.  2.6  by  eliminating  the  centroid  calculation 
and  replacing  it  with  an  operation  that  selects  the  largest  peak 
which  is  then  used  to  determine  frame  synch. 

■>v 

Satisfaction  of  this  inequality  is  required  to  guarantee  that 
the  output  of  the  discrete  matched  filter  consists  ideally  of 
repeated  segments  of  the  a periodic  correlation  function  with- 
out overlap  of  samples,  (see  Section  2.6). 
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Figure  2.6  Sync  Preamble  Processing 


The  sequence  of  events  defining  the  threshold  test  are  as 
follows:  Subsequent  to  the  wait  of  seconds  after  Doppler 

processing,  K„  samples  of  the  matched  filter  output  are  observed. 

^ 2 
The  receiver  then  formulates  the  decision  variable  |S(t^)|  /T 

S' 

i = 1,  where  S(t)  represents  the  received  signal  and  F is  an 

estimate  of  the  reference  noise  level  T given  in  (2.197).*  Each 
of  the  samples  are  tested  against  a threshold  p^.  The  largest 
peak  exceeding  the  threshold  is  chosen  as  the  correlation  peak 
to  be  used  for  synchronization. 

It  is  helpful  to  itemize  some  features  of  this  procedure: 

1)  Any  one  of  the  correlation  peaks  in  the  second  line  of 
Fig.  2.5  is  suitable  for  initiation  of  clock  signal  provided  the 
spacing  between  spikes  is  equal  to  an  integer  multiple  of  frame 
durations . 

2)  Consideration  of  Fig.  2.2  indicates  that  the  processing 

time  (subsequent  to  Doppler  processing  and  including  the  T^ 
second  "wait")  is  T ^ + (K^  x A)  + 6q,  where  6q  is  a safety 

margin  (to  allow  for  transients) that  is  small  enough  not  to  war- 
rant concern  in  these  preliminary  general  considerations. 

3)  The  actual  observation  (search  for  the  correlation  peak) 
is  made  over  x A seconds  subsequent  to  the  waiting  period. 

From  Fig.  2.5  it  is  clear  that  observation  over  samples  will 
yield  one  correlation  peak  in  the  matched  filter  output. 


Time  limitations  did  not  permit  an  evaluation  of  the  accuracy 
to  which  r can  be  estimated.  All  indications  point,  however, 
to  the  conclusion  that  this  quantity  can  easily  be  estimated  to 
sufficient  accuracy  to  have  negligible  impact  on  our  conclusions 
of  frame  sync  performance.  See  discussion  in  Section  2. 6. 2. 2. 
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4)  Because  of  periodic  repetition  n times  the  sync  preamble 
is  of  total  duration  n x (K^  x A)  and  (n-1)  x (K^  x A)  must 
exceed,  or  at  least  equal,  the  waiting  time  T ^ . This  require- 
ment is  clear  from  consideration  of  Case  II  in  Fig.  2.2  and  guar- 
antees, when  detection  requires  the  full  alotted  time  T ^ , that 
there  will  be  at  least  K2  samples  of  the  sync  preamble  available 
subsequent  to  the  "wait"  of  T ^ seconds,  i.e.,  there  will  be 
sufficient  time  available  for  observation  of  one  fully  developed 
spike  in  the  matched  filter  output. 

From  the  discussion  above  we  see  that,  though  the  total 
sync  preamble  is  long  due  to  the  n periodic  repetitions  of  the 
basic  probing  sequence,  the  observation  interval  is  shorter  by  a 
factor  of  n.  Reduction  of  the  observation  interval,  in  accord 
with  the  discussion  of  Section  2.6,  tends  to  minimize  the  effects 
of  noise. 

Considerations  discussed  above,  apart  from  guard  times  for 
transients,  indicate  that  data  transmission  can  begin  shortly 
after  the  end  of  the  sync  preamble.  The  requirement  that  the 
periodic  sync  preamble  be  synchronous  with  the  data  is  satisfied 


K2A  = mTp 


(2.7) 


where  m is  an  integer  and  T„  is  the  frame  duration  which  is 

r 

already  specified  to  be  22.5  ms. 
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With  the  choice  of  a Carley  sequence  we  have  = 

(see  Fig.  2.3  and  discussion  in  Section  2.6).  With  the 
further  choice  m = 2 and  i = 1 ms  we  have  = 45  which 
satisfies  the  requirement  > 2K^  - 1.  For  the  detec- 
tion window  of  duration  T . = 41  ms , item  4)  above 

win 

requires  that  (n  - 1)  x 45  > 41  which  is  satisfied  for 

n = 2 yielding  a sync  preamble  of  duration  2 x 45  = 90  ms 

in  length.  Note  that  this  also  satisfies  the  requirement 

that  at  least  = 45  matched  filter  samples  be  available 

for  processing  after  the  T . = 41  ms  "wait". 

v b win 

Thus,  the  total  duration  of  the  sync  preamble  is  given 
by 

Total  Duration  of  Sync  Preamble  = .09  sec  (2.8) 

The  above  considerations  and  consideration  of  Fig.  2.5 
indicate  that  data  transmitted  15  ms  after  the  last  non- 
zero pulse  in  the  sync  preamble  will  be  frame  synchronous 
with  the  clock  signal  initiated  by  the  output  of  the  dis- 
crete matched  filter. 

To  identify  the  initial  frame  boundary  requires  that 
the  correlation  peak  be  identified  as  being  due  to  the  first 
or  second  period  of  the  correlation  sequence.  The  ambiguity 
can  be  removed  by  processing  for  the  second  correlation 
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peak;  this  secondary  processing  can  take  place  by  threshold 
testing  a very  small  number  of  matched  filter  outputs 
(ideally  just  one)  because  of  the  a-priori  fixed  time  inter- 
val between  the  peaks. 

The  length  of  the  observation  interval  is  x A = 

45  ms.  The  length  of  this  interval  has  importance  rela- 
tive to  the  effects  of  noise  as  discussed  in  Section  2.6. 

The  numerical  parameters  describing  the  sync  preamble  and 
processor,  as  determined  from  the  considerations  of  Section 
2.6,  are  summarized  in  Table  2-3. 

To  be  specific  we  have  described  just  one  way  of 
implementing  the  sync  preamble  timing.  There  are  a var- 
iety of  different  approaches  one  could  use.  As  an  example 
one  can  conceive  of  signal  processing  aimed  at  detecting 
the  end  of  the  tone  transmission.  This  could  consider- 
ably reduce  the  forementioned  timing  uncertainties  and 
to  some  degree  perform  the  function  of  the  sync  preamble 
itself. 

2.3  Crest  Factor  Reduction 


The  parallel  tones  that  constitute  the  Doppler  pre- 
amble can  be  represented  mathematically  in  terms  of  the 
complex  envelope  D(t) . Since  the  N tones  are  harmonically 
related  we  have 


TABLE  2-3 


FRAME  SYNC  ESTIMATOR  DESIGN  PARAMETERS 

Signal-to-Noise  Density  Ratio: 

— = 40  dB-Hz 
N0 

Sampling  Interval: 

A = 1 ms 

Synchronization  Sequence: 

Length-15  Carley  Sequence  Repeated  for 
Two  Periods  with  Period  = 45 

Duration  of  Sync  Preamble  (Ignoring  Transients) : 
.09  second 

Correlation  Peak  Threshold  Setting: 

P0  = 15 

Start  of  Data  Transmission: 

15  ms  After  Last  Sync  Pulse 


It 


D(t)  = E ej9neJ2»nFt 


(2.9) 


The  set  of  phases  0 must  be  chosen  to  reduce  the  crest  factor 

n 

(peak/average  power)  of  the  resultant  transmitted  waveform. 

The  squared  envelope  of  D(t)  is  given  by 


0 N N j 0 -j0  , v 

|D(t>|2  - E S e ne  V2’(n-m)Ft 

n=l  m=l 


+(N-1) 

I>  Re 
p=- (N-l)  P 


j2irpFt 


(2.10) 


where 


N-p  J 0.  _ -J  0 
„r  J m+p  m n 

L e re  p > 0 


R = 
P 


(2.11) 


N j 0,  - j 0 

rn  i p m r, 

Zj  e re  p < 0 

m=l-p 


Note  that  R^  is  the  autocorrelation  function  of  the  phase  weight- 
ing sequence.  If  this  autocorrelation  function  had  zero  side- 
lobes,  the  crest  factor  would  be  unity. 
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The  case  of  interest  for  us  is  N=4.  A reasonable  simple 
estimate  for  the  desired  phasing  sequence  can  be  generated  by 
using  + 1 for  the  set  [e^n]  (see  Frank  [2.3]) 

[0n]  = [0,  0,  0,  ir]  (2.12) 

The  correlation  function  for  this  sequence  is  given  by 

(Rpl  = (-1,  0,  1,  4,  1,  0,  -1]  (2.13) 

and  the  crest  factor  p which  is  defined 


= Peak  \D(t) \ 
average  [D(t)  \‘ 


= peak 


Imp. 


is  given  by 


p = peak  (1  + ^ (cos  2irFt  - cos  6irFt)} 


(2.14) 


(2.15) 


After  some  manipulation  it  can  be  demonstrated  that  the  peak 
occurs  at 


2irFt 


(2.16) 


so  that 


p = 1.77 


(2.17) 
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A smaller  crest  factor  can  be  obtained  by  using  the  "best 
phases  discovered"  reported  by  Nuttall  12.1].  As  determined  by 
numerical  search  techniques  the  tone  phases  are  given  by 

0n  = [0,  1.81,  1.81,  0]  (2.18) 

and  result  in  a crest  factor  value 

p = 1.525  (2.19) 

The  resulting  waveform  envelope  is  shown  in  Fig.  2.7  (taken  from 

[2.1]). 

Though,  to  be  specific,  we  recommend  use  of  the  second 
sequence,  its  use  represents  a gain  of  .64  dB  over  use  of 
the  first;  in  view  of  the  simplicity  of  the  first  sequence 
(the  last  tone  is  simply  reversed  in  polarity  relative  to  the 
first  three)  consideration  of  its  use  should  not  be  dropped 
completely . 

2.4  Performance  of  Signal  Presence  Detector 

The  structure  of  the  signal  presence  detector  was  presented 
earlier  in  Fig.  2.3.  It  should  be  recalled  that  the  decision 
variable  was  chosen  to  be  the  ratio  of  tone  filter  energy  output 
to  noise  filter  energy  output  in  order  to  avoid  AGC  effects. 

When  tones  are  being  received,  the  output  of  the  k'th  tone 
filter  w^(t)  may  be  represented  as  the  sum  of  signal  and  noise 

wk(t)  = zk(t)  + nk(t)  (2.20) 
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where  z^(t)  represents  the  received  tone,  and  n^t)  is  the  noise 
at  the  tone  filter  output.  At  the  noise  filter  output  we  have 
noise  only  [designated  r>k(t)]  . 

The  threshold  test  is  given  by 


N K 
£ £ 
t=l  k=l 


|wk(^T) 


N 

£ 


K 


£ \%(IV 

1=1  k=l  K 


H1 

> 

< 

HO 


a 


(2.21) 


where  is  the  hypothesis  "signal  present",  HQ  is  the  hypo- 
thesis "signal  absent",  and  a is  a preset  threshold.  The  sum- 
mation in  the  numerator  is  over  the  K tone  filters  and  N time 
samples,  whereas  the  summation  in  the  denominator  is  over  the 
K noise  filters  and  N time  samples. 

From  (2.21)  the  probability  of  false  alarm  is  given  by 


Pr 


N K 

£ £ |n,(tT)l 

1=1  k=l 


N K 

£ £ | T).(  -IT ) | 

<1=1  k=l  K 


while  the  probability  of  incorrect  dismissal  is 


(2.22) 
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P = Pr 
ID 


N K 

E E |z.(<,T)  + n ,(4T)  |Z 

1=1  k=l 


S S |n.(W) 
-1=1  k=l 


(2.23) 


To  calculate  these  quantities  we  first  treat  the  signal  as  a 
deterministic  quantity,  i.e.,  we  evaluate  the  probabilities  P ^ 
and  P^-p  conditional  on  the  set  of  values  z^(-tT).  Some  simpli- 
fication is  obtainable  since  we  can  assume,  by  design,  that  our 
sample  interval  exceeds  the  correlation  time  of  the  noise,  i.e., 
the  noise  sample  covariance  takes  the  form 


|n,(fT)|2  = |r,.  (IT)\2  = P, 


n£(-tT)ns(rT)  = n*(fT)ns(rT)  =0  ; 'k  i s 


rj*(fT)ns(rT)  = 0 


(2.24) 


(2.25) 


(2.26) 


These  equations  say  simply  that  the  noise  samples  from  any 
two  distinct  filters  are  always  uncorrelated  but  of  equal  energy. 
Coupled  with  the  assumption  of  complex  Gaussian  statistics  for 
the  noise,  xve  have  shown  that  the  ratio  in  (2.22)  is  a special 
case  of  the  central  F distribution  while  with  z^(4T)  fixed  the 
ratio  in  (2.23)  is  a special  case  of  the  non-central  F distri- 
bution. The  specific  forms  of  interest  to  us  can  be  found  in 
[2.4,  pp.  9 and  13].  For  the  probability  of  false  alarm,  we  get 
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W 


Thus , we  concentrate  on  computing  the  density  function  of  X and 
averaging  P-^Ce)  with  respect  to  X to  obtain  P^.  With  our 
assumption  that  the  tone  does  not  fluctuate  during  the  signal 
detection  interval, 

X = 2 \zk \ 2 (2.30) 

cr  k=l 
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where  we  have  dropped  the  argument  (IT).  Assuming  that  the 
[z^l  are  complex  Gaussian,  (2.30)  represents  a quadratic  form 
in  complex  Gaussian  variables.  The  characteristic  and  density 
functions  of  quadratic  forms  in  complex  variables  have  been 
derived  by  Turin  [2.5]. 

We  denote  the  density  and  characteristic  functions  of  X 
by  W(X)  and 

F ( j t)  = ejtX  = J ejtXW(X)  dX  , (2.31) 

respectively,  In  carrying  out  the  average  of  (2.28)  over  X,  we 
note  that  the  following  typical  average  occurs: 

Ir  - J Xr  exp  j-  2(Aa)  j »(X)  dX  (2.32) 

Comparing  (2.31)  and  (2.32),  we  see  immediately  that 


= S F(1t) 


b(jt)J 


= F 


(r) 


jt  = -l/2(l+«) 


(-2<lW)  (2-33) 


Then,  using  (2.33)  in  the  average  over  (2.28),  we  see  that 


KN-1 

£ 

j-o 


KN-1 

£ 

l=j 


1 iKN  + l - l\ 

j!2j  ^ i I 


/ a _] 

j 

( 1 

r i 1 

\l  +ct> 

ii+J  F 

2(1  + a) J 

(2.34) 
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p = 1 11 
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For  a generalized  random  variable  with  characteristic  function 
F(jt)  this  is  the  general  result.  In  our  application,  however, 
Z^(t)  is  assumed  to  be  a complex  Gaussian  random  variable.  Thus, 
X in  (2.29)  is  a quadratic  form  in  complex  Gaussian  random  vari- 
ables. From  the  work  of  Turin  [2.5],  the  characteristic  function 
of  a quadratic  form  in  complex  Gaussian  variables  is  given  by 

= det  [I  - jtMQ]  (2‘35> 

where  Q is  the  matrix  of  the  quadratic  form  q 

| 

q = ZHQZ  (2.36) 

(the  superscript  H denotes  conjugate  transpose),  Z is  a column 
vector  of  complex  Gaussian  variables 

Z1 
Z2 

LZk 

and  M is  the  moment  matrix  of  the  complex  Gaussian  variables 


M = (Z  - Z)(Z  - Z)H 


(2.38) 


The  elements  of  the  matrix  Q follow  immediately  from  (2.30). 
They  are 


(2.39) 


2-32 


"1 


To  determine  the  moment  matrix  we  assume  that  the  HF  chan- 
nel is  characterized  by  a discrete  set  of  M independently- 
fluctuating  paths  having  delays  £m,  Doppler  shifts  y , and  com- 
plex amplitudes  8m(t) . The  transfer  function  corresponding  to 
this  channel  is 


M — j 2 7T  f 4 j2  7tv  t 

T(f,t)  = £ gm(t)e  me  m (2.40) 


Consequently,  the  complex  correlation  between  the  k'th  and  t'th 
received  tone  is  just 


M J27T(fk  " ft^m 

Z(.(t)z£(t)  ' S Gme  k S(fk-£t)  (2.41) 


where  f^  is  the  frequency  location  of  the  k'th  tone, 

Gm-|gm(t)l2  (2.42) 

is  the  strength  of  the  m'th  path,  and  S(f)  is  the  frequency  cor- 
relation function  of  the  channel, 


S(f) 


M 

E 

m=l 


G e 


j 2fff£ 


m 


m 


(2.43) 
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Thus,  the  elements  of  the  moment  matrix  M are  determinable 
from  knowledge  of  the  frequency  correlation  function  of  the  chan 
nel.  If  it  can  be  assumed  that  the  received  tones  fluctuate 
independently,  the  above  formulation  can  be  greatly  simplified; 
the  moment  matrix  then  becomes  diagonal 

M = PQI  (2.44 

where  I is  the  identity  matrix  and  is  the  average  received 
tone  strength.  In  such  a case, 

F(jt)  = det  [I-jtPQNI] 

Manipulation  of  this  expression  indicates  that  the  charac- 
teristic function  derivatives  to  be  used  in  (2.34)  are  given  by 

*.(n)  I I \ _ K(K  + 1)  •••  (K  + n-1) 

Npo  f + n 

2<j2  (1  +a) 

Utilizing  the  above  formula,  the  signal  presence  detection 
capability  of  the  preamble  was  determined.  Numerical  results,  dis- 
played as  the  probability  of  incorrect  dismissal  vs.  number  of 
time  samples  N are  provided  in  Figs.  2.8  to  2.11.  Each  curve  is 
parameterized  to  a fixed  false  alarm  probability  and  each  figure 
relates  to  one  particular  value  of  signal-to-noise  density  ratio, 
P/N0. 

To  get  some  idea  of  the  processing  time  required  consider  a 
sampling  interval  of  5 msec.  (This  particular  choice  approximates 
more  highly  motivated  choices  relating  to  the  Doppler  processing 
stage  as  discussed  in  Section  2.5.)  Note  that  20  time  samples 
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(i.e.,  N = 20)  corresponds  to  .1  second  of  processing  time.  For 
the  curves,  four  tones  (K  = 4)  have  been  used  and,  as  pointed  out 
above,  have  been  assumed  to  fluctuate  independently. 


The  assumption  of  independence  is  one  that  should  be  exam- 
ined more  closely.  This  can  be  done  by  appealing  to  the  represen- 
tation of  the  frequency  correlation  function  in  (2.44)  and  expand- 
ing the  exponential.  Retaining  terms  up  to  second-order  and  set- 
ting the  mean  delay 


M 


£ 

m=l 


G 4 
m^m 


(2.47) 


equal  to  zero  we  define  the  frequency  correlation  coefficient 


p ( f ) = §XJL1 

' S(0) 


and  the  rms  delay  spread 


/m  ?a/2 

S G f 

m 


V 


m=l 


m^m 


M 

E G. 


\ T m / 

' m=l  / 


(2.48) 


(2.49) 


and  obtain  after  some  simple  substitutions  the  approximation 


p ( A)  **  1 - 2tt2A2cj2  (2.50) 

which  holds  when  A is  not  much  greater  than  zero  (i.e.,  there  is 
still  substantial  correlation).  Thus,  to  have  enough  decorrela- 
tion such  that  the  inequality 

P < P0  (2.51) 
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is  satisfied  the  inequality 


I' 


i - irA  ~ p0} 


(2.52) 


For!  the  preamble  tone  separation  (about  .6  KC)  the  above 
expansion  indicates  that  decorrelation  to  a value  of  p less  than 
.75  requires  a multipath  spread,  2<j^,  in  excess  of  .36  ms.  This 
is  typical  of  multipath  spreads  on  HF  links.  In  fact,  on  a link 
for  which  we  have  data,  the  Fort  Monmouth-to-Palo  Alto  path,  [2.2] 
there  is  a 607,  chance  that  the  cited  value  will  be  met  or  exceeded. 


To  get  a complete  picture  of  the  consequences  of  correla- 
tion we  have  calculated  detection  results  for  the  situation  wherein 
all  the  tones  fade  in  step.  For  this  case  we  have 


X 


and 


(2.53) 


i 

(2.54) 


The  characteristic  function  of  the  one  term  quadratic  form  is 
obtained  simply  from  (2.46)  by  setting  K = l.  This  case  is  truly 
a worst  case  situation  and,  as  one  would  expect,  performance  is 
considerably  inferior  to  that  presented  earlier.  This  is  illus- 
trated in  Fig.  2.12.  For  this  degenerate  case  it  appears  to  be 
roughly  true  that  the  processing  time  required  to  obtain  a per- 
formance level  comparable  to  that  in  Fig.  2.9  is  about  four  times 
as  great  as  it  was  for  the  independent  tone  case. 
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The  true  case  is  somewhere  inbetween  the  two  extreme  cases 
presented  with  a definite  bias  towards  the  earlier  high  perform- 
ance situation.  It  would  be  of  interest  to  examine  the  possible 
case  wherein  there  is  a high  correlation  between  adjacent  tones 
only.  Time  does  not  permit  such  considerations,  here.  It  can  be 
reasonably  speculated,  however,  that  high  correlation  between 
adjacent  tones  would,  to  obtain  the  same  performance,  require  a 
processing  time  falling  somewhere  inbetween  the  time  Tq  to  pro- 
cess uncorrelated  tones  and  2Tq. 

Though  the  specific  threshold  settings  required  to  obtain 
the  performance  in  Figs.  2.8  to  2.11  are  of  no  interest  from  the 
viewpoint  of  performance,  they  are  of  interest  relative  to  design. 
The  values  of  a are  tabulated  in  Table  2-4. 

2.5  Doppler  Estimation  Performance 

The  four-tone  structure  described  for  the  detection  process 
in  the  previous  section  is  also  used  for  Doppler  estimation.  As 
described  in  Section  2.2.1  the  tones  are  on  for  a total  time 


T = T . + T 

tone  win  Dopp 


(2.55) 


where  T is  the  duration  of  the  moving  window  used  in  the 

detection  process*  and  Tjjq  is  the  fixed  time  interval  required 

to  do  the  Doppler  estimation  as  discussed  later  in  this  section. 

With  high  probability  detection  will  take  place  in  T seconds 

where  T < T . . The  time  T1  remaining  in  the  tone  is 

Det  win  & 


The  duration  referred  to  is  just  the  amount  of  time  it  takes  to 
collect  the  N time  samples  of  the  energy  comparison  detector  in 
Fig.  2.  . 
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TABLE  2-4 


THRESHOLD  VALUES  (a) 


Vfa 

2KNsv 

io"2 

10"3 

10-4 

io-5 

io-6 

8 

6.03 

12.05 

22.71 

41.64 

75.3 

16 

3.37 

5.20 

7.61 

10.79 

15.01 

24 

2.66 

3.74 

5.00 

6.51 

8.32 

32 

2.32 

3.09 

3.95 

4.92 

6.02 

40 

2.11 

2.73 

3.38 

4.09 

4.88 

48 

1.98 

2.49 

3.02 

3.59 

4.20 

56 

1.88 

2.32 

2.77 

3.24 

3.74 

64 

1.80 

2.19 

2.59 

2.99 

3.42 

72 

1.74 

2.09 

2.44 

2.80 

3.17 

80 

1.69 

2.01 

2.33 

2.65 

2.98 

88 

1.65 

1.95 

B 

2.53 

2.82 

* 
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T ' = T -T  = T -T  +T 

tone  Det  win  Det  Dopp 


(2.56) 


After  Doppler  processing  the  time  remaining  is 


T"  - T . - Tn  . < T . 

win  Det  — win 


(2.57) 


seconds  and  hence,  one  must  wait  this  length  of  time  before  com- 
mencing sync  estimation. 

Some  of  the  issues  raised  above  will  be  discussed  in  the 

later  section  (Section  2.7)  on  sync  estimation.  In  this  section 

we  concentrate  on  determining  performance  of  the  Doppler  shift 

estimator  and  on  the  fixed  amount  of  time  that  must  be 

Dopp 

allocated  to  the  Doppler  estimation. 

2.5.1  Doppler  Estimation  Error  for  Single  Stage 

Doppler  estimation  takes  place  immediately  after  signal 
presence  detection,  i.e.,  immediately  after  the  threshold 
exceedance  indicated  in  (2.3).  Based  upon  the  analysis  in  [2.6, 
Section  5]  one  may  postulate  near-optimal  estimation  algorithms 
for  utilizing  the  signal  filter  outputs  to  estimate  Doppler  shift, 
In  the  reference,  the  following  algorithm  was  found  to  be  near- 
optimal  : 


K . 

k=l  K 


(2.58) 
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where  and  0^  denote  the  envelope  and  phase  derivative  of  the 
signal  at  the  k'th  tone  filter  output,  respectively,  and  K is 
the  number  of  tones.  A form  equivalent  to  (2.58)  is: 


v = 


2ir 


K 

L 

k=l 


Vk 


Vk 


K 2 2 

E X|  + y. 

k=l  k k 


(2.59) 


where  and  y^  are  the  in-phase  and  quadrature  components  and 

• • 

x^  and  y^  are  the  corresponding  derivatives  of  the  k th  tone  fil- 
ter output,  respectively.  This  complex  representation  assumes  a 
reference  tone  frequency  for  each  transmitted  tone. 

Equations  (2.58)  and  (2.59)  yield  an  estimate  of  Doppler 
corresponding  to  a single  time  sample  of  envelopes  and  phase 
derivatives  (or  in-phase  components,  quadrature  components,  and 
their  derivatives)  at  the  output  of  the  signal  filter.  However, 
several  time  samples  could  be  used  in  a Doppler  estimate.  Since 
the  fading  will  be  slow  compared  to  the  estimation  time,  succes- 
sive tone  filter  output  samples  will  differ  primarily  in  that 
the  noise  components  will  fluctuate.  Thus,  the  utility  of  having 
several  output  samples  is  primarily  one  of  allowing  the  noise  to 
be  averaged  out.  Thus  we  may  modify  (2.58)  and  (2.59)  to 


v = 


N K 

S S \l6kl 
1=1  k=l 


2n 


N K 
L L 
1=1  k=l 


(2.60) 
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where  the  added  l subscript  denotes  the  -t'th  time  sample  and  N is 
the  number  of  time  samples.  Note  that  the  denominators  in  (2.60) 
and  (2.61)  would  be  identical  to  the  output  of  the  signal  energy 
filter  if  integrate-and-dump  implementation  were  used  over  N samples. 

Since  the  tone  filters  must  have  a bandwidth  of  the  order  of 
150  Hz  and  since  the  tones  may  be  ± 75  Hz  relative  to  nominal 
center  frequencies,  a sampling  rate  somewhat  faster  than  150  samples/ 
second  would  be  desirable  for  the  filter  output.  However,  there 
is  an  added  complication  in  the  fact  that  it  is  necessary  to  com- 
pute  {0^}  or  [x^  , y^}  from  the  filtered  tone.  This  would  seem 
to  require  a sampling  rate  high  compared  to  75  Hz,  or  at  least 
several  samples  close  together  in  time  at  a high  rate. 

For  example,  the  simplest  approximation  to  the  derivative  of 
x(t)  is  given  by 


i(t) 


x(t  + t) 

T 


x(t) 


(2.62) 


Here  r must  be  small  compared  to  the  correlation  time  of  x(t)  or 
small  compared  to  one  cycle  of  maximum  Doppler  of  75  Hz.  On  the 
other  hand,  it  must  not  be  made  too  small  or,  due  to  the  ambient 
noise,  the  derivative  x(t)  will  be  estimated  poorly.  More  accurate 
algorithms  than  (2.62)  are  available  to  estimate  the  derivatives 
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which  involve  several  delayed  samples.  In  any  case,  an  analysis 
and/or  simulation  is  required  to  fully  optimize  the  selection  of 
algorithms.  Figure  2.13  presents  a block  diagram  of  the  process- 
ing described  above  to  estimate  Doppler  shift. 

Generally  it  can  be  said  that  the  problem  here  is  that  the 
number  of  time  samples  N to  be  used  in  estimating  the  Doppler 
shift  v must  be  optimized.  On  the  one  hand,  making  N large  reduce 
the  effect  of  noise  fluctuations.  On  the  other  hand,  making  N 
too  large  can  result  in  an  averaging  time  over  which  the  tone 
fluctuation  becomes  noticeable.  This  would  violate  the  ground 
rules  on  which  the  algorithm  was  developed.  The  above  esti- 
mation procedure,  as  discussed  at  the  end  of  this  section, 
is  modified  to  become  adaptive  by  using  the  Doppler  estimate 
to  make  a Doppler  correction  and  subsequently  reducing  the 
bandwidth  of  the  tone  filters  to  improve  the  estimate.  Before 
making  this  step,  however,  future  developments  are  simplified 
by  determining  the  performance  of  a Doppler  estimator  based 
on  (2.60)  or  (2.61). 

Equations  (2.60)  and  (2.61)  present  two  equivalent  versions 
of  the  Doppler  estimator  of  interest,  and  Eqs.  (2.58)  and  (2.59) 
present  estimates,  previously  evaluated  in  [2.6],  based  upon  a 
single  sample  of  the  filter  outputs.  Due  to  the  difficulty  of 
an  exact  analysis,  we  present  here  an  approximate  performance 
analysis  of  the  multi-sample  algorithm  based  upon  a high  SNR 
approximation.  This  approximation  is  useful  here  because  Doppler 
estimation  takes  place  only  after  signal  detection  has  occurred, 
i.e.,  after 


IN-PHASE  AND 
QUADRATURE 
DIFFERENTIATORS 


llustration  of  a Doppler  Estimation  Technique 


(2.63) 


L 


N K 
£ £ 
1=1  k=l 


|wk(tT) 


N K 
> a £ £ 

l-l  k=l 


\vkwn 


where 


wk(t)  = zk(t)  + nk(fc) 


(2.64) 


In  this  expression  z^(t)  represents  the  received  tone  and  n^(t) 


is  the  noise  at  the  filter  output.  The  rms  bandwidth  (Doppler 
spread)  of  zk(t)  is  assumed  small  enough  so  that,  over  the  time 
interval  for  Doppler  estimation,  zk(t)  differs  little  from  a 


cissoid . 


Mathematically,  in  the  vicinity  of  the  time  instant  t^. 


we  may  represent  zk(t)  in  the  form 


, M _ f . j0k(tO)  j(1:  ‘ t0)9k(t0)  . 

zk(t)  = Rk(tQ)e  e e (2.65) 


where  R^t)  and  0k(t)  are  the  envelope  and  phase  of  zk(t),  respec- 
tively, ©k(t)  is  the  phase  derivative  or  instantaneous  angular 
frequency  of  zk(t)  with  respect  to  v^,  and  y„  is  the  centroid  of 
the  power  spectrum  of  zk(t).  Vq  may  also  be  called  the  mean 
Doppler  shift  on  the  "carrier"  frequency  of  zk(t). 

The  Doppler  estimate  to  be  analyzed  is  repeated  below  [in 
somewhat  altered  notation  from  (2.60)]: 


K 


v = 


N 

£ £ |w  (TT)  |Zy  (TT)/2tt 

1=1  k=l  k 


K 


N 

£ £ |w  (TT) 

1=1  k=l  k 


(2.66) 
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where 


rk(t)  = £ wk(t) 


Setting  tQ  = 0 in  (2.65)  for  simplicity,  we  see  that 


j <0k  + 2ny0)fc 

„k(t)  = V ke 

= zk(t) [1  + n^Ct)] 


+ "k(t) 


(2.67) 


(2.67) 


where  the  noise  is  given  by 


"k(t)  -jek  ~i('\  + 2^0)t 


/ . \ IN.  IN. 

e 8 


(2.68) 


For  large  SNR, 


o 'nk(t)' 

indoor  = — 2 — 1 


(2.69) 


Note  that  m^t)  is  essentially  a scaled,  frequency-shifted,  ver- 
sion of  the  additive  noise  nk(t) , where  the  frequency  shift  is 


vk  V0  + 2ff 


(2.70) 


Thus,  the  power  spectrum  of  mk(t) , Pmk(f),  is  related  to  that  of 


n,  (t)  and  P (f)  by  the  relationship 
k n 
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(2.71 


Vf>  = 3 Pn<f  + V 


where  we  have  considered  0^.  and,  thus,  fixed. 

If  we  represent  m^t)  in  terms  of  its  quadrature  components 


i\(t)  = uk(t)  + jvk(t) 


then  the  high  SNR  quadrature  approximation  is  given  by 


$ [1  + n^Ct)]  » vk(t) 


so  that 


Vk(t)  “ (0k  + 2irv0)t  + Vk(t) 


and 


\(t)  , \ A Vc) 

= Vn  + TT  + 


2v 


0 2v  2it 


Using  (2.75)  in  (2.66), 


(2.72) 


(2.73) 


(2.74) 


(2.75) 


N K N K 

L E |w.  (4T)  pe.  /2ir  E E |w,  (IT) \\ (lT)/2n 


n 


We  may  further  simplify  (2.76)  by  noting  from  a high  SNR 
approximation  that 

|wk(t)  |2  - r£[1  + 2Uk(t)]  (2.77) 

It  is  not  difficult  to  see  that  the  noise  component  in  the 
squared  envelope  of  wk(t)  constitutes  a second-order  effect  in 
the  output  noise  of  the  Doppler  estimate  (2.72).  Thus,  in  (2.76) 
we  replace  |w^(t) 


2 2 
by  Rk, 


yielding 


K 2.  K / N . 

E R^/2*  E Rk  U E Vk(tT)/2» 

A , k=l  , k=l  \ 1=1  y 

y + + 


E Rk 

k=l 


K 2 
E R. 

k-1  ^ 


(2.78) 


The  first  term  on  the  right-hand  side  of  (2.78)  represents 
the  true  Doppler  and  the  second  and  third  terms  are  noise  terms. 
The  second  term  is  the  estimation  error  due  to  the  fading  char- 
acteristics of  the  channel,  while  the  third  term  represents  the 
estimation  error  due  to  the  additive  noise. 

Computation  of  the  mean  value  of  the  Doppler  estimator 
(2.78)  shows  it  to  be  unbiased.  In  actual  fact,  the  mean  value 
of  the  unapproximated  estimator  (2.66)  is  biased,  but  the  bias 

2 2 

becomes  small  at  high  SNR.  This  bias  is  given  by  - u q ( | n 1 / \z  \ + 
2 

|n  | ),  and  would  be  around  170  of  for  a 20  dB  SNR.  Because  of 
our  high  SNR  approximation,  this  bias  does  not  appear  in  (2.78). 
However,  the  variance  of  the  estimate  (2.78)  should  be  acceptable 
as  a measure  of  the  variance  of  the  estimate  (2.66)  at  high  SNR. 
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We  focus  our  attention  on  determining  the  variance  of  the 
noise  terms  in  (2.74)  conditioned  on  the  set  k = l,  . . . , K) 

being  available.  From  Rice  [2.7]  the  joint  density  function  of 
0k,  and  ©k  is  given  by 


R, 


“(VVV-o 


0 iT  BqVttPq 


exp 


£ 

P 


0 


exp 


. 2 

_e_ 

2n 


VoK 


(2.79) 


where 


po  - - he'2 


(2.80) 


B0  = 2 


J O - Vq)  P(v)  dy 


J*  P(y)  dy 


(2.81) 


J yP(y)  dy 


v0 


J P(y)  dy 


(2.82) 


and  P(y)  is  the  power  spectrum  of  zk(t) . 

From  these  equations  we  see  that  the  conditional  density 


function  of  0,  , given  R,  , is  Gaussian,  zero  mean,  and  with  variance 

2r\  r\  K K 

t-\  r*  ^ /On  ^ 


iri'PQBQ/2Rk . In  addition,  if  we  assume  that  the  tones  fluctuate 


independently,  the  0k  will  be  independent.  As  a result, 
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(2.83) 


We  consider  now  the  last  term  in  (2.23).  Note  that  the 
sampling  interval  T has  been  selected  larger  than  the  correlation 
interval  of  the  complex  Gaussian  noise  process.  As  a result,  the 
set  (V^(-tT) ; k = 1 , . . . , K;  l = 1 , . . . , N]  of  random  variables  are 
statistically  independent  and  of  common  variance.  This  variance 
is  given  by 


tyt)!2  = \ 


00 

± J (2«£)2P  (£>  <J£ 

-»  K 


*4  f f2pn(£  + V d£ 

®k  -* 


(2.84) 


where  we  have  made  use  of  (2.71). 

Assuming  the  noise  has  a symmetrical  spectrum  and  defining 
the  rms  bandwidth  of  the  noise  as 


Bn  - 2 


; f pn(f>  df 

I Vf>  df 


(2.85) 


we  see  that  the  variance  of  V^(t)  conditioned  on  and  0^  is 
given  by 
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we  see  that  the  variance  of  Vk(t)  conditioned  on  and  0^  is 
given  by 


2 2tt2  (BNPN  . 2, 


[vk(or-^ 


(2.86) 


where 


PN  = I Pn(f)  df  = ,n(t) 


(2.87) 


To  simplify  the  calculations,  we  neglect  the  term  0k  con- 
tained within  vk,  and  thus  replace  i^k  by  to  obtain  the  approx- 
imation 


2 2”2pN  l BN  . 2 


Rk 


(2.88) 


The  neglect  of  0k  is  justified  because  it  may  be  shown  to  con- 
tribute a component  to  the  estimator  variance  which  is  small 
compared  to  the  second  term  in  (2.78),  as  long  as  the  SNR  is  high. 

Using  (2.33)  and  the  previously  noted  observation  that  the 


[V,  (tT)}  are  statistically  independent  leads  to  the  result 


[Rr  » k 1 , . . . , K] 


i pnbn/8+Vo/2 


N 


K 2 
L 1L 

k=l  k 


(2.89) 
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I *1*0.1. 


Thus,  the  variance  of  the  Doppler  estimator  is  given  approx- 


imately at  high  SNR  (conditioned  on 

i 


2 

E R,)  by 
k=l 
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L R, 


where 
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+ i!s 

N P 

0 


(2.90) 


(2.91) 


For  N = l,  this  expression  should  approximate  the  variance  of  the 
single-sample  estimation  at  high  SNR  evaluated  in  [2.6].  Indeed, 
examining  Eq.  (23)  of  Appendix  B in  [2.6]  [use  Eq . (2.6)  for 
= 0] , we  find  that  they  agree  at  high  SNR. 

2 ^2 

It  remains  now  to  average  a over  £ R,  subject  to  the 

V k=l 

decision  inequality  (2.8).  Because  of  the  high  SNR  approximation, 
we  replace  the  inequality  (2.53)  by 


K 

E 

k=l 


Rk>^ 


(2.92) 


where  the  random  variable  /3  is  the  quadratic  form 


N K 

0 - jjf-  E S 1VIT>  ' 

t=l  k=l  K 


(2.93) 
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The  random  variable 


K 


y = 


0 k=l 


(2.94) 


has  the  x probability  density  and  distribution  function  with  2K 
degrees  of  freedom: 


Vy)  ‘ <K^iyr yK'le"y  1 y>0 


(2.95) 


-Y 


K-l 


Fr(Y)  = P(y  < Y)  = 1 - e L — Y 


-L  vJ 


j-o  J- 


(2.96) 


The  desired  estimation  variance  averaged  over  y,  subject  to 
the  inequality  (2.37)  but  with  0 fixed,  is  given  by 

J dy 

E(o2/y  > 0;/S)  - A 2-^ (2.97) 

J Wk(y)  dy 
0 

But , 

7 wk<y)  ■ kTT  wK-i(y)  ; K 2 2 <2-98) 


so  that 


The  final  step  in  the  evaluation  of  the  estimator  variance 

2 

is  the  average  of  (2.44)  over  j3.  We  note  that  f3  has  a X den- 
sity function  with  2KN  degrees  of  freedom: 


/ NPf 

W(S> ' 
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(2.100) 


The  mean  and  variance  of  0 are: 


F = (5) KN 


(2.101) 


°l  - ^ - (s^)  <KN><KN  + D - (s^)  kV 


1 a? 
, NP 


KN 


(2.102) 


and  the  ratio  of  standard  deviation  to  mean  is: 


(2.103) 


Since  KN  will  be  much  bigger  than  1 in  applications,  it 

does  not  appear  that  significant  errors  will  be  introduced  in  the 

2 

final  average  over  ct  if  we  replace  0 by  its  average  value,  i.e., 
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(2.104) 


It  has  been  assumed  that  the  Doppler  shift  remains  strictly 
constant  during  the  measurement  period.  In  actual  practice,  the 
Doppler  shift  may  vary  at  a rate  of  3.5  Hz/second.  The  resultant 
variation  during  the  measurement  period  will  result  in  some 
averaged  Doppler  estimate.  Clearly  it  is  desirable  to  keep  the 
measurement  time  as  small  as  possible. 

The  result  given  by  (2.104)  and  (2.91)  is  of  primary 
importance  in  evaluating  the  operation  of  the  multi-stage  pro- 
cessor described  in  the  following  section. 

2.5.2  Multi-Stage  Doppler  Estimator 

Applying  the  reasonable  high  SNR  approximation  to  (2.104) 
and  combining  with  (2.91)  we  conclude  from  the  basic  theoretical 
considerations  of  the  previous  section  that  a single-stage  esti- 
mate of  Doppler  shift  based  on  reception  of  a 4-tone  preamble 
yields  an  estimate  with  variance  given  approximately  by 
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3NP, 


B 


+ — - 
1 n / 


0 

24 


(2.105) 


where  Pq  and  P^  are  the  signal  power  and  noise  power,  respec- 
tively at  the  tone  filter  output,  is  the  rms  bandwidth  of  the 

tone  filter,  is  the  Doppler  shift  on  the  incoming  signal,  Bq 
is  the  rms  Doppler  spread  of  the  channel,  and  N is  the  number  of 
independent  time  samples  of  each  output  tone  used  in  the  Doppler 
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estimate.  Denoting  the  tone  filter  transfer  function  by  H(f) 
we  make  use  of  the  definitions 


BF  = 


J )H(f)p  df 

\H(0)  |2 


BN  = 


4 J f2'H(f)l2 
J \H(f ) \2  df 


and  write 


PN  BFN0 


po  - \ 'H^o>  '2 


(2.106) 


(2.107) 


i^.ios) 

(2.109) 


where  Nq  is  the  single-sided  noise  power  density  and  P is  the 
total  input  power  in  the  four  tones. 

2 

In  the  subsequent  discussion  we  neglect  the  term  Bq/24  which 
will  usually  be  small  compared  to  noise  contributions. 

It  is  convenient  to  relate  BX1  and  B„  and  the  filter  time 

N F 

constant  tc  to  the  half-power  frequency  of  the  filter  H(f). 

To  momentarily  avoid  considering  a particular  filter  we  write 


BF  afl/2 


= bf1/2  tc  = c/f1/2 


(2.110) 
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where  a,  b,  and  c can  be  determined  for  the  particular  filter 
under  consideration.  For  the  double- tuned  filter  discussed 
later  t will  be  the  time  it  takes  the  impulse  response  of  one 
section  to  decay  to  e ^ times  its  original  value. 


Choosing  3tc  as  a reasonable  time  interval  between  samples 
to  achieve  independence  of  noise  samples  at  the  filter  output  we 
use  N = T/(3tc)  in  (2.105)  and  perform  the  other  obvious  sub- 
stitutions to  get 


2 

a 


sT|H(i/0)  | 


(b2f  1/2  + 2 
2 V — 4 — + •'0 


(2.111) 


We  now  consider  a worst-case  situation  wherein  the  half-power 
frequency  has  been  chosen  equal  to  some  maximum  design  Doppler 
shift  value,  and  the  Doppler  shift  on  the  channel  actually 
takes  on  that  value.  We  have 


2 _ 4ac 
° sT 


T + 1 


(2.112) 


As  a particular  possibility  we  consider  a double-tuned 


filter 


H(f ) = 


( j 2 fff  + a)' 


(2.113) 


Each  of  the  cascaded  sections  of  this  filter  has  an  impulse 
response 


hQ(t)  = ae"atU(t) 


(2.114) 
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The  overall  impulse  response  for  the  double-tuned  filter  is 


h (t)  = hQ(t)  ® hQ (t) 


= a^te  atU(t) 


(2.115) 


For  the  quantity  t in  (2.110)  we  choose  the  time  constant  of  a 
single  section,  i.e.. 


t = 1/a 
c 


(2.116) 


Though,  this  quantity  is  not  the  time-constant  of  the  overall 
filter  in  the  usual  sense,  it  serves  our  purposes  here  and  is 
referred  to  as  the  time  constant  in  the  remainder  of  our  present 
discussion . 

In  view  of  our  choice  of  3tc  as  the  sampling  interval  in 
our  derivation  of  (2.111)  it  is  of  interest  to  determine  the 
extent  of  correlation  between  samples  of  noise  at  the  filter 
output  taken  3t  seconds  apart.  Since  the  noise  is  flat  across 
the  band  the  correlation  coefficient  is  given  simply  by 


P (t)  = 


j h(r  + n)h*(n)  dn 
J \h(n)  \2  dn 


(2.117) 


Substituting  (2.115)  and  evaluating  the  integrals  yields 


p(t)  = [1  + a |t | ] e 


(2.118) 
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U ....  . . 


Evaluating  this  expression  at  3t^  yields 

P(3tc)  = .2  (2.119) 

Thus,  the  choice  of  3tc  as  a sampling  interval  seems  justified 
on  the  basis  of  providing  reasonably  uncorrelated  noise  samples 
at  the  filter  output. 

Now  the  constants  in  (2.110)  can  easily  be  determined  from 
(2.106)  and  (2.107).  We  get 

a = 2.441  b = 3.108  c = .1024  (2.120) 


Hence,  using  (2.102) 


2 3.414  r2 

° ' ~rr  fi/2 


(2.121) 


In  addition  to  the  usual  inverse  dependence  on  signal-to- 
noise  ratio  s and  observation  time  T,  (2.121)  shows  the  depen- 
dence on  the  half-power  frequency  f-j/2*  Obviously,  as  f^^  i-s 
increased  more  noise  is  allowed  to  contribute  to  the  estimate 
and  the  estimate  is  degraded. 

In  order  to  compensate  for  the  f^^  dependence  described 
above  we  have  considered  a multi-stage  estimation  process  where- 
in there  is  a Doppler  correction  made  at  the  end  of  each  stage 
with  a corresponding  reduction  in  filter  bandwidth,  i.e.,  £^j ^ . 
The  Doppler  correction,  of  course,  is  based  on  the  estimate  per- 
formed in  the  preceding  stage.  To  take  into  account  the  error 
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in  this  estimate  the  half-power  frequency  of  the  current  stage 
is  chosen  to  be  equal  to  three  standard  deviations  of  the  esti- 
mate in  the  previous  stage. 


For  a three-stage  estimator  the  sequence  of  filter  config- 


urations is  shown  in  Fig.  2.14.  We  have  used  f^  to  denote  the 


i'th  stage  half-power  frequency.  Note  that  the  a's  are  calcul- 


able from  (2.121)  and  the  f^  are  chosen  according  to  the  rule 


f . = 3ct  . ... 
1 l-l 


Figure  2.14  also  illustrates  the  tone  configuration.  At 


the  first  stage  a worst-case  Doppler  shift  of  75  Hz  is  assumed; 

2 

the  tone  is  estimated  with  variance  a^.  As  part  of  our  worst- 


case  consideration  a 3cr^  error  is  assumed  for  the  stage  1 esti- 
mate. This  estimate  is  used  to  perform  the  Doppler  correction 
which  just  carries  over  the  3a  ^ error  into  stage  2.  In  stage  2 
the  half-power  frequency  of  the  filter  has  been  changed  to 


f~  = 3a..  . In  this  stage  the  tone  frequency  is  estimated  with 
z i 2 

variance  a ^ and  pushing  through  the  worst-case  choice  we  again 


assume  a three-sigma  error  (this  time  302)- 

To  summarize  our  approach  to  getting  numerical  results  for 
multi-stage  estimation  of  Doppler  shift  we  have  assumed  a multi- 
stage operation  that  uses  the  previous  stage  estimate  to  make  a 
Doppler  correction.  During  the  Doppler  correction  preceding  the 
-t'th  stage,  the  half-power  frequency  is  changed  from  f ..  to 

v"  JL 


f = 3a.  n . With  this  fact  one  can  generate  a recursion  rela- 

"V  ‘b*  1 


tion  from  (2.121)  that  can  be  manipulated  to  show  that  the  esti- 
mate variance  at  the  output  of  the  L'th  stage  is  given  by 
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Figure  2.14  Three-Stage  Adaptive  Filter  for  Doppler  Shift  Estimation 
(Order  of  Occurrence  is  Top-to-Bottom)  "Worst-Case 
Performance  is  Illustrated 
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(2.122) 


where  T is  the  processing  time  for  the  t'th  stage. 

Ay 

Recall  that  our  design  goal  is  to  minimize  total  processing 
time  to  achieve  a given  level  of  estimator  performance.  In  doing 
this  we  fix  the  estimator  variance  to  some  desired  design  value 
and  encounter  the  mathematical  problem  of  solving  the  equation 


T1T2  ' ' ' TL  = const  (2.123) 

for  the  set  of  values  [T\;  i = 1,L]  that  minimizes  the  sum 


tot 


= T1  + T2  + 


(2.124) 


This  is  achieved  by  using  a basic  theorem  of  analysis  that  states 
that  the  geometric  mean  of  the  elements  in  a set  (a. a„  ...  a.l 
is  less  than  the  arithmetic  mean  of  the  elements  unless  the  ele- 
ments are  all  equal,  i.e.,  defining 


G(  a) 


n 


7 1 

i=l 


ai 


1/n 


and 


M(a) 


\ L 
° i=l 


(2.125) 


(2.126) 
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we  have 

/ 

G(a)  < M(a) 

with  equality  holding  only  if 


(2.127) 


al  ~ a2 


n 


(2.128) 


Since  G(a)  is  constant  in  our  context  and  since  minimizing  the 
sum  in  (2.124)  above  is  equivalent  to  minimizing  the  arithmetic 
mean,  the  inequality  (2.127)  implies  directly  that  the  sum  in 
(2.124),  subject  to  the  fixed  product  constraint  in  (2.123),  is 
minimized  by  choosing  all  the  T\  equal. 

Thus,  each  stage  processes  for  equal  time  intervals  of  dur- 
ation T.  From  (2.122)  this  value  is  given  by 


T - 9 ( 


/ f 

/ 3 . 414  \ 

1 

Is) 

|U/  J 

1/L 


(2.129) 


Sampling  rates,  of  course,  will  reduce  appreciably  as  we 
move  to  later  stages.  Using  the  fact  that  the  sampling  interval* 
in  each  stage  is  assumed  to  be  three  time  constants  we  have  from 
(2.110)  and  (2.120) 


Tf . 


3072 


(2.130) 


where 


h ' 30 <.-i 


(2.131) 


* 


By  sampling  interval  we  here  refer  to  the  time  between  samples 
used  in  the  computation.  This  would  be  implemented  as  an  integer 
multiple  of  the  fundamental  sampling  interval. 
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From  the  expressions  provided  above,  the  data  of  Table  2-5 
has  been  compiled  using  an  s value  (signal-to-noise  density  ratio) 
of  40  dB/Hz.  In  this  table  T represents  the  duration  of  each  stage 
as  determined  from  (2.129);  T is  just  this  quantity  multiplied 
by  the  number  of  stages.  The  number  of  samples  taken  during  each 
stage  is  determined  from  round-off  of  (2.130).  (For  those  stages 
in  which  round-off  has  resulted  in  a zero  value,  the  value  has 
been  replaced  by  unity  and  the  stage  number  marked  with  an  asterisk.) 

All  of  the  quantities  to  the  left  of  the  double  line  in 
Table  2-5  are  directly  determinable  from  the  expressions  pre- 
sented above.  To  the  right  of  the  double  line  two  corrections 
are  represented: 

a)  Integer  sample  correction 

b)  Transient  build-up  time  correction. 

The  first  correction  derives  from  the  fact  that  our  solution  for 
T does  not  necessarily  represent  an  integer  number  of  samples; 
to  correct  this  we  have  adjusted  T to  an  appropriate  value  T^ 
for  each  stage,  T^  being  the  sampling  period  for  that  stage  (three 
time  constants)  multiplied  by  the  integer  determined  from  round- 
off of  (2.130).  With  this  correction  we  have  calculated  [using 
(2.122)  now]  the  new  estimator  standard  deviation  and  total  pro- 
cessing time.  These  are  shown  in  the  first  two  columns  to  the 
right  of  the  double  line  in  Table  2-5.  In  the  third  column  we 
show  the  processing  time  that  is  obtained  when  the  transient 
build-up  time  between  stages  is  taken  into  account.  This  is 
obtained  by  adding  to  T t one  sampling  interval  (three  time  con- 
stants) for  each  stage  after  the  first. 

Among  other  things,  we  see  from  this  table  that  a three- 
stage  processor  can  be  used  to  reduce  estimation  error  to  sigma 
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values  of  1,  .75,  .5,  and  .25  Hz  in  .1305,  .157,  .2125,  and  .3181 
second,  respectively. 


Because  of  channel  Doppler  spread  (largest  typical  value 
approximately  1 Hz)  it  seems  illusory  to  attempt  estimates  of 
much  higher  precision  than  the  cases  cited  above. 

Some  indication  of  the  dependence  of  Doppler  estimator  per- 
formance on  SNDR  is  given  by  Fig.  2.15.  It  can  be  seen  that 
there  is  no  substantial  variation  of  estimator  performance  within 
a few  dB-Hz  (+  2)  of  40  dB-Hz. 

2 . 6 Frame  Sync  Estimation 

In  this  section  we  discuss  the  segment  of  the  preamble 
utilized  for  frame  synchronization.  The  structure  of  this  seg- 
ment, referred  to  as  the  sync  preamble,  will  be  described  in 
detail  in  Section  2.6.1  below.  Integration  of  the  sync  preamble 
into  the  preamble  as  a whole  was  discussed  in  Section  2. 

The  performance  of  the  frame  sync  estimator  can  be  degraded 
by  several  factors: 

1)  Doppler  correction  error 

2)  Terminal  equipment  distortion 

3)  Poor  threshold  setting 

4)  Additive  noise 

5)  Multipath. 

2.6.1  Basic  Elements  of  Sync  Preamble 

The  overall  structure  and  timing  of  the  sync  preamble  was 
described  in  Section  2.  The  basic  signal  element  of  the  sync 
preamble  was  seen  to  be  a pulse  train  designed  both  for  high  crest 
factor  and  low  autocorrelation  side  lobes.  It  is  assumed  that  in 
the  receiver  Doppler  correction  is  carried  out  prior  to  receipt 


Figure  2.15  Performance  of  Doppler  Estimator  for  Design 
Value  of  a = . 5 Hz  (a  s RMS  Doppler  Shift 
Estimation  Error) 


of  the  sync  preamble,  and  that  the  latter  preamble  duration  is 
small  compared  to  the  reciprocal  Doppler  error,  so  that  ambiguity 
function  design  of  the  sync  preamble  is  unnecessary.  To  test 
the  validity  of  this  assumption  the  effect  of  Doppler  error  on 
the  sync  estimation  is  evaluated  in  Section  2.6.3. 

The  pulse  train  is  assumed  generated  by  a sequence  of  two 
operations.  First,  a discrete  sequence  is  generated  which  has  a 
desirable  discrete  autocorrelation  function  (i.e.,  low  side  lobes) 
This  discrete  sequence  is  fed  to  a pulse-shaping  filter  which  is 
designed  to  have  the  following  properties: 

1)  The  spectrum  of  the  pulse  is  kept  within  the 
"flat"  passband  of  the  channel  terminal  equip- 
ment to  prevent  ringing  or  pulse  broadening. 

2)  The  pulse  shape  is  designed  to  preserve,  as 
much  as  possible,  the  good  crest  factor  of  the 
discrete  sequence,  and  to  provide  low  side 
lobes  at  the  receiver  matched  filter  output. 

In  practice,  the  pulse  shaping  filter  will  include  any  zero- 
order  holds  or  anti-alias  filters  present  in  the  D/A  conversion 
process.  The  major  part  of  the  pulse  shaping  can  be  done  in  the 
computer  or  by  an  external  analog  filter. 

We  consider  first  possible  discrete  sequences  with  good 
autocorrelation  properties.  Selection  of  the  discrete  sequence 
is  influenced  in  part  by  the  anticipated  receiver  processing. 

To  simplify  the  latter,  the  discrete  sequence  should  ideally 
contain  numbers  that  allow  the  matched  filtering  in  the  receiver 
to  use  few  multiplies.  A sequence  containing  only  the  numbers 
± 1 would  avoid  the  use  of  a multiplier  entirely.  Use  of  ± j 
weights  should  also  avoid  the  use  of  a multiplier  since  multiplica 
tion  by  ± j involves  interchanging  the  real  and  imaginary  parts  of 
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words  with  possible  sign  changes.  The  inclusion  of  multiples 
of  45°  in  the  phase  weights  would  introduce  the  requirement  to 
have  available  the  input  in-phase  and  quadrature  components 
weighted  by  ± 1/72  . 


If  weights  of  ± 1 and  ± j are  used,  only  relatively  short 
sequences  are  known  which  have  low  side  lobes.  From  the  paper 
by  Frank  [2.3],  one  may  find  a 64  sample  sequence  containing 


1 i 

+ 1 , + ] , + — + j — weightings  with  very  low  side  lobes  in 
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its  autocorrelation  function.  This  sequence  is  the  entry  N = 8 
in  Table  I of  [2.3].  The  ratio  of  autocorrelation  peak  to  side 
lobe  for  this  waveform  is  25.  A sequence  of  length-15  is  given 
by  Carley  [2.9],  containing  ± 1,  ± j weightings,  which  has  an 
autocorrelation  function  peak-to-sidelobe  ratio  of  15.  The 
length-64  and  -15  waveforms  are  listed  here  in  Fig.  2.16  together 
with  the  longest  Barker  sequence  known  [2.10]  (length-13)  con- 
taining ± l's  as  elements.  This  latter  sequence  has  an  auto- 
correlation function  peak-to-sidelobe  ratio  of  13. 


The  transmitted  pulse  train  complex  envelope  z(t),  ^ay  be 
represented  as  the  convolution  of  an  impulse  train  having  areas 
corresponding  to  the  selected  discrete  sequence  with  a continuous 
filter  p(t)  providing  the  pulse  shape,  i.e., 


M-l 

z(t)  = p(t)  ® E y 6(t  - mA) 

0 m 

M-l 

= E y p(t  - mA)  (2.132) 

^ m 


where  (y^ ; m = 0,1,  ...,  M-l]  is  the  discrete  sequence,  6(>)  is 
the  unit  impulse,  and  A is  the  time  between  samples. 
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LENGTH- 13 
SEQUENCE 


Figure  2.16  Sequences  of  Length-13,  -15,  and  -64  with  Low  Autocorrelation  Sidelobes 


Since  the  have  unit  magnitude,  the  discrete  sequence  may 
be  said  to  have  a unity  crest  factor.  The  pulse  shape  p(t)  is 
assumed  to  be  of  the  Nyquist  type,  having  the  normalized  value 
unity  as  t = 0 and  zero  value  at  time  samples  equal  to  multiples 
of  A.  Such  a pulse  shape  assures  that  z(t)  has  constant  magnitude 
at  the  time  instants  t = mA,  m an  integer.  The  values  inbetween 
will  generally  not  be  constant  envelope.  However,  if  p(t)  is 
chosen  to  drop  to  zero  rapidly  for  t outside  the  interval  -A<  t < A, 
the  values  inbetween  will  be  primarily  some  linear  combination  of 
adjacent  values,  and  thus  quite  limited  in  their  fluctuations. 
Another  requirement  on  p(t)  is  that  it  be  bandlimited  since  the 
spectrum  of  the  pulse  must  lie  within  the  flat  portion  of  the  HF 
terminal  equipment  channel.  A raised  cosine  spectrum  of  band- 
width 2/A  satisfies  the  Nyquist  requirement.  The  low  pass  spec- 
trum that  would  be  used  is  shown  in  Fig.  2.17.  The  shaded  regions 
are  frequency  guard  bands  placed  where  the  transceiver  passband 
is  not  flat.  If  a 2 kHz  wide  flat  passband  is  available,  one  may 
select  A = 1 ms.  A length-64  discrete  sequence  would  occupy  64  ms 
and  the  analog  pulse  train  an  additional  few  milliseconds,  or 
around  three  frames  for  a typical  22.5  ms  frame  duration. 

The  raised  cosine  spectrum  is  only  given  as  an  illustration, 
and  other  shapes  are  possible.  It  must  be  realized  that  the  so- 
called  "flat"  part  of  the  transceiver  will  have  some  "bumps"  and 
produce  pulse  shape  deterioration.  The  following  section  will 
include  a discussion  of  the  interaction  between  receiver  process- 
ing and  pulse  shape  in  producing  low  side  lobes  at  the  matched 
filter  output. 


Signal  processing  of  the  synch  preamble  is  illustrated  in 
Fig.  2.18.  Estimation  of  frame  sync  takes  place,  of  course, 
after  signal  detection,  Doppler  estimation,  and  Doppler  correc- 
tion. The  estimation  of  frame  sync  is  carried  out  by  utiliza- 
tion of  a "matched"  filter  followed  by  a multipath  centroid 
detector  which  locates  all  peaks  above  a given  threshold  and 
forms  the  centroid  of  their  locations. 

As  indicated  in  Fig.  2.18,  the  general  structure  of  the 
matched  filter  would  be  the  cascade  of  a pulse  matched  filter 
matched  to  the  transmitted  pulse  shape  and  a discrete  FIR  filter 
matched  to  the  M-length  discrete  sequence  formed  in  the  trans- 
mitter. Referring  to  (2.1),  m = 0,1,  ...,  M-l]  denotes  the 

discrete  sequence  and  p(t)  denotes  the  transmitted  pulse  yield- 
ing a transmitted  signal  with  complex  envelope  equal  to  z(t)  in 
(2.1). 

Since  the  received  pulse  filter  may  not  be  matched  exactly 
to  the  transmitted  pulse,  we  denote  the  pulse  filter  impulse 
response  by  q*(-t),  and  the  convolution  of  p(t),  q*(-t)  becomes 

r(t)  = J q*(r)p(r  + t)  dr  (2.133) 

When  the  receive  pulse  filter  is  truly  matched,  q(t)  = p(t) 
yielding 

r (t)  = J*  p*(t)P(t  + t)  dr  (2.134) 

At  the  receive  pulse  filter  output,  the  pulse  train  is  given 
by 


Figure  2.18  Sync  Preamble  Processing 


I 


I 

I I 


w^(t)  = z(t)  ® q*(-t) 


M-l 


= £ y r (t  - mA) 

_A  ^ 

m=0 


The  discrete  matched  filter  has  the  impulse  response 


£ y*  6(t  - mA) 


(2.136) 


(apart  from  an  unimportant  delay)  which  yields  at  its  output  the 
pulse  train 


w(t)  = £ Cmr(t  " mA) 


(2.137) 


where  C is  the  discrete  autocorrelation  of  the  discrete  sequence 
m 

tv 


cm  ■ x 


(2.138) 


As  discussed  above,  [y^]  is  selected  so  that  [C^]  has  low  side 
lobes.  If  q*(-t)  is  selected  so  that  r(t)  is  of  the  Nyquist  form 
with  zeroes  at  t = mA  and  normalized  unit  value  at  t = 0,  then  the 
good  side  lobes  in  will  be  preserved  in  w(t)  at  the  sampling 
instants  t = mA,  i.e., 

w(tA)  = Ct  (2.139) 

Moreover,  if  r(t)  is  small  outside  the  interval  -A  < t < A,  there 
will  be  limited  fluctuation  between  sample  points. 
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To  make  the  complexity  of  the  receiver  pulse  filter  as 
small  as  possible,  a recursive  filter  realization  should  be  * 

used.  However,  such  a realization  conflicts  with  the  previous 
specification  that  the  transmitted  pulse  be  bandlimited.  An 
approximation  is  clearly  necessary.  One  may  also  entertain  the 
idea  of  removing  the  receive  pulse  filter  entirely  with  probably 
only  a couple  of  dB  loss  in  SNR  performance  for  a considerably 
simpler  receiver  processing. 

With  proper  filter  and  sequence  design  and  a close  enough 
Doppler  correction,  the  matched  filter  output  from  an  ideal 
single  path  channel  will  consist  of  a single  narrow  pulse  plus 
low-level  side  lobe  hash.  Multiple  paths  will  yield  multiple 
pulses  with  low-level  side  lobe  hash.  Because  of  the  previous 
signal  detection  and  signal  and  noise  energy  measurement  opera- 
tions and  the  known  relation  between  the  transmitted  Doppler 
preamble  and  sync  preamble  strengths,  it  should  be  possible  to 
set  a threshold  well  above  the  rms  noise  level  (and  certainly 
above  the  hash  level)  so  that  resolvable  path  components  will  be 
deary  visible.  Algorithms  can  be  developed  for  finding  the 
location  and  strength  of  maxima,  (£  ,a  },  above  the  threshold, 
and  computing  the  centroid 


i = E a i /E  a (2.140) 

s man  m 

From  this  centroid  one  may  determine  the  frame  sync  to  within  an 
error  determined  primarily  by  the  multipath  spread  but  also  by 
the  width  of  r(t)  and  the  SNR. 

There  are  clearly  certain  pitfalls  in  the  above  operations. 
In  determining  the  pulse  threshold,  one  must  consider  the  action 
of  the  AGC  when  the  Doppler  preamble  is  followed  by  the  sync 
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preamble,  because  this  procession  may  produce  gain  changes.  In 
selecting  maxima,  one  must  consider  that  two  or  more  paths,  which 
are  close  together  and  not  resolved , may  produce  subsidiary  maxima 
that  would  confuse  the  computation.  In  fact,  to  simplify  the  pro- 
cessing, it  may  be  sufficient  to  use  only  the  largest  peak  in 
determining  frame  sync.  This  estimation  procedure  should  produce  a 
timing  error  within  the  range  of  protection  of  the  time  guard  bands. 

2.6.2  Performance  of  Sync  Preamble 

In  this  section  we  discuss  some  sources  of  degradation  in 
sync  preamble  operation.  These  are  seen  to  be  Doppler  correction 
error  (Section  2. 6. 2.1),  poor  threshold  setting  and  additive  noise 
(Section  2. 6. 2. 2).  Terminal  equipment  distortion  and  effects  of 
multipath  sidelobes  are  not  treated  here.  These  latter  disturb- 
ances will  require  an  increase  in  the  threshold  setting  to  achieve 
the  required  false  peak  detection  probability  in  comparison  with 
the  absence  of  these  disturbances. 

2. 6. 2.1  Effect  of  Doppler  Correction  Error 

A sync  preamble  designed  using  the  simple  criterion  that  the 
set  m = 0,1,  ...,  M-l]  should  be  chosen  to  provide  desired 

correlation  properties  discussed  in  connection  with  (2.137)  may 
not  provide  suitable  protection  against  Doppler.  Because  of  the 
inevitable  (but  predictably  small)  Doppler  shift  remaining  after 
Doppler  correction,  sequences  of  one  length  may  be  preferable 
over  sequences  of  another  for  reasons  that  do  not  derive  strictly 
from  correlation  considerations. 

To  analyze  this  situation  we  assume  that  a Doppler  correc- 
tion error  denoted  F exists.  In  such  a case  the  complex  repre- 
sentation of  the  received  sync  preamble  for  a unity  gain  channel 
is  given  by 


L 
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(2.141) 


z(t)  - E ymP(t  - mA) 

m=0 

This  may  be  rewritten  in  the  altered  form: 


M-l 

- v j2irmFA  j2TTF(t-mA)  /0  ,/0x 

z(t)  = E y eJ  p(t  - mA)eJ  (2.142) 

m=0  m 

which  may  be  interpreted  as  the  transmission  of  an  impulse  train 
with  complex  amplitudes  exp  (j2irmFA)]  through  a pulse-shaping 

filter  with  impulse  response  p(t)  exp  (j2frFt).  The  dominant 
effect  concerns  the  "mismatching"  of  the  impulse  complex  ampli- 
tudes from  their  correct  values  [v  ] . The  slight  frequency  shift 
of  the  pulse  filter  produced  by  the  factor  exp  (j2rrFt)  multiply- 
ing p(t)  will  produce  negligible  mismatch. 

The  receiver  may  be  regarded  as  the  cascade  of  a pulse 
matched  filter  and  a discrete  filter  matched  to  the  transmitted 
discrete  sequence.  Using  the  same  notation  as  in  the  previous 
section,  we  see  that  the  receive  pulse  filter  output  is  given  by: 

w1(t)  = z(t)  ® q*(-t) 

= E y ej2fffmAr(t -mA,F)  (2.143) 

m=0 

where 


r(t,F)  = J*  p(r)q*(T  - t)ej2nlT  dr  (2.144) 

But  since  p(t)  is  of  duration  very  much  less  than  the  product 

r 

2ttFt  will  be  very  small,  and  we  can  use 
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r(t,F)  - r(t,0)  = r(t) 


(2.145) 


Thus , 


w- (t)  E y e^ff^m^r(t  - mA) 
m=0  m 


(2.146) 


and  the  discrete  matched  filter  has  the  output 


M-l 

w(t)  = E A (F)r(t  - mA , F ) 
m=- (M-l) 


(2.147) 


where 


'm-i 
,E  Vt- me 

t,=m 


j2irF4A 


Am(p)H 


M-l+m 

L y . y*  e 

{1=0  1 l~m 


j 2trF't  A 


0 < m < M-l 


(M-l  < m < 0 


(2.148) 


is  the  ambiguity  function  of  the  discrete  sequence  (y^;  m = 0,  M-l}. 
When  the  Doppler  error  F is  zero,  this  expression  reduces  to  the 
autocorrelation  function  of  the  discrete  sequence,  i.e., 


Am<°>  " Cm  (2.149) 

* 

If  F is  too  large,  the  side  lobes  of  w(t)  will  be  degraded;  and 
if  F is  sufficiently  large,  even  the  main  pulse  lobe  will  be 
reduced.  To  determine  sensitivity  to  Doppler  correction  error, 
Am(F)  must  be  evaluated  for  any  candidate  pulse  train.  A con- 
servative requirement  for  small  degradation  is  that  2ntFL  be  much 
less  than  unity  for  the  maximum  value  of  l = M-l,  i.e., 
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2w(M  - 1 ) FA  « 1 


(2.150) 


then  the  exponential  in  (2.148)  will  differ  little  from  unity. 

For  A = 1 ms,  M-l  - 63,  we  see  that  (2.150)  leads  to 

F « 2.53  Hz  (2.151) 

which  is  a stringent  specification  on  Doppler  correction  error. 

When  the  channels  contain  L paths  with  relative  delays  £ , 
relative  Doppler  shifts  v and  amplitudes  g^,  the  matched  fil- 
ter output  becomes 

L 

w(t)  = E g,  S Ak(F  + i^)r(t  - kA  - 4k)  (2.152) 

1=1  ' k 

Note  that  the  relative  Doppler  shifts  between  multipath  compon- 
ents can  result  in  output  pulse  degradation  due  to  the  ambiguity 
function  being  evaluated  at  F + u . This  degradation  persists 
even  if  F = 0,  and,  for  a given  waveform,  can  only  be  reduced  by 
shortening  the  duration  of  the  sync  preamble.  An  alternate  route 
is  to  use  Doppler-shift  tolerant  discrete  sequences.  The  develop-  ■ 
ment  of  such  sequences  requires  the  use  of  computer-aided  numer- 
ical design  techniques,  and  is  not  attempted  here  because  of  time 
limitations . 

To  pin  down  the  effect  of  the  Doppler  correction  error  on 
frame  sync  estimation  we  have  calculated  values  of  the  ambiguity 
function  for  three  sequences  of  interest: 
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1) 

Length- 13  — Barker 

[2.10 

2) 

Length-15  — Carley 

[2.9] 

3) 

Length- 64  — Frank 

[2.31  • 

Results  are  illustrated  in  Figs.  2.19  through  2.21  for  F in 
the  range  0 to  2 Hz.  For  each  sequence  we  note  from  the  F = 0 
entry  that  the  desired  autocorrelation  properties  (low  side 
lobes)  are  present.  For  the  length-13,  length-15,  and  length-64 
sequences  the  values  of  main-to-side  lobe  ratio  p are  13,  15,  and 
24.5,  respectively  in  the  absence  of  Doppler  offset.  Clearly, 
with  F = 0 the  third  choice  would  be  the  best  of  the  three  were 
there  no  residual  Doppler  shift. 

In  view  of  the  achievable  precision  in  estimating  Doppler, 
as  described  earlier,  and  the  Doppler  spread  on  HF  channels  (a 
reasonable  maximum  would  be  about  1 Hz)  it  seems  reasonable  to 
design  for  F = 1 Hz.  In  this  case  the  length-13,  length-15,  and 
length-64  codes  have  main-to-side  lobe  ratios  of  13,  14.99,  and 
11.97,  the  length-15  code  being  clearly  the  most  suitable  choice. 

A further  examination  of  the  figures  indicates  that  a design 
choice  based  on  F = 1 Hz  is  still  the  leading  contender  at  F = 2 
Hz  and  it  would  seem  that  the  crossover  point  (where  the  length- 
13  code  begins  to  surpass  the  length-15  code)  is  far  beyond  the 
range  of  Doppler  offsets  that  need  to  be  considered. 

2. 6. 2. 2 Effect  of  Additive  Noise  and  Threshold  Setting 

Before  directly  addressing  the  problem  of  the  effect  of  addi- 
tive noise  we  develop  a necessary  formalism  for  the  received  sig- 
nal. Recall  from  (2.141)  that  the  transmitted  frame  sync  pre- 
amble is  a pulse  train  with  complex  envelope 

z(t)  = E ymp(t  “ mA)  (2.153) 
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The  channel  is  modeled  by  a small  set  of  discrete  paths  having 
amplitudes  [gp;  p = 1,2,  ...,  P]  and  delays  {£  ; p = 1,2,  ...,  P] . 
Consequently,  the  received  signal  is  given  by 

w(t)  = E gpz(t  - £p)  + n(t)  (2.154) 

where  n(t)  is  an  assumed  additive  white  noise. 

An  ideal  matched  filter  receiver  is  assumed.  The  output  of 
the  matched  filter  is  thus  given  by 

u(t)  = J z*(t  - t)w(r)  dr 


= L gpC(t  - 4p)  + N(t)  (2.155) 

where  C(r)  is  the  aperiodic  autocorrelation  function  of  z(t), 

C (t)  = J z*(t)z(t  + t)  dt  (2.156) 

and  N(t)  is  the  additive  noise  at  the  matched  filter  output, 

N(t)  = J z*(t  - t)n(r)  dT  (2.157) 

The  autocorrelation  function  of  n(t)  is  given  by 

n*(t)n(t  +t)  = 2Nq6(t)  (2.158) 


where  Nq  is  the  one-sided  power  density  of  the  real  noise.  Using 
(2.158),  we  readily  find  that 


lN(t)|2  = 2Nq  J |z(t)l2  dt  (2.159) 


2-90 


J 


The  preamble  autocorrelation  function  has  been  designed  to 
have  low  side  lobes,  as  discussed  earlier.  Using  the  15  pulse 
length  Carley  sequence,  the  side  lobe  level  will  be  1/15  the 
level  of  the  peak.  This  peak  value  is  given  by 

C(0)  = J | z ( t ) |2  dt  (2.160) 

We  now  turn  to  a direct  consideration  of  the  additive  noise. 
Additive  noise  can  have  two  effects;  it  can  randomly  move  the 
location  of  the  correlation  peaks,  and  it  can  produce  new  peaks 
exceeding  the  threshold.  We  consider  each  of  these  effects  in 
order . 

To  determine  the  effect  of  noise  on  movement  of  the  corre- 
lation peak,  it  is  sufficient  first  to  consider  a single  path 
with  zero  delay  and  complex  amplitude  g.  Then  the  received  sig- 
nal is 


u(t)  = gC (t)  + N(t) 


(2.161) 


In  the  vicinity  of  the  tone  peak  location  in  the  absence 
of  noise,  t = 0,  we  can  represent  C(t)and  N(t)  by  Taylor  series: 


C(t)  « C(0) 


i + C£0JL  t + C£01 
C(0)  C(0)  2 


(2.162) 


N(t)  » N(0)  + N(0)t 


(2.163) 


where  we  have  retained  terms  to  second  order  only  (the  noise  is 
assumed  small  compared  to  signal) . The  envelope  of  the  received 
pulse  is  given  by 
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ln(t)  1 = |gC(t)  + N(t) I 


A 

« !g|c(0) 

i + N(°)  + 

[c<0,  + M 

g 

t . C(0)  t2 

1 + gC(0)  + 

C(0) 

+ C(0)  2 

- 

8 


|C(0)  [ 


1 + Re  l } + Re 


gC(O) 


N(0)  L , C£01  tl 

gC(O)  + C(0)  2 


(2.164) 


where  we  have  assumed  the  perturbation  terms  are  small  compared 
to  unity  and  have  made  use  of  the  fact  that  C(0)  is  imaginary  and 
C(0)  is  real.  In  particular. 


'cm 2 1 f df 


C(0) 


(2ir)‘ 


j iz(f)r  df 


(2.165) 


where  Z(f)  is  the  preamble  spectrum.  Defining  the  rms  bandwidth 
of  z(t)  as 


W = 2 


J 


'j  f 2 \Z  ( f ) | 2 df 


J \Z(f)  I df 


(2.166) 


we  see  that 


cm 

C(0) 


2,  ,2 

-7T  W 


(2.167) 


Completing  the  square  in  the  last  approximation  in  (2.164), 


jx(t)  ~ \g \C(0)  1 + Re 


N(0) 

gC(0) 


+ c 


o 2,  ,2 

2 it  W 


2.  ,2  , 

v W N 2 

o (t  - c) 


(2.168) 


where 


e = 


1 Re>-^ 


2 2 
n w 


gc(0) 


(2.169) 


In  the  rewritten  form  (2.168),  it  is  seen  by  inspection 
that  the  maximum  has  been  shifted  from  t = 0 to  t = c . 


The  rms  value  of  this  shift  will  now  be  computed: 


.2  = 1 l 'N(O) i 

4.,4  2 , ,2_2,n. 
v W \g  \ C (0) 


(2.170) 


where  use  has  been  made  of  the  following  relation  for  a complex 
Gaussian  variable  N: 

Re2(N]  = \ IN\2  (2.171) 


From  (2.157)  and  (2.158)  we  see  that 


Utilizing  Parseval's  theorem,  we  find  that 


|N(0) l2  = 2N  (2n)2  J f2\Z(f) l2  df 


= 2N0ir2W2C(0) 


(2.173) 


Substituting  this  expression  in  (2.170)  gives 


2 

c 


N, 


0 


ir2W2\gl2C(0) 


(2.174) 


The  received  signal  energy  over  the  single  path  is  given  by 


= \ \S \2  J |z(t)  l2  dt 


= \ \g \2C(0) 


(2.175) 


In  terms  of  this  energy, 


€ 


2,  ,2  25 

* w tT 
0 


(2.176) 


Defining  an  average  power  for  the  received  signal, 


p = — 

T 


(2.177) 


where  T is  the  duration  of  the  preamble, 


1 

J2  irW 


(2.178) 


0 

2-94 


To  be  specific,  we  assume  a raised  cosine  spectrum  for 
the  preamble: 


Z ( f ) = cos2  (nf/B)  ; \£\  < -| 


The  rms  bandwidth  W for  this  spectrum  is  given  by 


4 J f cos  (irf/B)  df 

w2  = - ■-*/* 

w R/2 

cos^  (irf/B)  df 

-B/2 


*/2  2 4 
f y:  cos  x dx 

4^  o 

2 tt/2 


J cos  x dx 


Upon  evaluating  the  integrals  in  (2.180),  we  find 


W2  (.197) 


J2  ffW  = 1.26  B 


so  that 


f?.T' 


_1 

26  B 
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(2.179) 


(2.180) 


(2.181) 


(2.182) 


(2.183) 


The  dependence  of  this  expression  on  P/Nq  is  illustrated  in 
Fig.  2.22  over  a range  of  interest.  Note,  at  40  dB-Hz  that  the  rms 
jitter  of  the  peak  is  only  .03  ms,  which  is  negligible.  It  may 
thus  be  concluded  that  the  jitter  of  the  correlation  peaks  caused 
by  additive  noise  in  the  ranges  of  interest  is  small  enough  to  be 
of  second-order  effect  since  the  guardband  available  will  be  of 
the  order  of  4 ms  and  the  multipath  spread  (separation  between 
correlation  peaks  for  all  paths)  can  be  a substantial  fraction 
of  4 ms. 

We  now  turn  to  consideration  of  the  spurious  threshold 
exceedances  that  may  occur  due  to  noise. 

There  exists  a finite  probability  that  noise  peaks  will 
exceed  the  threshold  used  to  define  correlation  peaks.  It  is  of 
interest  to  determine  the  height  of  the  threshold  required  above 
the  additive  noise  in  order  to  make  the  probability  very  small 
that  one  noise  peak  will  exceed  the  threshold  in  the  vicinity  of 
the  expected  a priori  region  of  correlation  peak  arrivals. 

To  make  the  subsequent  calculations  simple,  we  shall  assume 
that  the  spectrum  of  the  preamble  is  flat  over  B Hz  instead  of 
being  a raised  cosine  shape.  This  will  result  in  a probability 
of  threshold  exceedance  somewhat  larger  than  would  be  true  for 
the  raised  cosine  shape,  which  is  a conservative  approach  to 
system  design. 

We  assume  that  the  matched  filter  output  (complex  envelope) 
is  sampled  at  a rate  of  B samples  per  second  (a  somewhat  higher 
sampling  rate  would  probably  be  used  in  the  actual  system,  but 
this  will  not  appreciably  change  our  results).  An  uncertainty 
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region  of  duration  Tq  is  assumed  to  exist  wherein  the  correla- 
tion peaks  from  the  various  paths  would  be  expected  to  arrive. 

The  size  of  this  uncertainty  region  has  impact  on  the  sync  pre- 
amble's vulnerability  to  noise;  the  tendency  is  toward  degrada- 
tion of  performance  as  the  region  is  widened.  This  was  one  of 
the  major  motivating  factors  for  the  particular  preamble  timing 
discussed  in  Section  2.2.  We  recall  from  arguments  given  there 
that  the  actual  observation  (search  for  the  correlation  peak) 
is  made  over  x i seconds  (also  see  Figs.  2.2  and  2.5). 

From  the  specifications  of  that  earlier  section  (K^  = 45  and 
A = 1 ms)  we  will  have  occasion  to  later  use  the  specific  value 
Tq  = 45  ms.  For  the  time  being,  however,  we  will  retain  Tq  as  a 
general  unspecified  variable. 

With  a flat  spectrum  of  duration  B for  the  matched  filter, 
the  output  noise  will  also  have  a flat  spectrum  of  duration  B 
and  samples  taken  at  a rate  of  B/second  will  be  statistically 
independent.  In  a duration  Tq  there  will  be 

M = BTq  (2.184) 

independent  noise  samples. 

We  are  interested  in  the  probability  p that  one  or  more 
noise  instantaneous  power  samples  exceed  a threshold.  Along 
these  lines  we  define 

y = \N(t)  \2 

T = y = |N(t) |2 

= 2Nq  J \z(t)\2  dt  (2.185) 
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; 


the  last  equality  following  from  (2.159).  We  are  interested  in 
the  probability  that  one  sample  or  more  of  y/T  exceeds  a fixed 
threshold  p^.  This  probability  must  satisfy  the  inequality 
p < Pg,  where 

M / y(t.) 

po  ■ s.  *r[— r“>po 

m=l  ' 


- MPr(r>po) 


(2.186) 


The  bound  is  very  tight  as  long  as  p^  1 which  may  be  assumed 
here. 

The  value  of  p^  is  easily  determined  by  noting  that  the 
assumption  of  complex  Gaussian  statistics  for  the  noise  implies 
that  y is  exponential ly  distributed,  i.e.,  its  pdf  is 


P(y)  = j:e  1 
which  yields 


(2.187) 


(2.188) 


A numerical  example  of  interest  here  utilizes  B = 2 kHz, 

Pq  = 15,  and  T^  = 45  ms,  the  values  currently  under  consideration 
for  the  preamble.  From  (2.184)  we  have  M = 90.  Thus, 
pQ  = 2.75  x 10'5. 
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The  above  behavior  is  not  surprising;  as  one  increases  the 
threshold,  the  probability  of  threshold  exceedance  by  noise  alone 
decreases  exponentially.  The  threshold  cannot  be  placed  too  high, 
however,  because  the  correlation  peaks  of  the  paths,  at  least  the 
strongest  path,  must  exceed  the  threshold.  Recall  that  the 
decision  variable  in  the  absence  of  sync  preamble  is  simply  the 
normalized  noise  variable  y/T.  It  is  of  interest  to  determine 
the  value  of  the  decision  variable  at  the  correlation  peak  in  the 
absence  of  noise.  Along  these  lines  we  note  that  the  squared  cor- 
relation peak  for  the  p'th  path  is  given  by 

2 

Sp  = \gp \2  J \z(t)  l2  dt  (2.189) 

so  that 


S 

_el 

N 


'gpl2  J 'z(t) j2  dt  _ 3^  _ 


2N 


0 


N 


0 


(2.190) 


where  <?  is  the  preamble  energy  received  over  the  p'th  path,  and 

Pp  is  the  corresponding  average  preamble  power.  Assuming  T = 15  ms 

and  P /Nn  = 40  dB-Hz,  S /N  = 150.  The  peak  power  in  the  side  lobes 
P 0 p 2 

of  the  output  pulse  will  be  down  by  a factor  of  15  = 225,  and 

thus  well  in  the  ambient  noise.  For  the  example  cited  in  setting 

the  noise  threshold,  the  correlation  peak  will  exceed  the  threshold 

by  150  - 15  = 135,  which  is  large  enough  to  make  the  probability 

of  noise  cancelling  the  peak  negligible.  Reducing  the  value  of 

P^/Nq  by  a factor  of  5 (to  around  33  dB)  results  in  a correlation 

peak/noise  power  ratio  of  30  and  a threshold  exceedance  of  15,  which 

is  still  large  enough  to  make  the  probability  of  noise  cancelling 

a signal  pulse  negligible. 
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The  considerations  described  above  indicate  that  the  effects 
of  noise  (jitter  of  the  correlation  peak  and  decision  error) 
will  be  minimal.  Additionally,  they  have  presented  us  with  some 
reasonable  preliminary  design  specifications. 

2 ~ 

1)  Formulate  decision  variables  \S ( t ^ ) \ /T 

i = 1,  M where  S(t)  represents  the  sync  pre- 
amble received  in  noise  and  r is  an  estimate 
of  the  reference  noise  level  T given  in  (2.185). 

2)  Test  each  of  the  M samples  against  a thresh- 
old Pq.  As  an  initial  choice  we  use  the 
specific  value  pQ  = 15. 

There  are  three  questions  relative  to  the  above  comments  that 
we  have  not  been  able  to  pursue  in  depth  because  of  time  limitations. 

The  first  relates  to  the  estimate  of  reference  noise  level 
in  1)  above.  Presumably,  this  estimate  could  derive  from  the  mea- 
surements made  during  the  processing  for  detection  of  signal  pre- 
sence. It  could  also  derive  from  measurements  made  subsequent  to 
that  point  if  collection  of  more  noise  measurements  were  found  to 
provide  a decided  advantage.  In  any  event,  it  should  be  clear 
from  the  nature  of  the  decision  variable  that  the  basic  prob- 
abilistic conclusions  described  above  would  not  be  greatly  affected 
by  any  but  very  large  fractional  errors  in  estimating  T-  There 
seem  to  be  several  avenues  for  guaranteeing  that  such  large  errors 
will  be  extremely  improbable. 

The  second  consideration  relates  to  the  AGC  which  must  be 
kept  from  changing  appreciably  when  the  sync  portion  of  the  pre- 
amble enters  the  receiver.  If  necessary,  to  hold  the  receiver  AGC 
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constant,  one  may  transmit  out-of-band  tones  that  do  not  get 
through  the  matched  filter,  a procedure  that  may  also  be 
employed  following  the  Doppler  estimation  procedure. 


The  third  question  relates  to  the  effect  of  multipath  and 
equipment  distortion  on  the  sync  preamble  peak-to-sidelobe  ratio 
and  on  threshold  setting.  The  existence  of  two  or  more  paths 
may  result  in  superposition  of  sidelobes  and  spurious  threshold 
exceedances.  In  such  a case  it  is  felt  that  the  high  perform- 
ance levels  described  here,  with  some  modification  of  operating 
parameter  values,  could  still  be  attained.  This  would  be 
achieved  via  several  avenues.  As  an  example  we  point  out  that 
the  length  of  the  synchronization  sequence  has  in  no  way  been 
optimized.  We  have  chosen  from  a very  small  set.  There  is 
every  reason  to  believe  that  going  to  longer  sequences  would 
allow  greater  flexibility  in  threshold  setting  and  substantially 
increase  pfrformance.  Further  examination  of  multipath  effects 
(and  methods  for  overcoming  them)  is  recommended. 

Though  of  importance,  the  AGC  and  noise  level  questions 
raised  above,  once  corrected  in  the  straightforward  manner  sug- 
gested, are  not  expected  to  impact  preamble  performance  to  the 
same  extent  as  the  fundamental  considerations  emphasized  in  this 
section . 
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SECTION  3 


CODING  FOR  A KG  SYNCHRONIZATION  SEQUENCE 

j 

In  this  section,  good  block  codes  are  evaluated  for  the 
protection  of  the  KG  synchronizing  sequence.  This  KG  sync 
sequence  may  actually  be  composed  of  the  combination  of  a KG 
synchronizing  sequence  and  a small  number  of  bits  used  for  net 
control.  Sequences  of  lengths  100,  128,  and  150  are  assumed. 

The  number  of  tones  used  to  code  the  sequence  varies  from  250 
to  4000  which  (assuming  approximately  40  tones  per  frame)  covers 
the  range  of  6 to  100  frames.  With  22.5  ms  per  frame,  coding 
delays  in  the  range  of  .135  second  to  2.25  seconds  are  considered. 
While  the  actual  specification  limits  the  KG  decoding  delay  to 
68  frames,  or  1.53  seconds,  performance  curves*  for  delays  up 
to  2.25  seconds  are  still  of  interest. 

Since  the  modulation  is  4-phase  DPSK  with  2 bits  per  tone, 

L • N 

the  number  of  tones  used  is  given  by  — - — , where  N is  the  block 
code  length  and  L is  the  order  of  diversity  combining  for  each 
coded  bit.  The  best  BCH  codes  are  chosen  for  the  simple  reason 
that  they  can  be  decoded  easily.  The  codes  and  diversities 
used  are  shown  in  Table  3-1.  Note  that  the  total  number  of 
bits  used  to  protect  the  128-bit  KG  sequence  is  a multiple  of 
252  rather  than  a multiple  of  250  which  is  used  for  the  100-bit 
KG  and  150-bit  KG  sequences. 


In  this  section,  performance  curves  for  bit  error  probabilities 
are  presented.  The  analyses  for  computing  the  bit  and  block 
error  rate  are  given  in  Appendix  D. 
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TABLE  3-1 


CODES  FOR  SYNCHRONIZATION  SEQUENCE  (4  = 1,  2,  4,  8) 


Length  of 
Sequence 
in  Bits 

Number  of  Tones 

Code 

Order  of 
Diversity 

100 

250 

(125 , 100 ;8)* 

4 

100 

250 

(250,100;46)+ 

2 

100 

250 

(500,100;120) 

1 

100 

5004 

(125, 100; 8)* 

84 

100 

5004 

(250,100;46)+ 

44 

100 

5004 

(500 , 100 ; 120) 

24 

100 

5004 

(1000,100; 341) 

4 

128 

252 

(252 , 128 ; 37) 

2 

128 

252 

(504, 128;  110  )* 

1 

128 

5044 

(252,128 ; 37) 

44 

128 

5044 

(504 , 128 ; 110) 

24 

128 

5044 

(1008 ,128 ; 253) 

4 

150 

250 

(250,150;27) 

2 

150 

250 

(500 , 150 ; 96)** 

1 

150 

5004 

(250,150;27) 

44 

150 

5004 

(500, 150;  96)** 

24 

150 

5004 

(1000 , 150 ; 245) 

4 

Obtained  from  (125,103;8)  code 
Obtained  from  (250,101;46)  code 
^Obtained  from  (504 , 131 ; 110)  code 
Obtained  from  (500,155 ; 9 6 ) code. 


** 
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To  illustrate  the  effects  of  fading  bandwidth  and  frame 

rate  on  the  KG  sequence  performance,  the  hard  decoding  results 

are  used  for  comparison.  The  results  plotted  in  Figs.  3.1 

through  3.7  are  limited  to  500  tones  and  hard  decoding.  Table 

3-2  shows  the  different  values  for  these  seven  figures.  As 

can  be  seen  from  these  figures,  the  best  codes  are  the  codes 

with  length  250  or  252  using  fourth-order  diversity  for  all 

cases.  The  signal-to-noise  ratio  per  received  tone  (E  /N„)  in 

-7  c u 

dB  for  each  code  to  reach  a bit  error  rate  of  10  for  two  frame 
rates  and  three  different  rms  fading  bandwidths  are  given  in 
Table  3-3.  In  the  worst  case,  i.e.,  B = 5 Hz  at  a frame  rate 
of  44.44  frames/second,  the  (250,150;27)  code  needs  an  addi- 
tional 1.7  dB  to  reach  10  when  compared  to  the  case  B = 0. 
Therefore,  the  effect  of  the  rms  fading  bandwidth  B on  the  per- 
formance is  quite  small.  This  is  especially  true  for  the 
stronger  codes. 

In  Fig.  3.8,  the  1000  or  1008  tone  KG  codes  are  compared 
for  the  slowly  fading  channel  only.  It  is  expected  that  the 
performance  is  fairly  independent  of  the  rms  fading  bandwidth  B. 
In  this  case,  the  codes  with  lengths  of  250  and  252  again  yield 
the  best  performance. 

In  the  remainder  of  this  section,  the  results  are  extended 
to  cover  both  soft  and  hard  decodings,  and  the  number  of  tones 
per  block  range  from  250  (or  252)  to  4000  (or  4032) . From  the 
results  given  in  Figs.  3.1  through  3.7  we  note  that,  for  a 
fixed  fading  bandwidth,  the  results  for  a frame  rate  = 44.44 
frames/second  are  inferior  to  those  that  are  obtained  for  a frame 
rate  = 75  frames/second . Hence,  in  the  following  discussion, 


3-3 


TABLE  3-2 


v 
r 


t 


FRAME  RATES  AND  RMS  FADING  BANDWIDTH 
FOR  FIGURES  3.1  THROUGH  3.7 


Figure 

Frame  Rate 
(frames/sec) 

RMS 

Fading  Bandwidth 
(Hz) 

3.1 

44.44  or  75 

0 

3.2 

44.44 

1 

3.3 

44.44 

3 

3.4 

44.44 

5 

3.5 

75 

1 

3.6 

75 

3 

3.7 

75 

5 
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TABLE  3-3 


SNR  PER  RECEIVED  TONE  REQUIRED  TO  REACH  AN 
ERROR  RATE  OF  10~7  USING  FOURTH-ORDER  DIVERSITY 
AND  HARD  DECODING 


RMS 

Fading 

Bandwidth 

(Hz) 

Frame  Rate 
(frames/ sec) 

(250,100;46) 

Code 

(252 , 128 ; 37) 
Code 

(250 , 150 ; 27) 
Code 

0 

44.44  or  75 

4.6  dB 

5.4  dB 

6.5  dB 

1 

44.44 

4.7  dB 

5.5  dB 

6.6  dB 

3 

44.44 

5.0  dB 

5.9  dB 

7.1  dB 

5 

44.44 

5.8  dB 

6.7  dB 

8.2  dB 

1 

75 

4.7  dB 

5.4  dB 

6.5  dB 

3 

75 

4.8  dB 

5.6  dB 

6.7  dB 

5 

75 

5.0  dB 

5.9  dB 

7.1  dB 

1 


Figure  3.1  Performance  of  500- or  504 -Tone  KG  Codes  over  a Slowly 

Fading  HF  Channel  with  4-Phase  DPSK  Modulation  (Hard  Decoding) 
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Figure  3.2  Performance  of  500-  or  504-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames/sec,  B = 1 Hz,  Hard  Decoding) 
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3.3  Performance  of  500-  or  504-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames/sec,  B = 3 Hz,  Hard  Decoding) 
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Figure  3.4  Performance  of  500-  or  504-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames /sec,  B = 5,  Hard  Decoding) 
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Figure  3.6  Performance  of  500-  or  504-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 75  frames/sec,  B = 3 Hz,  Hard  Decoding) 
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Figure  3.7  Performance  of  500-  or  504-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 75  frames/sec,  B = 5 Hz,  Hard  Decoding) 
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Figure  3.8  Performance  of  1000-  or  1008-Tone  KG  Codes  over  a Slowly 

Fading  HF  Channel  with  4-Phase  DPSK  Modulation  (Hard  Decoding) 
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only  the  worst  case  will  be  considered.  Also,  when  the  number 
of  tones/code  block  becomes  larger,  a large  amount  of  diversity 
combining  is  used  for  each  bit  in  a code  block.  As  a result, 
the  performance  is  fairly  independent  of  the  rms  Doppler  band- 
width B.  Thus,  only  the  case  B = 5 Hz , the  worst  case,  is 
shown  when  the  number  of  tones /block  > 1000. 

Table  3-4  lists  all  the  pertinent  information  for  Figs.  3.9 
through  3.26. 

For  hard  decoding,  the  three  codes  of  code  length  250  or 

252  yield  the  best  performance.  This  is  also  true  for  most  soft 

decoding  cases.  [The  noted  exceptions  are  the  soft  decoding 

estimates  for  the  (1000, 100; 341)  code.  For  the  case  of  4000  or 

4032  tones,  the  (1000 , 100 ; 341)  code  outperforms  the  (250,100;46) 

code  by  0.5  dB,  and  the  (500,150;96)  code  edges  the  (250, 150 ;96) 

code  by  0.3  dB  at  a bit  error  rate  of  10  ^ . However,  it  is  much 

easier  to  soft  decode  the  shorter  code  than  the  longer  c^de. 

Also,  a better  algorithm  can  be  applied  to  the  shorter  codes  to 

compensate  for  the  0.5  to  0.3  dB  estimated  loss  in  performance.) 

Therefore,  the  three  codes  of  length  250  or  252  are  recommended 

to  be  used  for  KG  codes.  In  Table  3-5,  the  signal-to-noise 

ratios  per  tone  (E  /Nn)  required  to  reach  a bit  error  rate  of 
-7  t u 

10  are  listed.  The  required  signal-to-noise  ratios  per  received 

tone  (E  /N„)  for  the  KG  sequences  to  reach  a bit  error  rate  of 
-7  C U 

10  vs.  the  number  of  tones  employed  for  the  KG  sequences  are 
plotted  in  Fig.  3.27  on  the  top  horizontal  scale.  For  the  same 
delay,  a savings  of  2 dB  in  signal-to-noise  ratio  results  if 
soft  decoding  is  employed  instead  of  hard  decoding.  Note  that 


TABLE  3-4 
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Figure  3.9  Performance  of  250- or  252-Tone  KG  Codes  over  a Slowly  Fading 
HF  Channel  with  4-Phase  DPSK  Modulation  (Hard  Decoding) 
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Figure  3.10  Performance  of  250- or  252-Tone  KG  Codes  over  a Time-Varying 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 44.44 
frames/sec,  B = 1 Hz,  Hard  Decoding) 
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Performance  of  250-  or  252-Tone  KG  Codes  over  a Time-Varying 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 44.44 
frames /sec,  B = 3 Hz,  Hard  Decoding) 
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Figure  3.12  Performance  of  250- or  252-Tone  KG  Codes  over  a Time-Varying 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 44.44 
frames /sec,  B = 5 Hz,  Hard  Decoding) 

3-19 


128 ;V10) 


Figure  3.13  Estimated  Performance  of  250-  or  252-Tone  KG  Codes  over  a 


Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(Soft  Decoding) 
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Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
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Figure  3.15  Estimated  Performance  of  250- or  252-Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
(Frame  Rate  = 44.44  frames/sec,  B = 3 Hz,  Soft  Decoding) 
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Figure  3.16  Estimated  Performance  of  250- or  252-Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
(Frame  Rate  = 44.44  frames/sec,  B = 5 Hz,  Soft  Decoding) 

3-23 


^ \ 

« \ \ 
v-M  \ \ 


\ \\ 


1G  1 1 r i r~l~r  i i~r  j n rrjrn  rpri  i j rr  i t j~r  t r i | i i r i | i t r i j r i 1 i pi  r i i | t rvi 


-1  \"  ' O \\ 

10  A ' \ \?  W 

r.  \ ' \ ' 


, ft  \ \ « \v 

10  ft  \V-  v»  \* 


10 


Figure 


0 1 2 3 4 5 S / 8 9 10  11 

TOTPL  SNR  PER  RECEIVED  TONE  CET/NOJ  IN  DB 

3.17  Estimated  Performance  of  500-  or  504-Tone  KG  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(Soft  Decoding) 
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Figure  3.18  Estimated  Performance  of  500- or  504-Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
(Frame  Rate  = 44.44  frames/sec,  B = 1 Hz,  Soft  Decoding) 

3-25 


r 


i n riTT  i r | I n r n r rri  t it i 1 i t ti  i n i i i i r t r r -r  i i i i t r i t i i r i i i i i i 


_V  "n  * - - -:y. 

\ " ' \ \ 

\ \\  \\  i 

\ \ S ' oT>  ( 


Bit  Err 

-"rorRat, 


\'  \ \«* 

'S,  \ «. 

>\  V \ \ \o 

-{s  \\  \ \YvO 

\i  i \ \<r 

'W  '■  v ^ 

- ■.  'V  ' "o 

- \ I \ \ ' WO 

: •>  v'  \ * mo* 

- \ \ \\ 

- \\\  'V 

• \ \ \w  f> 

\ v \y° 

i ■ \ \ 'i 
: \ \\  \ \ 


> '■  o 
..  \ o 


to  ' ' 
itO-  '■  l 

'.ro  '■  \ 


\ \ ' 
\ \ \ 


\\  \ \o 


i w\\r*  . 

I ',ro  ,uo  ' 


• V*  | V '• 

1 I.  n f. 


\ 1 \ro  ♦ 

' SCP  -,o  I' 


i E?  V 6 

’ i u* 


I'C  ' Si  \S  I o\  ~ 

rri  ! i , i _ \\  oo  i o 


to 

/~'l  1 £~ 

i-*  i on 

o\ 

Pi  '' 


rn  > • ob\ -o 

\ tO,  l!,  - 1 to 

i’,°  ■,  o ' ° r-v  t-1 ! 1 
If  . '.cl.'  '£■  ro  10  t® 

. \\g  » fils 

§\£  \ Wdb.v  g! ^ 


oo r>  i it"  on  - . >t  f r;t  , II 
ot'  iV* ; t'.fo  pt  i 


gi'o  a o',  ■ 

f?  I 8 U,  % \ ^ i I 


fi  rp  \ llo  *;i ; 

1 xr*  ft-*'  ° i l°*  ^ ' i 

t * O o ' "*  U1  1 


S ' » o|  o \ V 
f M o i' 


g.  t's 


TO  l ' CL  It*  ' , 

* 'S  I TO  V'«  ' , 

_ 'o.'1  i*  v ^ 1 I 

y V"i'r  ; ; 

fe  . • t 


« 1 1 r*_.'  i ' , t 

t fr1'  ■,  n V ' . 

1 1 I*1  \\  l 

-L-O.I-1--L-1_1-J_l  !1.1|_I  -1  -iHUlL  Jt-L.,1. . I ..LAjL-1  _ I J J1  I I I I L.J 


L = Order  of  Diversity 


-J—l — I— ! J .l._l_l  J-.L..l-.l-J-1-J.  1 I—  I I LJ-l- 


0 11 


TOTPL  SNR  PER  RECEIVED  TONE  I.ET/N0)  IN  DB 


Figure  3.19  Estimated  Performance  of  500- or  504-Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
(Frame  Rate  = 44.44  frames/sec,  B = 3 Hz,  Soft  Decoding) 
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Figure  3.20  Estimated  Performance  of  500-  or 504 -Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation 
(Frame  Rate  = 44.44  frames/sec,  B = 5 Hz,  Soft  Decoding) 
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Figure  3.21a  Performance  of  1000-  or  1008-Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames/sec,  B = 5 Hz,  Hard  Decoding) 
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Performance  of  1000- or  1008 -Tone  KG  Codes  over  a Time- 
Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
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Figure  3.23  Performance  of  2000- or  2016-Tone  KG  Codes  over  a Time-Varying 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 44.44 
frames /sec,  B = 5 Hz,  Hard  Decoding) 
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Figure  3.24  Estimated  Performance  of  2000- or  2016  Tone  KG  Codes  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames/sec,  B = 5 Hz,  Soft  Decoding) 
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Figure  3.25  Performance  of  4000- or  4032-Tone  KG  Codes  over  a Time- 
HF  Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 
frames/sec,  B = 5 Hz,  Hard  Decoding) 
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Figure  3.26  Estimated  Performance  of  4000-  or  4032-Tone  KG  Coding  over  a 
Time-Varying  HF  Channel  with  4-Phase  DPSK  Modulation  (Frame 
Rate  = 44.44  frames/sec,  B = 5 Hz,  Soft  Decoding) 
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TABLE  3-5 


REQUIRED  SIGNAL-TO-NOISE  RATIO  PER  TONE  (Et/N0)  IN  dB 

TO  REACH  A BIT  ERROR  RATE  OF  10-7  FOR  A TIME- VARYING  HF  CHANNEL 
WITH  4-PHASE  DPSK  MODUTATION  (FRAME  RATE  = 44.44  FRAMES/SEC, 
RMS  DOPPLER  BANDWIDTH  = 5 Hz) 


• 

Number  of 
Tones  Per 
Code  Block 

(250 , 100 ;46)  Code 

(252 ,128 ; 37  ) Code 

(250 ,150; 27 ) Code 

Hard 

Decoding 

Soft 

Decoding 

Hard 

Decoding 

Soft 

Decoding 

Hard 

Decoding 

Soft 

Decoding 

250  or 

252 

11.9 

6.5 

* 

8.1 

* 

10.5 

500  or 

504 

5.8 

2.6 

6.7 

3.7 

8.2 

5.1 

1000  or 
1008 

2.2 

-0.2 

2.9 

0.6 

3.9 

1.7 

2000  or 
2016 

-0.5 

-2.6 

0.1 

-1.8 

0.8 

-1.1 

4000  or 
4032 

-2.8 

-4.6 

-2.3 



-4.0 

-1.6 

-3.1 

Cannot  reach  a bit  error  rate  of  10  7 due  to  irreducible  errors. 
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DELAY  FOR  40  TONES /FRAME  FORMAT  (SEC) 
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Figure  3.27  Required  Signal-to-Noise  Ratio  for  the  KG  Sequence  to 
Reach  a Bit  Error  Rate  of  10"^  for  a Time-Varying  HF 
Channel  with  4-Phase  DPSK  Modulation  (Frame  Rate  = 
44.44  frames /sec,  Bandwidth = 5 Hz) 
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NUMBER  OF  TONES  EMPLOYED  FOR  KG  SEQUENCES 


SECTION  4 


MODEM  SIMULATION 


Conventional  parallel  tone  modems  operating  on  the  HF  channel 
generally  suffer  significant  performance  degradation  under  condi- 
tions in  which  large  multipath  spread  can  be  expected.  These 
modems  do  not  take  advantage  of  in-band  signal  diversity  during 
instances  of  large  multipath.  As  a result,  the  information  on 
a tone  which  is  selectively  eliminated  during  a fade  cannot  be 
recovered  from  the  information  in  other  tones  which  have  not  been 
so  catastrophically  affected. 

"Jc 

Modems  which  employ  coding  attempt  to  explicitly  use  in-band 
diversity  by  adding  redundant  information  to  the  transmitted  tone 
array.  In  order  to  demonstrate  the  usefulness  of  coding  021  the 
HF  channel,  three  2400-bit  per  second  modems  will  be  defined  to 
form  the  basis  of  a comparative  study: 

(1)  A 16  parallel  tone,  4-phase  DPSK  modem  will  be 
used  when  transmissions  not  involving  coding  are 
to  be  tested.  Each  tone  is  separated  by  110  Hz 
from  adjacent  tones.  The  modem  has  a frame 
duration  of  13.33  ms. 

(2)  A 39  parallel  tone,  4-phase  DPSK  modem  will  be 
used  for  simulations  involving  the  Golay  code 
(referred  to  as  ANDVT-G).  Each  tone  is  separated 
from  adjacent  tones  by  55.56  Hz.  The  modem  has 

a frame  duration  of  22.5  ms. 


Information  concerning  code  selection  is  given  in  Appendix  B. 
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(3)  A 41  parallel  tone,  4-phase  DPSK  modem  will  be 
used  for  simulations  involving  a multiple  rate 
coding  structure  (referred  to  as  ANDVT-MRC) . 

Each  tone  is  separated  from  adjacent  tones  by 
35.56  Hz.  The  modem  has  a frame  duration  of 
22.5  ms . 

Each  of  these  modems  will  be  evaluated  using  two  HE  channel 
models.  (A  detailed  description  of  the  HF  channel  model  can  be 
found  in  Appendix  A.) 

Model  1 - The  HF  channel  is  represented  by  two  equi- 

strength  paths  separated  by  1 ms . Each  path 
is  assumed  to  have  a double-sided  Doppler 
bandwidth  of  1.0  Hz.  This  channel  is  repre- 
sentative of  channel  behavior  when  operating 
near  the  maximum  usable  frequency  (MUF). 

Model  2 - The  HF  channel  is  represented  by  a 6-path 
equi-strength  split  ray  model  with  a total 
multipath  spread  of  2.6  ms.  Each  path  is 
assumed  to  be  resolvable  into  two  components 
separated  by  0.1  ms.  The  double-sided 
Doppler  bandwidth  per  path  is  1.5  Hz.  This 
channel  is  representative  of  channel  behavior 
when  operating  significantly  below  the 
maximum  usable  frequency  (MUF). 

A functional  block  diagram  of  the  HF  channel  simulation 
appears  in  Figure  4.1.  After  the  input  data  sequence  has  been 
encoded,  the  resulting  vector  is  modulated,  transmitted  over  a 
synthetic  channel,  and  demodulated.  The  exclusive-OR  (®)  of  the 
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Figure  4.1  Simulation  Block  Diagram 


transmitted  vector  and  received  vector  yields  a bit-by-bit 
indication  of  the  errors  incurred  during  data  transmission. 

An  exclusive-OR  operation,  using  the  decoded  error  vector  and 
the  input  data  sequence,  yields  the  output  data  sequence. 

The  simulation  program  generically  diagrammed  in  Figure  4.1 
is  capable  of  emulating  several  coding  schemes  (e.g.,  no  coding, 
Golay  code,  BCH  code,  multiple  rate  code)  as  well  as  most  parallel 
tone  modems  of  practical  interest.  Throughout  the  development  of 
the  simulation,  considerable  emphasis  has  been  placed  on  main- 
taining flexibility  in  order  to  insure  the  usefulness  of  these 
programs  as  a research  tool. 

The  data  to  be  transmitted  over  the  HF  channel  consists  of 
LPC-10  quantized  parameters.  The  parameter  content  of  each  frame 
analyzed  by  the  LPC-10  processor  appears  in  Figure  4.2.  There 
are  54  information  bits  for  each  analyzed  frame.  For  the  LPC-10 
analyzer,  a frame  consists  of  22.5  ms  of  data.  Since  the  16-tone 
modem  uses  4-DPSK  modulation,  only  32  bits  of  information  are 
conveyed  by  a given  set  of  tones.  As  a result,  a frame  of  LPC-10 
data  is  not  synchronous  with  the  16  parallel  tone  modem.  Each  of 
the  other  modems  (ANDVT-G  and  ANDVT-MRC)  uses  a tone  complement 
which  insures  frame  synchronism  for  the  particular  coding  structure 
employed . 

4 . 1 ANDVT  Modem  with  Golay  Code  Bit  Protection 

A 39  parallel  tone  modem  using  4 DPSK  can  convey  up  to  78  bits 
per  frame.  Since  54  information  bits  comprise  the  LPC-10  quantized 
parameters,  there  are  24  residual  bits  which  can  be  used  for  the 
parity  bits  of  the  Golay  codes  employed.  A Golay  code  with  overall 
parity  has  dimensions  (24,12;8).  By  using  two  Golay  codes  to 
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protect  24  bits  of  the  LPC-10  quantized  parameters,  the  full 
complement  of  available  bits  (54  + 2 x 12)  in  the  ANDVT-G  modem 
is  utilized.  : 

For  the  purposes  of  simulation,  the  LPC-10  quantized  para- 
meters have  been  restructured.  An  attempt  has  been  made  to  order 
the  bits  comprising  these  parameters  according  to  their  relative 
importance  in  synthesized  speech  quality.  The  resulting  ordering 
is  portrayed  in  Figure  4.3.  The  two  Golay  codes  employed  by  the 
modem  are  given  the  task  of  protecting  the  first  24  bits  of  the 
ordered  bit  sequence.  A summary  of  the  parameters  to  which  the 
protected  bits  belong  appears  as  Case  1 in  Figure  4.2.  Case  2 
was  attempted,  however,  speech  tapes  indicated  that  Case  1 offers 
superior  performance. 

Once  the  bits  to  be  protected  have  been  encoded,  all  the  bits 
must  be  placed  on  the  39-tone  array.  In  order  to  benefit  from  any 
available  in-band  diversity,  care  is  exercised  in  maintaining  a 
large  tone  separation  between  bits  representing  the  same  codeword. 
Of  particular  importance  is  the  assurance  that  two  bits  from  one 
codeword  do  not  occupy  the  same  tone  (since  both  bits  may  be  lost 
during  a fade) . 

The  mapping  between  the  78  transmitted  bits  and  the  39  tone 
array  is  detailed  in  Figure  4.4.  For  convenience,  this  mapping 
will  be  referred  to  as  the  modulation  format  bit  assignments.  The 
first  12  bits  of  the  ordered  bit  sequence  (Figure  4.3)  are  pro- 
tected by  the  Golay  code  in  FIELD  1 of  Figure  4.4.  The  second  set 
of  12  bits  in  the  ordered  bit  sequence  is  protected  by  the  Golay 
code  in  FIELD  2 of  Figure  4.4.  The  remaining  30  bits  in  the 
ordered  sequence  are  unencoded.  Collectively,  they  comprise 
FIELD  3 in  Figure  4.4. 
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ed  Sequence  of  Bits  Comprising  the  LPC-10  Quantized  Parameters 


12  BITS  , 12  BITS  12  BITS  , 12  BITS  30  BITS 


Figure  4.4  Modulation  Format  Bit  Assignment  (Golay  Code) 


4 . 2 ANDVT  Modem  with  Multiple  Rate  Code  Bit  Protection 

Simulations  have  also  been  conducted  for  a modem  (ANDVT-MRC) 
which  uses  a multiple  rate  coding  structure.  This  coding  approach 

I 

places  redundant  information  across  time  frame  boundaries  as  well 
as  the  tone  array.  A detailed  discussion  of  this  time  domain 
coding  technique  can  be  found  in  [4.1]. 

For  the  purposes  of  this  simulation,  it  is  only  necessary  to 
understand  the  impact  of  the  codeword  structure  on  the  LPC-10 
bit  protection  choices.  Toward  that  end,  a cursory  description 
of  a multiple  rate  coding  technique  seems  appropriate. 

In  Figure  4.5,  a typical  multiple  rate  codeword  is  depicted. 
In  general,  each  multiple  rate  codeword  contains  the  modulo-2 
(exclusive-OR)  addition  of  several  inner  codewords.  With  little 
loss  in  generality,  assume  that  each  inner  codeword  is  the  result 
of  a BCH  encoding  of  a set  of  information  bits.  The  inner  code- 
words forming  a multiple  rate  codeword  will  have  been  encoded  at 
more  than  one  rate  (the  ratio  of  the  number  of  information  bits 
in  the  codeword  to  the  total  length  of  the  codeword) . This  simu- 
lation only  considers  code  structures  in  which  two  code  rates  are 
needed . 

In  Figure  4.5,  there  is  a total  of  K1  + K2  information  bits 
in  the  codeword.  The  information  section  which  is  to  receive  the 
greatest  protection  contains  Kl  bits.  The  other  information 
section  contains  K2  bits;  it  is  afforded  a somewhat  smaller  pro- 
tection level.  In  an  ordinary  single  rate  implementation,  the 
information  bits  are  encoded  to  form  an  N-bit  codeword  with  rate 
(Kl  + K2)/N.  In  this  case,  both  information  sections  receive  the 
same  level  of  protection. 
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Figure  4.5  Contents  of  a Multiple  Rate  Codeword 


The  multiple  rate  code  structure  incorporates  the  Kl  most 
important  bits  into  a codeword  with  rate  (2Kl  + K2)/N.  This 
higher  rate  encoding  is  accomplished  by  appending  the  Kl  most 
important  bits  to  Kl + K2  zeros.  The  resulting  conglomerate 
becomes  the  2K1  + K2  information  bits  of  the  high  rate  codewords. 
After  encoding,  the  high  rate  vector  is  added  (modulo-2).  to  a 
predefined  set  of  future  codewords.  Observe  that  the  Kl  most 
important  information  bits  which  have  been  included  in  high  rate 
codewords  appear  in  the  parity  section  of  the  low  rate  encoding 
of  future  information  blocks.  The  number  of  times  that  the  high 
rate  code  is  fed  forward  is  equivalent  to  the  burst  error-correcting 
capability  of  the  multiple  rate  code.  The  description  in  Figure  4.5 
and  the  structure  depicted  in  Figure  4.6  typify  the  feed-forward 
character  of  this  time  domain  coding  approach. 

The  resulting  admixture  of  different  rate  codewords  can  be 
decoded  in  the  low  rate  mode  in  the  absence  of  errors  or  in  the 
high  rate  mode  when  a past  block  of  Kl  bits  is  unknown.  Of  course, 
the  low  rate  code  (N,  K1  + K2)  must  be  a subcode  of  the  high  rate 
code  (N , 2K1  + K2).  This  is  generally  not  a problem  for  BCH  codes 
since  they  can  be  selected  to  insure  that  the  generator  polynomial 
for  the  high  rate  code  is  a multiple  of  the  generator  polynomial 
for  the  low  rate  code. 

During  the  decoding  process,  the  Kl  most  important  informa- 
tion bits  can  potentially  be  recovered  from  any  of  the  codewords 
which  they  influence.  If  all  the  information  from  previous  code- 
words has  not  been  removed,  then  the  decoding  must  be  attempted 
in  the  high  rate  mode.  When  the  information  from  all  but  one 
previous  frame  has  been  removed  from  a codeword,  and  when  a 
successful  high  rate  decode  can  be  achieved,  the  Kl  most  important 
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Figure  4.6  Structure  of  a Triple  Burst  Correcting  Multiple  Rate  Code 
(K^  information  bits  are  principally  protected  by  the 
high  rate  code  and  the  multiple  repeats;  K2  information 
bits  are  not  protected  by  multiple  repeats;  H bits  are 
the  parity  check  bits  of  the  high  rate  code;  and  there 
are  U unprotected  bits,  t is  a delay  constrained  to  be 
an  integral  multiple  of  a frame  duration.) 

4-12 


issxsm 


information  bits  from  the  previous  frame  can  be  recovered.  A 
low  rate  decode  can  be  tried  only  when  the  information  from  all 
previous  codewords  has  been  removed.  The  K2  information  bits 
in  a previous  codeword  can  never  be  recovered  if  they  have  not 
been  correctly  decoded  initially.  Thus,  they  are  not  as  well 
protected  as  the  Kl  most  important  information  bits. 

As  previously  mentioned,  each  frame  of  the  LPC-10  quantized 
parameters  consists  of  54  information  bits.  If  an  LPC-10  frame 
is  to  remain  synchronous  with  a modem  frame  when  a multiple  rate 
code  structure  is  employed,  then  54-(Kl  + K2)  bits  must  be  trans- 
mitted unencoded.  The  presence  of  unencoded  bits  is  reflected  in 
the  structure  of  Figure  4.6.  With  two  bits  per  tone  in  a 4-DPSK 
modem,  [N  + 54-  (Kl  + K2)]/2  tones  are  required.  Since  the  ANDVT-MRC 
is  presently  conceived  as  a 41  parallel  tone  modem,  the  number  of 
parity  bits  in  the  low  rate  code  [N  - (K1  + K2)]  is  constrained  to 
be  28. 

Initially,  two  multiple  rate  codes  are  being  considered. 

Each  code  will  be  defined  by  a descriptor  of  the  form  (N,k  ;d  ) — 
(N,k^;d^);  where  (k^)  are  the  number  of  information  bits  in 
the  low  (high)  rate  code,  and  where  d^  (d^)  i-s  minimum  dis- 

tance of  the  low  (high)  rate  code.  The  values  of  Kl  and  K2  for 
a (43 , 15 ; 12)  -* (43 , 30 ;6)  code  are,  respectively,  15  and  0.  Simi- 
larly, the  values  of  Kl  and  K2  for  a (64, 36 ; 12)  ^ (64 , 51 ;6)  code 
are,  respectively,  15  and  21.  A summary  of  the  parameters  to 
which  the  protected  bits  belong  appears  in  Figure  4.7. 

After  each  multiple  rate  codeword  has  been  formed,  the 
resulting  vector  along  with  the  unencoded  bits  must  be  placed 
in  the  41  tone  array.  The  modulation  format  bit  assignments 
for  the  (43 , 15 ; 12)  — (43 , 30 ;6 ) multiple  rate  code  are  clearly 


I 


4-13 


Figure  4.7  List  of  Protected  Bits  for  LPC-10  Multiple  Rate  Code  Structure 


defined  in  Figure  4.8.  Similar  bit  assignments  for  the  (64, 36; 12) 
(64,51;6)  multiple  rate  code  are  contained  in  Figure  4.9. 


Ordinarily,  more  than  one  bit  from  the  same  codeword  is  not 
placed  on  a single  transmitted  tone.  This  restriction  is  ob- 
viously not  applicable,  however,  when  the  length  of  the  codeword 
exceeds  the  number  of  tones  available.  The  procedure  followed  in 
the  preparation  of  Figures  4.8  and  4.9  involves  maximizing  the 
tone  separation  of  the  bits  within  each  of  the  defined  fields. 

When  a successful  decode  does  not  occur,  the  procedure  minimizes 
the  variance  of  the  number  of  bit  errors  per  frame  for  a given 
coding  field. 

4 . 3 Simulation  Results 

For  each  of  the  two  HF  channel  models  described  in  the  intro- 
duction to  Section  4,  simulations  have  been  conducted  using  the 
following:  1)  a 16  parallel  tone  modem  without  coding;  2)  the 

ANDVT  modem  with  the  Golay  code;  and  3)  the  ANDVT  modem  with  the 
(64, 36 ; 12) — (64 , 51 ;6)  multiple  rate  code.  A principal  product  of 
these  simulations  has  been  the  production  of  a set  of  voice  tapes 
which  document  the  performance  of  these  modems  as  a function  of 
E^/N^  (s:*-8na-*--t:o“no:‘-se  rati°  Per  information  bit).  In  Appendix  C 
the  relationship  between  E^/Nq  and  E /Nq  (signal-to-noise  ratio 
per  tone)  has  been  explored. 

\ 

Of  more  immediate  interest  are  the  statistical  summaries 
detailing  modem  performance  as  a function  of  E^/Nq.  In  Figure  4.10 
comparative  bit  error  rates  are  presented  for  the  16  parallel  and 
the  ANDVT  modem  using  the  Golay  code;  a two-path  HF  channel  is 
assumed.  Similar  results  appear  in  Figure  4.11  for  a 6-path  split 
ray  HF  channel  model.  The  raw  error  rate  before  decoding  for  the 
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ANDVT  modem  applies  to  the  54  bits  which  comprise  the  LPC-10 
quantized  parameters.  The  curves  which  present  results  for 
the  Golay  hard  decoder  and  the  Golay  soft  decoder  (Chase  algo- 
rithm 2)  obviously  apply  only  to  the  24  bits  which  are  protected 
by  the  Golay  code.  Coding  gains,  relative  to  raw  bit  error  rate, 
are  significantly  larger  for  the  6-path  split  ray  model  than  for 
the  2-path  model.  This  is  directly  attributable  to  the  greater 
in-band  diversity  for  the  6-path  split  ray  model. 


Bit  error  rates  for  a 16  parallel  tone  modem  and  the  ANDVT 
modem  employing  a multiple  rate  code  are  plotted  in  Figure  4.12; 
a 2-path  HF  channel  is  assumed.  Figure  4.13  contains  comparative 
bit  error  rates  for  conditions  under  which  a 6-path  split  ray  HF 
channel  model  is  appropriate.  In  all  the  multiple  rate  code 
simulations  conducted,  a burst  error-correcting  capability  of  1 
has  been  assumed  for  the  code.  In  addition,  the  coding  has  been 
implemented  with  a decoding  delay  equivalent  to  the  duration  of 
ten  LPC-10  frames.  With  reference  to  the  structure  in  Figure  4.6, 
the  delay  parameter  T has  a value  equal  to  the  duration  of  ten 
LPC-10  frames. 

Statistics  are  kept  separately  for  the  bits  which  receive 
first  and  second  levels  of  protection  from  the  multiple  rate  code. 

The  annotation  on  the  curves  of  Figures  4.12  and  4.13  is  consis- 
tent with  the  notation  of  Section  4.2;  Kl  denotes  the  set  of  bits 
which  receive  the  first  level  of  protection,  while  K2  denotes  the 
set  of  bits  which  receive  the  second  level  of  protection.  Again, 
relative  coding  gains  are  larger  for  the  6-path  split  ray  HF  model 
than  for  the  2-path  HF  model.  This  is  dramatically  exemplified  by 
comparing  the  performance  of  the  soft  decoder  acting  on  the  bits 
receiving  the  second  level  of  protection  with  the  performance  of 
the  hard  decoder  acting  on  the  bits  receiving  the  first  level  of 
protection . 
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It  is  of  some  interest  to  determine  which  of  the  two  prin- 
cipal coding  structures  (Golay,  multiple  rate)  performs  best  in 
conjunction  with  the  ANDVT  modem.  From  Figures  4.9  through  4.13 
it  is  clear  that  the  Golay  code  outperforms  the  multiple  rate 
code  at  low  E^/Nq.  However,  at  higher  E^/N^,  soft  decoding  the 
bits  which  receive  the  multiple  rate  codes'  first  level  of  pro- 
tection produces  lower  bit  error  rates  than  soft  decoding  the 
bits  protected  by  the  Golay  code.  The  crossover  point  occurs  at 
a lower  E^/Nq  for  the  2-path  HF  model  than  for  the  6-path  split 
ray  model.  No  conclusions  can  be  drawn  from  these  observations, 
however,  since  the  two  coding  structures  protect  different  numbers 
of  bits.  The  final  measure  of  relative  value  must  be  speech- 
listening quality. 
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SECTION  5 


ANDVT  IMPLEMENTATION  STUDY 

5 . 1 Introduction 

This  section  describes  the  results  of  a study  of  the  imple- 
mentation of  the  ANDVT  modem  by  a programmable  processor.  This 
study  was  undertaken  to  satisfy  the  contract  requirement  that  the 
proposed  modem  algorithms  be  checked  for  implementation  prac- 
ticality and  has  influenced  the  algorithms  selected.  The  study 
focusses  on  obtaining  estimates  of  the  processor  time  required  to 
perform  each  of  the  modem  tasks.  The  lack  of  timely  information 
concerning  candidate  processors  necessitated  the  use  of  a hypo- 
thetical processor  model.  This  model  was  selected  to  be  simple 
enough  so  that  it  would  not  exceed  the  capabilities  of  the  candi- 
date processors  and  yet  to  be  general  enough  so  that  the  perform- 
ance of  the  candidate  processors  could  be  easily  related  to  that 
of  the  model.  A description  of  this  model  may  be  found  in 
Section  5.3. 

The  lack  of  knowledge  of  the  input/output  and  operating  sys- 
tems precludes  the  assumption  of  a sophisticated  multi-tasking 
operating  environment.  The  study  has  thus  been  limited  to  a single 
task  mode  appropriate  to  a half-duplex  device.  Since  the  modu- 
lator portion  of  the  modem  is  much  less  complex  than  the  demodu- 
lator, the  peak  processing  time, and  thus  the  limits  to  performance, 
will  lie  with  the  demodulator  implementation.  Because  of  this  and 
the  limited  amount  of  time  for  the  accomplishment  of  this  study, 
only  the  demodulator  implementation  has  been  studied  in  detail. 
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In  the  sequel  the  modulator  will  be  mentioned  only  in  regards 
to  the  demodulator  functions  of  which  it  performs  the  inverse 
function . 

The  remaining  sections  will  briefly  review  the  operation 
of  the  modem,  describe  the  assumptions  made  in  the  study,  and 
present  the  results  of  the  study.  Section  5.2  provides  the 
review  of  the  modem.  Section  5.3  describes  the  processor  model, 
and  Section  5.4  presents  the  study  of  the  demodulator  algor- 
ithms. A summary  of  the  results  and  an  estimate  of  the  memory 
requirements  is  given  in  Section  5.5. 

5 . 2 Modem  Review 

5.2.1  Signaling  Format 

The  ANDVT  modem  is  a multi-tone  modem  with  a frame  time 
compatible  with  the  LPC-10  speech  processor.  The  tone  modu- 
lation is  four-phase  DPSK.  The  transmitted  phase  at  the  begin- 
ning of  a frame  is  defined  as: 


where 


0,  =0,  - + 0 

k k-i  s 


0^.  = tone  phase  at  beginning  of  current  frame 
0^._^  = tone  phase  at  beginning  of  previous  frame 


0 = phase  shift  due  to  data, 

s 


The  values  for  0 given  the  two  data  bits  x, xn  modulating  the  tone 
are  specified  in  Table  5-1.  The  tone  library  may  consist  of 
about  40  tones  (depending  on  the  error  correcting  coding  format) 
and  an  empty  frequency  slot  among  the  active  tones.  This  empty 
slot  is  used  by  the  receiver  to  maintain  proper  frame  synchroni- 
zation. No  Doppler  reference  tone  is  present  because  the  receiver 
will  recover  Doppler  tracking  information  from  the  data  phase. 

In  initiating  communications,  the  transmitter  will  transmit  a 
preamble  to  allow  synchronization  of  the  receiver.  The  preamble 
consists  of  three  sequences:  a four-tone  Doppler  reference  sequence, 
a frame  sync  sequence  of  wideband  pulses,  and  a phase  reference 
frame  with  about  40  tones.  The  four-tone  Doppler  reference  sequence 
announces  the  presence  of  a signal  to  the  receiver  and  allows  the 
receiver  to  estimate  the  Doppler  offset  of  the  incoming  signal. 

The  frame  sync  sequence  consists  of  a wideband  signal  modulated 
by  a sequence  whose  autocorrelation  function  exhibits  a high  peak- 
to-sidelobe  ratio.  The  detection  of  these  peaks  allows  the 
receiver  to  find  the  proper  locations  of  the  frame  boundaries. 

The  phase  reference  frame  establishes  the  starting  phase  of  the 
DPSK  modulation.  The  differential  phase  changes  of  the  successive 
frame  begins  the  transmission  of  data. 

5.2.2  Sampling  Rate  Considerations 

For  a digital  implementation  of  a multi-tone  HF  modem  it  is 
desirable  that  the  chosen  sampling  rate  result  in  an  integral 
number  of  samples  over  both  the  frame  time  and  the  integration 
time.  For  the  ANDVT  modem  frame  time  of  22.5  ms,  the  acceptable 
range  of  integration  times  is  approximately  16  ms  (a  1.5  dB  loss 
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in  signal  power)  to  20  ms  (a  2.5  ms  guard  time).  An  additional 
constraint  due  to  the  planned  use  of  a 64  point  DFT  is  that  the 
integration  time  contain  64  samples. 

Let  the  number  of  samples/frame  be  64  + n 


T = Integration  time 


64  + n 


22.5  ms 


16  ms  < T-j.  < 20  ms 


8 < n < 26 


The  following  relations  can  also  be  developed: 


T = Guard  time  = 77 — ; 22.5  ms 

G 64  + n 


f = Sampling  rate  of  analytic  baseband 


1 64  + n 

signal  = ■ - - 

Z Z • j m s 


Af  = Tone  spacing  = fg/64 


N = Number  of  tones  = BW/Af 


BW  = Channel  bandwidth. 


The  quantities  are  ennumerated  in  Table  5-2  for  the  allowable 
values  of  n and  three  values  of  BW. 


5.2.3  Modulator  Functions 

A flow  chart  of  the  modulator  is  given  in  Fig.  5.1.  At  the 
start  of  a message  four  signal  sequences  must  be  sent  to  allow 
the  demodulator  to  acquire,  synchronize,  and  properly  decode  the 
message.  These  are: 

• Doppler  preamble 

• Frame  sync  preamble 

• Phase  reference  frame 

• Cryptographic  key  and  net  control  frames. 

These  are  followed  by  a sequence  of  data  frames. 

5. 2. 3.1  Doppler  Preamble 

The  Doppler  preamble  consists  of  four  unmodulated  tones.  It 
announces  the  start  of  a message  and  provides  a reference  whereby 
the  amount  of  Doppler  offset  can  be  measured  at  the  receiver. 

One  period  of  this  sequence  is  generated  by  an  inverse  DFT  opera- 
tion on  the  desired  spectrum.  This  is  then  repeated  for  the  dura- 
tion of  the  Doppler  preamble. 

5. 2. 3. 2 Frame  Sync  Preamble 

The  exact  time  at  which  the  receiver  declares  signal  pre- 
sence will  vary  with  the  channel  conditions  and  the  AGC  parameters 
of  the  radio.  Consequently  the  signal  must  contain  a reference 
from  which  the  receiver  can  recover  the  boundaries  of  the  succeed- 
ing frames.  This  function  is  accomplished  by  the  frame  sync  pre- 
amble. The  frame  sync  preamble  consists  of  a complex  sequence 
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TABLE  5-2 


ALLOWABLE  SAMPLING  RATES 

I 


n 

TI 

f 

s 

Af 

2400 

Hz 

2700 

Hz 

3000 

Hz 

8 

20 

2.5 

3200.00 

50.00 

48 

54 

60 

9 

19.73 

2.77 

3244.44 

50.69 

47 

53 

59 

10 

19.46 

3.04 

3288.89 

51.39 

46 

52 

58 

11 

19.20 

3.30 

3333.33 

52.08 

46 

51 

57 

12 

18.95 

3.55 

3377.78 

52.78 

45 

51 

56 

13 

18.70 

3.80 

3422.22 

53.47 

44 

50 

56 

14 

18.46 

4.04 

3466.67 

54.17 

44 

49 

55 

15 

18.23 

4.27 

3511.11 

54.86 

43 

49 

54 

16 

18.00 

4.50 

3555.56 

55.56 

43 

48 

54 

17 

17.78 

4.72 

3600.00 

56.25 

42 

48 

53 

18 

17.56 

4.94 

3644.44 

56.94 

42 

47 

52 

19 

17.35 

5.15 

3688.89 

57.64 

41 

46 

52 

20 

17.14 

5.36 

3733.33 

58.33 

41 

46 

51 

21 

16.94 

5.56 

3777.78 

59.03 

40 

45 

50 

22 

16.74 

5.76 

3822.22 

59.72 

40 

45 

50 

23 

16.55 

5.95 

3866.67 

60.42 

39 

44 

49 

24 

16.36 

6.14 

3911.11 

61.11 

39 

44 

49 

25 

16.18 

6.32 

3955.56 

61.81 

38 

43 

48 

26 

16.00 

6.50 

4000.00 

62.50 

38 

43 

48 

1 
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modulating  a wideband  pulse.  The  modulating  sequence  consists  of 
a few  repetitions  of  a sequence  selected  to  have  an  autocorrela- 
tion function  with  a high  peak-to-sidelobe  ratio.  The  receiver 
uses  a pulse  filter  matched  to  the  wideband  pulse  shape  to  recover 
the  sequence.  The  correlation  can  then  be  computed  to  locate  the 
frame  sync  boundary. 

5. 2. 3. 3 Phase  Reference  Frame 

The  modulation  format  described  in  Section  5.2.1  is  four- 
phase  DPSK.  Since  the  information  is  carried  by  the  frame-to- 
frame  phase  difference,  the  first  frame  must  provide  the  refer- 
ence for  the  next  frame.  The  64-point,  complex  FFT  routine  of 
the  demodulator  is  used  to  perform  an  FFT  ^ and  produce  128  real 
samples  in  the  integration  time  of  the  frame.  This  sequence  is 
extended  partially  to  fill  out  the  frame.  The  initial  values  for 
tone  phase  are  selected  to  minimize  the  ratio  of  peak-to-average 
power  in  the  transmitted  signal. 

5 . 2 . 3 . 4 Cryptographic  Key  and  Net  Control  Frames 

A binary  sequence  specifying  the  cryptographic  key  and  the 
net  control  designator  are  heavily  encoded  for  error  correction 
purposes  and  modulate  a number  of  frames  as  per  Section  5.2.1. 

The  codes  used  are  described  in  Section  3 of  this  report.  They  are 
easily  implemented  using  a linear  feedback  shift  register  to 
encode  the  parity  section. 

5. 2. 3. 5 Data  Frames 

Each  data  frame  lasts  for  one  frame  time  (22.5  msec)  and 
contains  the  encrypted  speech  parameters  for  the  corresponding 
speech  frame.  The  most  critical  bits  of  the  speech  parameters 


are  protected  by  a block  code  prior  to  the  modulation  of  the  tone 
phases.  The  code  improves  the  quality  of  the  reconstructed  speech 
in  the  presence  of  channel  errors.  After  each  frame  the  modulator 
checks  for  end  of  message  and  decides  whether  or  not  to  continue 
transmitting  the  data  frames. 

5. 2. 3. 6 FFT~1  Implementation 

The  modulator  shares  the  FFT  routine  of  the  demodulator  and 
uses  it  to  translate  from  the  desired  tone  phases  to  the  real 
samples  representing  the  signal  in  the  time  domain.  While  the 
demodulator  uses  an  FFT  routine,  the  duality  of  the  discrete 
Fourier  transform  indicates  that  an  inverse  FFT  can  be  implemented 
by  transforming  the  reverse  sequence  of  the  frequency  components 
and  dividing  by  N,  the  number  of  samples. 


X(k)  = DFT  [x(n) } 


where 


X(k)  = frequency  sample 
x(n)  = time  samples. 


Then 


x(n)  = DFT"1  (X(k)}  = | DFT  (X(-k)} 


Figure  5.2  shows  a flow  graph  of  the  means  by  which  the  128 
time  samples  are  produced  using  the  64-point  DFT  The  variables 
in  this  graph  are: 
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= Frequency  samples 


I 

X£  = Conjugate  of 


-i2tr/2N 

e 


x = Time  samples 
n 


i = J-l 


This  method  uses  the  fact  that  the  spectrum  of  a real  sequence 
(of  length  2N)  has  conjugate  symmetry  (X2N-k  = to  reduce  the 

DFT  size  requirement  from  2N  points  to  N points.  In  this  fashion 
both  the  real  and  imaginary  parts  of  the  output  contain  informa- 
tion. This  technique  is  used  to  increase  the  commonality  of  the 
modulator  and  demodulator  and  reduce  the  storage  requirements. 

5.2.4  Demodulator  Functions 

A flow  chart  of  the  demodulator  is  given  in  Fig.  5.3.  The 
boxes  indicate  the  various  tasks  the  demodulator  performs.  These 
are: 

• Signal  presence  detection 

• Doppler  estimation 

• Frame  sync  estimation 

• Phase  reference  frame  processing 

• Cryptographic  key  and  net  control  frame  processing 

• Data  frame  processing. 


j 


li 


i 
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Signal  Presence 
Detection 


T! 


While  quiescent  the  demodulator  monitors  the  spectrum  of  the 
incoming  signal  to  determine  the  presence  of  a signal.  When  a 
signal  is  present  the  demodulator  begins  the  processing  of  the 
preamble.  The  first  portion  of  the  preamble  consists  of  four 
unmodulated  tones.  These  tones  are  processed  to  provide  an  esti- 
mate of  the  Doppler  offset  of  the  signal.  This  estimate  is  used 
in  the  following  stages  to  remove  the  Doppler  offset  prior  to 
processing.  The  second  stage  of  the  preamble  is  used  to  estimate 
the  location  of  the  frame  boundaries  in  the  frame-oriented  pro- 
cessing which  follows.  This  allows  the  modulated  tones  in  these 
frames  to  be  orthogonal  over  the  integration  time.  The  first 
frame-oriented  signal  is  the  phase  reference  frame.  This  frame 
provides  the  tone  phase  reference  needed  to  demodulate  the  next 
frame  in  the  DPSK  format.  The  next  set  of  frames  contain  a 
cryptographic  key  encoded  for  error  correction.  These  frames  are 
demodulated  and  decoded  and  sent  to  the  deciphering  unit.  The 
frames  which  follow  are  of  uniform  format  and  contain  the  speech 
parameters.  The  most  important  bits  are  protected  by  an  error 
correcting  code. 

5. 2. 4.1  Signal  Presence  Detection 

Signal  presence  detection  is  accomplished  by  looking  for  the 
Toppler  preamble.  The  signal  presence  detector  consists  of  eight 
filters  and  a threshold  detector.  Four  of  the  filters  are  cen- 
tered on  the  transmitted  tones  (with  enough  bandwidth  to  pass  the 
Dopplered  tone)  and  four  are  centered  to  pick  up  only  noise.  The 
total  energy  at  the  outputs  of  the  "noise"  filters  is  summed  over 
successive  time  instants  and  compared  with  "signal"  energy  similarly 
collected  at  the  outputs  of  the  tone  filters.  Signal  presence  is 
declared  when  the  ratio  of  the  two  sums  ("signal"-to-"noise") 
exceeds  a design  threshold. 
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5. 2. 4. 2 Doppler  Estimation 


The  Doppler  estimator  uses  a three-step  iterative  procedure 
to  generate  an  estimate  of  the  Doppler  offset.  The  four  noise 
filters  of  the  signal  presence  detector  are  discontinued  while 
the  four  tone  filters  remain.  The  rate  of  change  of  the  phase 
of  each  tone  indicates  the  Doppler  offset.  To  minimize  the 
effect  of  noise,  an  energy  weighted  average  over  time  and  the 
tones  is  formed.  In  the  first  stage  the  filters  are  centered  on 
the  transmitted  frequencies,  with  a large  enough  bandwidth  to 
enclose  all  expected  Doppler  offsets.  The  estimate  of  the  first 
stage  is  used  to  compute  new  filter  coefficients.  These  coef- 
ficients specify  filters  whose  center  frequencies  match  those 
measured  in  phase  1,  but  with  a smaller  bandwidth.  These  filters 
are  used  in  phase  2 to  re-estimate  the  Doppler  offset  error. 

This  estimation  causes  the  filter  to  again  be  repositioned  and 
narrowed  for  phase  3.  The  result  obtained  from  phase  3 is  taken 
as  the  final  result  of  the  Doppler  estimation  process.  In  the 
remaining  processing  this  estimate  is  used  to  remove  the  Doppler 
offset  by  multiplying  a rotating  phasor  times  the  complex  samples 
of  the  received  signal. 

5. 2. 4. 3 Frame  Sync  Estimation 

The  frame  sync  preamble  consists  of  a complex  sequence  modu- 
lating a wideband  pulse.  The  frame  sync  estimator  first  filters 
the  incoming  signal  by  a filter  matched  to  the  wideband  pulse. 
This  recovers  the  complex  sequence.  The  incoming  sequence  is 
correlated  with  a local  replica  of  the  transmitted  sequence. 


The  peaks  of  this  correlation  function  indicate  the  positions  of 
the  frame  boundaries  in  the  following  signal.  A threshold  is 
compared  to  the  correlation.  All  values  above  the  threshold  are 
used  to  compute  a magnitude  weighted  average  of  the  correlation 
peak  location.  The  sequence  is  repeated  twice  and  the  two 
averages  obtained  are  combined  to  produce  the  estimate  of  frame 
boundary  position. 

5. 2. 4. 4 Frame-Oriented  Processing 

The  frame-oriented  processing  encompasses  the  phase  refer- 
ence frame,  the  cryptographic  key  and  net  control  frames,  and 
the  data  frames.  Figure  5.4  presents  a flow  chart  of  the  opera- 
tions which  occur  during  each  frame. 

The  input  samples  are  acquired  in  either  an  interrupt-driven 
background  mode  or  are  stored  in  a FIFO  in  the  A/D  system.  These 
are  complex  samples  [s(kT)  + js(kT)]  of  the  input  process  s(t). 

It  has  been  assumed  throughout  this  study  that  the  Hilbert  trans- 
form s(t)  is  done  external  to  the  processor.  The  complex  samples 
have  been  chosen  because  Doppler  correction  can  be  accomplished 
simply  by  multiplying  by  a rotating  complex  phasor  at  the  Doppler 
offset  frequency.  Following  this,  the  64  samples  occupying  the 
integration  time  are  transformed  by  a 64-point  complex  FFT  algor- 
ithm to  measure  the  data  tone  phases. 

If  the  current  frame  is  the  phase  reference  frame  then  the 
phases  are  stored  and  the  processing  for  this  frame  ends.  Other- 
wise, the  measured  tone  phases  are  projected  onto  the  phases 
from  the  last  frame  and  thresholded  by  zero  to  accomplish  the  tone 
demodulation.  As  the  tone  demodulation  is  accomplished  the 

/ 
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Figure  5.4  Demodulator  Frame  Activity  Flow  Chart 


r 


decisions  are  used  to  compare  the  received  phases  to  the  expected 
received  phase.  The  mean  of  the  angle  difference  indicates  an 
error  in  the  Doppler  estimate.  An  amplitude  weighted  average  of 
the  Doppler  estimate  error  is  formed  and  is  used  to  up-date  the 
Doppler  estimate  used  by  the  next  frame. 

The  demodulated  bits  are  then  processed  by  an  error  correct- 
ing decoder.  Two  decoders  are  present  in  the  demodulator,  one 
for  the  cryptographic  key  and  net  control  sequence,  the  other 
for  the  speech  parameter  data  frames.  The  cryptographic  key  is 
given  more  protection  than  the  speech  parameters  because  its  cor- 
rectness is  essential  to  the  intelligibility  of  the  entire  message. 

The  final  task  is  the  slot  sync  tracking  function.  The  DFT 
is  used  to  measure  the  energy  in  a blank  tone  slot  for  both  an 
early  and  a late  integration  gate.  (These  early  and  late  gates 
overlap  both  the  guard  time  and  the  integration  time.)  The  energy 
difference  provides  a measure  of  the  centering  of  the  integration 
time  between  the  adjacent  phase  transition  points.  This  value  is 
thus  averaged  and  used  to  control  the  sample  clock  to  maintain 
the  optimum  position  of  the  integration  time. 

5. 2. 4. 5 End-of -Mess age 

Following  each  data  frame  the  demodulator  must  check  for 
end-of-message.  If  found,  the  demodulator  is  returned  to  the 
signal  presence  detection  mode.  A method  of  accomplishing  this 
is  to  compare  the  energy  in  the  active  tone  slots  of  the  FFT  to 
the  energy  in  the  empty  tone  slots.  A low  ratio  can  be  taken  as 
an  indication  of  the  loss  of  signal.  The  exact  point  at  which 
the  loss  is  declared  will  be  influenced  by  the  fading  character- 
istics of  the  channel  and  the  AGC  characteristics  of  the  radio, 
and  deserves  further  study. 
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5 . 3 Processor  Model 

This  section  presents  the  processor  model  used  in  the  imple- 
mentation study.  The  model  was  chosen  to  be  simple  and  general 
and  was  based  on  the  AMD  2900  micro-processor  family  to  match  the 
Quintrell  processors  under  consideration.  Figure  5.5  presents 
the  block  diagram  of  the  processor  model.  The  AMD  2900  processor 
accepts  input  from  the  D-buss  and  places  its  output  on  the  Y-buss. 

The  D-buss  handles  data  from  the  input  devices,  the  data  memory, 
and  the  fast  multiplier.  The  Y-buss  delivers  data  to  the  output 
devices,  the  data  memory,  and  the  fast  multiplier.  The  micro- 
program sequencer  contains  the  program  memory  and  controls  the 
operation  of  the  entire  system.  The  following  assumptions  have 
been  made: 

1)  Modem  functions  are  micro-programmed 

2)  The  data  memory  contains  both  RAM  and  ROM 

3)  Constants  can  be  supplied  to  the  ALU  from 
the  micro-store 

4)  Overlapped  operation: 

a)  A memory  reference  and  an  ALU  operation 
can  be  hidden  behind  a multiplication 

b)  An  ALU  operation  can  be  hidden  behind 
a memory  reference 

5)  Samples  of  the  signal  and  its  Hilbert  transform 
are  available. 


Flow  charts  of  the  demodulator  algorithms  are  presented  in 

Section  5.4.  The  time  required  to  perform  these  flow  charts  will 

be  expressed  in  terms  of  the  three  quantities  t.TI7,  t , and 

1 ALU  memory 

tmult’  t*le  t:^mes  f°r  t^le  execution  of  an  ALU  operation,  a memory 
operation,  and  a multiply  operation,  respectively.  Each  block  in 


the  flow  chart  is  labelled  with  the  4-tuple  (R;  n.T7I,  n , 

ALU’  memory 


Figure  5.5  Processor  Model 


R = Number  of  times  the  block  is  executed  in  the 
flow  chart 


nALU  = Number  °f  cycles  (shifts,  adds,  branching, 

etc.)  for  one  execution  of  the  block 


n = Number  of  data  memory  cycles  in  one  execu- 

memory  tion  Gf  the  block 


nmult 


Number  of  multiplications  in  one  execution 
of  the  block. 


In  assigning  (R;  nAUJ,  nmemory,  nmult>,  overlapped  operation  is 
accounted  for  by  not  counting  memory  cycles  that  can  be  hidden  in 
a concurrent  multiplication.  The  total  execution  time  is  thus: 

t = tATIT  £ R.n.T1I  . + t £ R.n  . + t . £ R.n  . . 

ALU  ^ i ALU,i  memory  ^ 1 memory, 1 mult  ^ 1 mult,i 

5 . 4 Demodulation  Algorithms 

5.4.1  Signal  Presence  Detection 

Figure  5.6  shows  the  proposed  structure  for  determining  sig- 
nal presence.  The  energy  in  the  frequency  slots  covering  the  four 
Doppler  acquisition  tones  is  compared  to  a threshold  set  by  the 
energy  in  four  empty  frequency  slots.  The  filters  are  described 
by  the  transfer  function 


This  filter  can  be  realized  by  the  block  diagram  of  Fig.  5.7.  An 
estimate  of  the  number  of  operations  needed  to  implement  a single 
filter  is  7 ALU  operations,  6 memory  operations,  and  8 multiply 
operations.  This  assumes  concurrency  of  operation  in  that  a memory 
operation  and  an  ALU  operation  can  be  hidden  behind  a multiply 


I 
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3.6  kHz 
Baseband 


Figure  5.6  Signal  Presence  Detector 


operation,  and  that  an  ALU  operation  can  be  hidden  behind  a memory 
operation.  The  further  assumption  is  made  that  one  sub-routine 
will  perform  all  eight  filters  by  operating  on  successive  offsets 
in  arrays  containing  the  filter  coefficients  and  state  variables. 
The  time  to  execute  the  filters  for  each  sample  value  is 


t = 56  t 


ALU 


+ 48  t 


'memory 


+ 64  t 


mult 


An  estimate  of  the  number  of  operations  needed  to  perform 
the  magnitude  operation  on  a filter  output  is  4 ALU  operations, 
one  memory  operation,  and  2 multiply  operations.  The  number  of 
operations  to  perform  the  eight  magnitude  operations  and  the  two 
summations  is  approximately  40  ALU  operations,  8 memory  opera- 
tions, and  16  multiply  operations.  A simple  threshold  operation 
which  requires  the  threshold  to  be  exceeded  N times  consecutively 
before  declaring  signal  presence  can  be  done  in  about  10  ALU 
operations  and  2 memory  operations.  Thus,  the  total  time  required 
per  8 magnitude  operations  and  one  threshold  operation  is 


t 


50  t 


ALU 


+ 10  t 

memory 


+ 16  t 


mult 


The  filter  operations  must  be  done  for  each  sample  but  an 
M to  1 sample  rate  reduction  is  possible  after  the  filters. 

Thus  the  magnitude  and  threshold  operations  can  be  performed  on 
every  M'th  filter  output  and  the  required  computations  spread 
over  the  time  remaining  after  the  filter  operations  during  the  M 
samples.  The  total  execution  time  per  sample  for  the  detection 
of  signal  presence  is  then 


1 CALU  + 1 

( 48  + —■  \ t 4-  (64  + 

\ M / memory  \ M 1 
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5.4.2  Doppler  Estimation 


I 


The  structure  of  the  Doppler  estimator  is  shown  in  Fig.  5.8. 
The  four  filters  each  isolate  one  of  the  Doppler  reference  tones. 
Each  filter  output  is  mixed  with  a reconstructed  version  of  the 
transmitted  tone  to  produce  a signal  which  contains  only  a tone 
at  the  Doppler  frequency.  From  these  four  signals  an  amplitude 
weighted  average  of  the  rate  of  change  of  the  phase  is  computed. 
This  value  can  be  scaled  to  produce  the  Doppler  frequency.  These 
estimates  can  be  averaged,  if  desired,  to  reduce  the  effects  of 
noise. 

The  processing  to  determine  the  Doppler  frequency  requires 
the  values  of  two  derivatives  per  tone.  These  derivatives  must 
be  approximated  digitally.  An  implementation  in  which  the  deriva- 
tive is  approximated  by 

x(k)  = x(k)  - x(k-l) 

has  been  studied.  The  input  to  the  processor  is  modeled  as 
x(k)  = cos  (kA0) 
y(k)  = sin  (kA0) 

The  output  of  the  processor  is  A0,  the  estimate  of  A0-  The 
values  of  A0  and  A0/A0  have  been  tabulated  for  various  values 
of  A 0 in  Table  5-3. 
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TABLE  5-3 

DOPPLER  ESTIMATION  ERROR 


A 0 

A 

A 0 

A0/A0 

0° 

0° 

(1.000) 

5° 

4.994° 

.999 

o 

o 

1— f 

9.949° 

.995 

15° 

14.829° 

.989 

20° 

19.596° 

.980 

25° 

24.214° 

.969 

30° 

28.648° 

.955 

35° 

32.864° 

.939 

40° 

36.829° 

.921 

45° 

40.514° 

.900 

/N 

A A 

The  values  of  — - 
A A 


are  fit  excellently  by 


L 


A 0 , .1  ^2 

— - 1 2 (A6) 

(45)2 


i •} 

A 0 = A 6 - ' v (A0) 

(45)2 


To  keep  the  estimation  error  small,  the  sampling  rate  at  the 
input  to  the  processor  must  be  high  enough,  relative  to  the  maximum 
expected  Doppler,  so  that  A0  is  in  the  range  of  accurate  estimation. 
This  restricts  the  amount  of  sampling  rate  reduction  that  can 
occur  after  the  multipliers  in  Fig.  5.8.  No  sampling  rate  reduc- 
tion and  a 3600  Hz  sampling  rate  results  in  a .21  Hz  error  in  the 
estimation  of  a 75  Hz  Doppler  frequency. 
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In  implementing  the  Doppler  estimator  of  Fig.  5.8  it  is 
necessary  to  up-date  the  filters  with  each  sample.  The  remainder 
of  the  processing  can  be  done  in  a background  manner  by  storing 
two  successive  outputs  of  each  filter  and  processing  them  in  the 
time  leftover  from  sampling  and  filtering  operations.  The  four 
filters  are  the  same  as  the  upper  four  filters  of  Fig. 5.6  and  can 
be  implemented  in  28  ALU  operations,  24  memory  operations,  and 
32  multiply  operations.  The  time  required  is 


C - 28  (ALU  + 


24  t + 32  t 

memory  mult 


Since  the  derivatives  required  in  the  estimation  process  are 
approximated  by  the  difference  of  two  samples,  only  two  successive 
output  values  from  each  filter  are  necessary  to  produce  an  esti- 
mate of  the  Doppler  frequency.  Also,  since  the  absolute  phase  of 
each  tone  after  the  multiplier  is  of  no  consequence,  only  one  of 
each  pair  of  filter  outputs  needs  to  be  multiplied  by  a complex 
phasor  to  remove  the  tone  center  frequency.  If  the  output  of  the 
i'th  filter  is  represented  as  x^(k)  + jy^(k),  then  the  following 
equations  define  the  Doppler  estimator: 


x^(k)  = x^k)  - x^(k-l)  cos  u)J  + y^(k-l)  sin  u\T 
yi(k)  = y^k)  - x^k-1)  sin  10.T  - yi(k-l)  cos  u^T 


2’V0T 


£ (x1(k)y1(k)  - y1(k)x1(k)) 
i 

L (x^(k)  + y*(k)) 
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An  estimate  of  the  number  of  operations  required  to  implement 
this  is  146  ALU  operations,  28  memory  operations,  and  32  mul- 

I 

tiplier  operations.  The  time  required  to  execute  these  instruc 
tions  is 


t = 146  tA1II  + 28  t + 32  t .. 

ALU  memory  mult 


If  a Doppler  frequency  estimate  is  made  every  P samples, 
the  execution  time  per  sample  for  the  filters  and  the  estimator 
is  given  by 


t ...  . = 

estimate 

l28  + 'p6)  CALU  + (M  + 'pi 

t + 

memory 

32  + t . 

P / mult 

5.4.3  Doppler  Correction 

The  input  to  the  processor  is  samples  of  the  baseband  audio 

signal,  s(t),  and  its  Hilbert  transform  s(t).  The  sum  r(t)  = s(t)  + 

js(t)  represents  the  positive  frequencies  of  s(t).  The  phase-keyed 

tones  present  in  r(t)  contain  a Doppler  offset  which  must  be  removed 

prior  to  demodulation.  This  can  be  accomplished  by  multiplying  by 

-j2irvQt 

the  complex  phasor  e , where  v„  is  the  Doppler  frequency.  The 

-j2wv0t 

signal  rc(t)  = r(t)e  is  used  as  the  received,  Doppler- 

corrected  signal  in  the  frame  sync  interval  and  in  the  demodula- 
tion interval. 

-j2irv0t 

The  complex  phasor  e is  constructed  from  a look-up 

table  containing  64  samples  of  one  cycle  of  a cosine.  The  table 
is  addressed  by  a pointer  representing  the  product  Vgt.  This 


5-29 


pointer  is  incremented  each  sample  interval  by  a quantity  related 

to  the  Doppler  estimate.  Because  the  table  length  is  a power  of 

2,  the  pointer  can  be  maintained  within  proper  bounds  simply  by 

ignoring  the  carry  out  of  the  most  significant  bit.  The  top  six 

bits  of  the  pointer  can  be  used  to  address  an  entry  in  the  table 

and  its  neighbor.  The  remaining  bits  (10  in  a computer  with 

16  bit  words)  are  used  to  perform  a linear  interpolation  between 

-j2irvnt 

the  two  table  values.  This  results  in  the  real  part  of  e u . 

The  process  is  repeated  with  a quarter  cycle  shift  in  the  table 

“j  2ifvQt 

look-up  to  find  the  imaginary  part  of  e . The  remaining 

-j2*v0t 

step  is  the  calculation  of  rc(t)  = r(t)e 

Figure  5.9  shows  a flow  chart  for  the  Doppler  correction 
algorithm.  The  variables  used  and  their  meanings  are  as  follows: 

6 = Equal  to  the  current  angle  of  the  complex 
phasor  used  for  Doppler  correction. 

I 


A 9 = The  sample-to-sample  change  in  Q.  A 6 is 
proportional  to  the  Doppler  frequency. 

ITP  = The  fractional  part  of  6.  ITP  is  used  to 
interpolate  between  values  from  the  cosine 
look-up  table. 

LI  and  L2  = The  indices  for  the  cosine  look-up  table 
which  bracket  the  desired  value. 

-j2irv0t 

A = Real  part  of  e 

-j2TTVQt 

B = Imaginary  part  of  e 
r = Complex  sample  of  the  input. 


r = Doppler  corrected  sample. 


Beside  each  block  of  the  flow  chart  is  a 3-tuple  (n^^,  nmernory’ 
n ) giving,  respectively,  the  number  of  ALU  operations,  the 
number  of  memory  operations,  and  the  number  of  multiply  opera- 
tions to  execute  the  block.  In  the  worst  case  all  blocks  may  be 
executed  for  a particular  sample,  requiring  27  ALU  operations, 

7 memory  operations,  and  8 multiply  operations.  The  processor 
time  to  execute  these  operations  is  the  per  sample  Doppler  cor- 
rection time  t^„: 


27  tAIII  + 7 t + 8 t ... 

ALU  memory  mult 


Doppler  correction  must  take  place  during  the  sync  pre- 
amble and  during  data  transmission.  This  means  that  t^  is  a 
per  sample  overhead  to  the  other  processing  occurring  during 
these  times. 

5.4.4  Frame  Sync  Matched  Filter 

Frame  sync  is  acquired  by  forming  a weight  average  based  on 
the  correlation  of  the  received  signal  with  a local  sequence.  The 
elements  of  this  sequence  are  ± 1 and  + j . The  sequence  is 
1 -1  1 j - j -1  j - j - j 1 j j 1 1 1 . The  correlation  is  formed  by 
passing  the  received  sequence  through  a FIR  filter  whose  taps  are 
the  complex  conjugates  of  the  sequence  terms.  The  algorithm  takes 
advantage  of  the  fact  that  no  multiply  operations  are  necessary 
when  multiplying  by  ± 1 or  ± j . 


I 
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Figure  5.10  shows  a block  diagram  of  the  frame  sync  detector. 
The  matched  filter  is  implemented  by  a FIR  filter  with  N taps. 

The  magnitude  squared  of  the  output  of  the  FIR  filter  is  compared 
to  a threshold.  If  the  threshold  is  exceeded,  the  present  term  is 
used  in  the  centroid  calculation  to  determine  the  time  of  peak 
matched  filter  output.  The  phase  reference  frame  will  start  a 
fixed  number  of  samples  following  the  peak  of  the  matched  filter 
output . 


Figure  5.11  shows  a flow  chart  implementing  this  function. 
This  flow  chart  incorporates  an  ending  criterion  not  mentioned  in 
the  introduction.  The  procedure  ends  and  produces  its  result  when 
a limit  on  the  variable  L is  exceeded.  L indicates  the  number  of 


times  the  threshold  has  not  been  exceeded  following  the  last  time 


the  threshold  was  exceeded.  The  variables  used  and  their  meanings 
are  as  follows: 


k 


L 


Time  index  relative  to  which  location  of  the 
matched  filter  peak  is  found,  k does  not  start 
until  the  threshold  on  the  matched  filter  out- 
put is  first  exceeded 

A negative  value  of  l indicates  that  the 
threshold  has  not  yet  been  exceeded,  a non- 
negative value  of  l indicates  the  number  of 
times  the  threshold  has  not  been  exceeded 
since  the  last  time  it  was  exceeded 


A = Numberator  of  centroid  sum 
B = Denominator  of  centroid  sum 
I = Loop  counter  for  FIR  filter 
F = Filter  output 

r(I)  = Delayed  sample  values  used  in  the  FIR  filter 

M - Magnitude  squared  of  FIR  filter  output 

A = Number  of  samples  between  the  peak  of  the  FIR 
filter  output  and  the  start  of  the  phase 
reference  frame 


T = The  number  of  samples  to  be  skipped  before 
the  start  of  the  phase  reference  frame. 
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The  blocks  in  Fig.  5.11  along  the  longest  per  sample  path  are 

labelled  with  the  4-tuple  (R;  n.TIt,  n , n lA_).  The  number 

K v ALU  memory  mult 

of  operations  required  is  164  ALU  operations,  60  memory  operations, 
and  3 multiply  operations.  The  time  required  per  sample  is 


CMF  164  tALU  + 


60  t + 3 t 

memory  mult 


5.4.5  64-Point  FFT 

Two  algorithms  for  computing  a 64-point  DFT  have  been  studied. 
The  first  is  a radix-2,  decimation-in-frequency  algorithm  and  the 
second  is  a radix-4,  decimation-in-frequency  algorithm. 

5.4. 5.1  Radix-2  FFT 

The  flow  chart  of  Fig.  5.12a  and  Fig.  5.12b  gives  an  algorithm 
for  the  computation  of  a 64-point  DFT  by  the  radix-2,  decimation- 
in-frequency  technique.  The  algorithm  is  in  place  and  accepts 
input  in  normal  order  and  produces  output  in  a bit  reversed  order- 
ing. The  (R;  nATIT,  n , n ...)  assignments  shown  assume  a 

fixed-point  representation  of  the  numbers  and  a memory  assignment 
of  one  word  for  each  real  quantity  and  two  words  for  each  complex 
quantity . 

T 

The  looping  structure  of  this  flow  chart  allows  the  detection 
of  all  instances  of  multiplication  by  1 or  j and  replaces  them 
with  addition  to  save  execution  time.  The  innermost  loop  (Fig. 
5.12b)  successively  performs  adjacent  butterfly  operations.  The 
middle  loop  successively  sequences  through  groups  of  butterflies 
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f ( J2  + K)  = W*fA-B 


f ( J2  + K 


»(63;0,2,0) 


f(J2  + K)  = A-B 


Figure  5.12b  Radix-2  FFT 
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and  the  outer  loop  counts  the  six  iterations  necessary  to  compute 
the  DFT.  This  algorithm  requires  a table  of  15  samples  of  a sine 
wave.  Tabulating  the  values  in  the  flow  chart  results  in  an  execu- 
tion time  of 


t = 1509  t.T1I  + 1634  t + 392  t , „ 

ALU  memory  mult 


5. 4. 5. 2 Radix-4  FFT 

The  flow  chart  of  Fig.  5.13a  and  Fig.  5.13b  gives  an  algorithm 
for  the  computation  of  a 64-point  DFT  by  the  radix-4,  decimation- 
in-  frequency  technique.  The  input  and  output  ordering  and  memory 
allocation  are  the  same  as  for  the  radix-2  algorithm. 

The  looping  structure  of  the  flow  chart  allows  the  replace- 
ment of  all  multiplication  by  1 with  addition.  The  remainder  of 
the  looping  structure  is  the  same  as  for  radix-2,  with  minor  mod- 
ifications due  to  butterfly  size.  This  algorithm  requires  a table 
of  64  samples  of  a sine  wave.  Tabulating  the  values  in  the  flow 
chart  results  in  an  execution  time  of: 


t = 2477  t . _ + 768  t + 324  t - 

ALU  memory  mult 


5.4.6  Data  Demodulation  From  Tone  Phases 

The  4-phase  DPSK  modulation  planned  for  the  ANDVT  is  diagrammed 
in  Fig.  5.14  below.  The  solid  vector  r(k-l)  is  the  tone  phase 
measured  during  frame  k-1.  The  dotted  vectors  indicate  the  pos- 
sible expected  tone  phases  during  frame  k,  conditioned  on  the 
values  of  the  two  data  bits  x^x^.  The  two  bits  can  be  demodulated 
by  examining  the  signs  of  the  dot  and  cross  products  of  r(k)  and 
r(k-l).  Specifically: 
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KW  = K/4  (21 ; 2 , 0 , 0) 


K-l 

f(J  +K) 
f(J  + K) 
f(J  +K) 
f(J  +K) 

A + f ( J + LIT  + K) 

B - jf(J  + LIT  +K) 

C - f(J  + LIT  + K) 

D + j f ( J + LIT  + K) 

A + f(J  + 2*LIT  + K) 

B - f (J  + 2*LIT  + K) 

C + f (J  + 2*LIT  +K) 

D - f (J  + 2*LIT  +K) 

A + f(J  + 3*LIT  + K) 

B + j f ( J + 3*LIT  + K) 
C - f(J  + 3*LIT  +K) 

D - jf  (J  + 3*LIT  + K) 
+ K)  = A 


(48 ; 44 , 10 , 0) 


K = 0 
? 


(48;1,0,0) 
SNV  Yes 


i (27 ;2,6,12) 


KW  = KW  - NIT 
f(J  + LIT  + K)  = 8*W(KW) 
f(J  + 2*LIT  + K)  = C*W(2*KW) 
f(J  + 3*LIT  + K)  = D*W(3*KW) 


f(J  + LIT  + K)  = B 

f(J  + 2*LIT  + K)  = C (21 ; 6 , 6 , 0) 

f(J  + 2*LIT  + K)  = D 
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Figure  5.13b  Radix-4  FFT 
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Figure  5.14  Receiver  Tone  Phases 
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x1  = sgn  (r(k-l)  x r(k)) 
xn  = sgn  (r(k-l)  • r(k)) 


where 


sgn  (x)  = 


x > 0 


x < 0 


Figure  5.15  shows  a simple  flow  chart  for  demodulating  the 

data  from  the  tone  phases.  The  flow  chart  also  incorporates  moving 

the  new  tone  phases  into  the  memory  occupied  by  the  old  tone  phases 

to  be  ready  to  demodulate  the  next  frame.  Each  flow  chart  block 

is  labelled  with  (R;  n.TTI,  n , n - ) . For  demodulation  of 

ALU  memory  mult 

L tone  phases,  the  processor  time  required  is  approximately: 

tnpMm  = (2  + 11L)  tATTT  + 4L  t + 4L  t 

DEMOD  ALU  memory  mult 

5.4.7  Doppler  Tracking  During  Data  Transmission 

One  method  of  tracking  changes  in  the  Doppler  offset  frequency 
is  to  include  an  unkeyed  reference  tone  in  the  tone  library.  The 
tone  is  surrounded  by  frequency  guard  bands  to  allow  the  receiver  to 
separate  it  from  the  data  tones  and  determine  its  Doppler  offset. 

This  method  suffers  two  drawbacks.  First,  a frequency  selective 
fade  could  attenuate  the  reference  tone  and  cause  Doppler  tracking 
to  be  inaccurate  for  the  duration  of  the  fade.  Second,  the  number  of 
tones  proposed  for  the  ANDVT  may  not  leave  sufficient  bandwidth  for 
the  reference  tone  and  its  guard  bands.  A solution  to  the  first  draw- 
back is  to  include  multiple  reference  tones  and  their  guard  bands. 
However,  this  makes  the  second  drawback  even  more  critical. 

Another  method  of  tracking  changes  in  the  Doppler  offset 
frequency  uses  the  information  present  in  the  data  tone  phasors. 

The  frame-to-frame  phase  change  of  a tone  includes  a component 
due  to  the  data  and  a component  proportional  to  the  error  in  the 

5-44 


Doppler  measurement.  The  Doppler  component  of  the  phase  change 
is  the  same  for  all  tones  in  a given  frame  and  causes  the  tone 
phasors  to  precess  from  frame- to- frame.  Two  techniques  of 
extracting  the  Doppler  component  are  discussed  below.  The  first, 
referred  to  as  data  averaging,  measures  the  phase  change  on  each 
tone  and  averages  a sufficient  number  so  that  the  data  component 
in  the  average  is  small.  The  second,  referred  to  as  data  feed- 
back, uses  the  results  of  the  demodulation  process  to  cancel  the 
data  component  of  the  phase  change  and  to  produce  a measurement 
requiring  less  averaging. 

Data  Averaging 

The  frame-to-frame  phase  difference,  at  the  transmitter,  of 
a data  tone  can  be  either  135°,  45°,  -45°,  or  -135°.  If  the 
transmitted  data  is  uniformly  distributed  then  each  of  the 
possible  phase  differences  will  occur  with  probability  .25.  The 
mean  of  the  phase  difference  A 6 is  A0  = 0 and  the  variance  is 
a2  Q = 10125(°)2. 

If  the  receiver  does  not  have  a correct  estimate  of  the 

Doppler  then  the  frame-to-frame  phase  difference  at  the  receiver 

is  = A0  + 6 , where  0^  is  an  angle  proportional  to  the  error 

in  the  receiver  Doppler  estimate.  The  mean  and  variance  of  A«4 

2 2 

are  A 6 = 0„  and  <j  A / = 10125(°)  . If  the  A 6 from  n independent 
D A v 

tones  are  averaged,  the  variance  of  the  average  is  reduced  to 


For  Doppler  frequency  errors  of  1 Hz  or  less,  the  maximum 
value  of  0D  is  8.1°.  If  the  3-sigma  point  of  an  average  over  n 
tones  is  set  to  8.1°,  the  required  value  of  n is  1389.  If  there 
are  43  tones/frame,  then  the  averaging  must  occur  over  33  frames. 


5-46 


Since  the  Doppler  may  change  at  a rate  of  3.5  Hz/sec,  the  actual 
Doppler  error  may  have  changed  by  2.6  Hz  during  the  measurement 
interval.  Approximately  half  of  the  2.6  Hz  change  would  be 
indicated  by  the  averaged  measurement.  It  can  be  seen  that  the 
amount  of  averaging  required  to  produce  what  is  still  a noisy 
estimate  of  the  Doppler  error  has  limited  the  rate  of  change  of 
Doppler  frequency  which  can  be  tracked  by  this  method.  If  the 
3-sigma  point  of  the  average  is  set  to  1°  (a  much  less  noisy 
measurement),  then  n = 91125  or  2120  frames,  an  impractical  amount 
of  averaging. 

Data  Feedback 

The  data  feedback  technique  uses  the  result  of  the  demodu- 
lation process  to  remove  the  data  component  of  the  tone  phase 
change.  Averaging  is  still  desirable,  however,  to  remove  the 
effects  of  additive  noise,  round-off  noise,  tone  fading,  and 
demodulation  error. 

Figure  5.16  shows  the  four  possible  phasors  received  during 
frame  k with  an  error  in  the  Doppler  estimate,  relative  to  the 
received  phasor  from  frame  k-1.  The  data  feedback  technique 
uses  the  demodulated  data  to  rotate  the  phasor  r(k  - 1)  to  the 
position  expected  for  r(k)  if  there  were  no  Doppler  error.  The 
angle  0^  is  then  measured  and  an  amplitude  weighted  average  com- 
puted over  the  set  of  data  tones.  If  we  define: 


A.  = Re  {r.(k)r*(k-  1)} 


B.  = Im  {r..(k)r*(k-l)} 


C.  = r.(k)r*(k) 
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a.  = 
1 


Ai  > 0 

A±  < 0 


B.  > 0 

l — 


B.  < 0 

l 


Then  if  is  the  frame  time,  the  Doppler  error  is 


A v = 


E (a.Bi  - biAi] 

1 ^ 
l 


It  should  be  noted  that  the  derivation  of  (6)  involved  the 
approximation  of  a derivative  by 


x(k) 


x(k)  - x(k  - 1) 
T 


An  interesting  feature  of  the  estimator  defined  in  (6)  is 
that  the  A^  and  B^  are  the  decision  variables  used  to  demodulate 
the  data.  The  notation  used  here  differs  from  that  of  Section 
5.4.6  so  the  comparisons  for  the  demodulator  are  restated  as: 
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xli = 


A.  > 0 
1 — 


A.  < 0 

l 


B.  > 0 

l — 


Bi  < 0 


Since  the  quantities  in  the  numerator  of  (6)  are  computed  in  the 
course  of  demodulation,  the  additional  computations  required  to 
complete  the  estimator  are  simply  those  necessary  to  perform  the 
summations,  compute  the  terms  in  the  denominator,  and  perform 
the  division.  An  estimate  of  this  is  (104  + 8L)  ALU  operations, 

4 memory  operations,  and  2L  multiply  operations  for  L data  tones. 
The  amount  of  processor  time  in  excess  of  the  demodulation  time 


t.  , = (104  + 8L)  t.T„  + 4t  +2Lt 

track  ALU  memory  mult 
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5.4.8  Slot  Sync 

I j 

To  maintain  a proper  frame  sync  position  in  the  presence 

of  clock  drift  and  changes  in  the  multipath  environment,  it  is 
necessary  to  track  the  position  of  the  phase  keying  and  adjust 
the  receiver  clock  to  keep  it  within  the  guard  time.  This 
serves  to  ensure  that  the  keyed  tones  are  orthogonal  during  the 
integration  time.  One  method  of  accomplishing  this  is  slot 
sync.  If  an  integration  gate  covers  a phase  keying  point,  the 
phase  discontinuity  causes  crosstalk  between  the  tones.  The 
slot  sync  technique  measures  the  energy  in  an  empty  frequency 
slot  of  an  early  integration  time  and  also  measures  the  energy 
in  the  same  frequency  slot  of  a late  integration  time.  The 
difference  of  these  energies  provides  a discriminator  function 
to  indicate  the  centering  of  the  data  integration  time  within 
the  frame  time. 

This  technique  requires  the  measurement  of  a tone  amplitude 
in  each  of  two  integration  times.  The  FFT  is  not  efficient  in 
this  application  because  the  amplitude  of  only  one  tone  is 
desired  per  set  of  samples.  An  efficient  method  of  obtaining  a 
small  number  of  tone  amplitudes  is  the  Goertzel  algorithm  for 
the  DFT.  The  Goertzel  algorithm  has  two  advantages  over  the 
direct  form  of  the  DFT.  First,  the  number  of  real  multiply  oper- 
ations is  half  that  of  the  direct  DFT,  and,  second,  only  a single 
coefficient  is  required  rather  than  the  N required  by  the  direct 
DFT.  This  saves  table  look-up  time  during  program  execution. 

I 
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Algorithm  Description 


The  Goertzel  algorithm  is  a digital  filter  approach  to 
obtaining  the  tone  phasor.  Figure  5.17  shows  a block  diagram 
of  the  filter  used.  The  N samples  are  input  to  the  filter 
and  the  output  at  t = NT  is  y(NT)  = F(kfi),  the  value  of  the  DFT 
for  the  k'th  tone.  Since  the  output  of  the  filter  is  needed 
only  at  n = N,  the  operations  to  the  right  of  the  dotted  line 
need  only  be  performed  for  that  sample  time. 

Figure  5.18  shows  a flow  chart  implementing  the  Goertzel 
algorithm  with  N = 64.  The  blocks  are  labelled  with  (R;  n , 

ALU 

n , n . ).  The  number  of  operations  to  execute  the  flow 

chart  once  is  646  ALU  operations  and  132  multiply  operations. 
This  must  be  done  twice  for  slot  sync  and  a few  additional  fin- 
ishing operations  to  implement  the  discriminator,  approximately 
4 memory  operations  and  6 multiply  operations.  The  total  time 
for  the  slot  sync  operation  is  given  approximately  by: 


'slot 


= 1292  t 


ALU 


+ 4 t 


memory 


+ 270  t 


mult 


5.4.9  Golay  Soft-Decoder 

This  section  describes  an  implementation  of  a soft-decision 
decoder  for  the  (24,12)  Golay  code  with  overall  parity.  This 
decoder  selects  the  four  least  reliable  bits  in  the  code  word 
and  performs  a hard-decode  operation  for  each  of  the  16  trial 
error  patterns  generated  by  these  four  bits.  The  decoded  word 
with  the  lowest  analog  weight  is  selected  as  the  output  of  the 
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decoder.  The  absolute  values  of  the  decision  statistics 
defined  in  the  first  two  equations  of  Section  5.4.7  are  used 
in  the  selection  of  the  least  reliable  bits  and  in  the  computa- 
tion of  the  analog  weights. 

The  hard-decoding  algorithm  chosen  is  a table  lookup 
technique.  The  syndrome  addresses  a 1-bit  wide  table  whose 
content  indicates  the  presence  or  absence  of  an  error  in  the 
rightmost  bit  position.  The  syndrome  is  successively  modified 
to  correspond  to  a 1-bit  cyclic  shift  of  the  code  word  and  the 
table  indicates  the  presence  or  absence  of  an  error.  This  is 
performed  23  times  to  check  the  validity  of  the  entire  Golay 
code  word.  The  parity  of  the  decoded  word  is  used  to  check  the 
validity  of  the  overall  parity  bit.  The  modified  error-trapping 
procedure  of  Kasami  was  compared  to  the  table  lookup  method  and 
found  to  require  either  a much  larger  table  or  to  require  sig- 
nificantly more  processor  time. 

Figure  5.19  shows  a flow  chart  of  the  soft-decoding  pro- 
cess. Four  arrays  of  24  words  and  two  lookup  tables  are  used 
in  this  implementation.  The  four  arrays  are  the  received  bit 
array,  the  weight  array,  the  error  pattern  array,  and  the  old 
error  pattern  array.  The  received  bit  array  and  the  weight 
array  are  set  by  the  demodulation  routines.  The  error  pattern 
array  holds  the  error  pattern  when  each  of  the  trial  error  pat- 
terns is  decoded.  The  old  error  pattern  array  holds  the  error 
pattern  which  resulted  in  the  lowest  weight.  One  of  the  lookup 
tables  is  used  by  the  hard  decoder  to  determine  the  presence  of 
an  error  in  the  rightmost  bit  position.  The  remaining  lookup 
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Figure  5.19  Soft  Decoding  Flow  Chart 
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table  contains  the  syndrome  modifications  used  to  compute  the 
starting  syndrome  from  the  received  syndrome  given  the  trial 
error  pattern.  ) 

The  flow  chart  starts  by  locating  the  four  least  reliable 
bits,  computing  the  syndrome  of  the  received  word,  and  setting 
the  variable  which  holds  the  weight  of  the  best  error  pattern 
(OWT)  to  a maximum  (insuring  that  the  first  hard  decode  is  stored 
in  the  best  error  pattern  array  (OEA)).  The  loop  counter  N is 
also  set  to  zero.  Each  pass  through  the  loop  causes  a new  trial 
error  pattern  to  be  generated,  and  a hard  decode  performed  using 
it.  The  weight  (WT)  of  the  error  pattern  is  compared  to  OWT  to 
determine  if  the  present  contents  of  the  EA  should  replace  the 
contents  of  OEA.  Following  the  sixteenth  pass  through  the  loop, 
the  decoded  word  is  taken  as  the  bit-by-bit  exclusive-OR  of  the 
received  word  with  OEA. 

Each  block  in  Fig.  5.19  is  labelled  to  indicate  the  number 
of  times  the  block  is  executed,  the  number  of  ALU  operations,  the 
number  of  memory  operations,  and  the  number  of  multiply  opera- 
tions, respectively.  Tabulating  these  numbers  results  in  a count 
of  13456  ALU  operations  and  2580  memory  operations.  The  time 
required  to  execute  this  operation  is 


CSD  " 13456  6 ALU  + 2580  6 


memory 


Various  stopping  rules  may  be  employed  to  speed  up  the  above 
processing  time.  These  stopping  rules,  which  yield  no  loss  in  per- 
formance, may  be  used  in  the  actual  implementation,  if  required. 


- 


I 
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5. A. 10  BCH  Decoder  for  KG  Sequence 


The  cryptographic  key  and  net  control  sequences  are  protected 
by  a (249,101;45)  BCH  code.*  "'’he  single  code  word  has  a total  of 
249  bits,  101  of  which  are  information  bits  containing  the  crypto- 
graphic key  and  net  control  sequences.  The  minimum  distance  of 
the  code  is  45,  allowing  22  or  fewer  errors  to  be  corrected.  A 
flow  chart  for  the  decoding  of  this  code  is  given  in  Fig.  5.20. 

The  decoding  operation  consists  of  three  main  steps:  syndrome 
calculation,  calculation  of  the  error  location  polynomial  cr(X),  and 
the  Chien  search  to  locate  the  errors.  If  the  order  of  a(X)  is 
< 22  and  the  number  of  errors  found  in  the  Chien  search  equals  the 
order  of  cr(X)  then  the  decoded  output  V(X)  is  the  modulo-2  sum 
of  the  received  sequence  R(X)  and  the  error  sequence  E(X)  from  the 
Chien  search.  Otherwise  we  take  R(X)  as  the  decoded  output. 


The  decoding  process  begins 
dromes  S^,  i = 1,2,3,  ...,  44. 
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S.  = S r a 
j=0  J 


ij 


with  the  computation  of  the  syn- 
The  odd  syndromes  are  computed  as: 


The  even  syndromes  can  be  determined  more  easily  by  the  relation: 


S 


2i 


(S.) 


2 


i i 8 

The  terms  a J are  field  elements  in  the  Galois  field  GF(2  ). 

A flow  chart  for  the  calculation  of  the  syndromes  is  given 
in  Figs.  5.21a  and  5.21b.  This  algorithm  uses  a table  look-up  to 
find  the  a1"'  used  in  the  summation  for  the  odd  syndromes.  Figure 


The  complexity  of  this  code  is  equivalent  to  that  of  the 
(250,100;46)  code  considered  in  Section  3,  and  the  performance 
is  essentially  equivalent.  The  above  code  may  be  selected  in 
place  of  the  (250,100;46)  code. 
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Figure  5.20  BCH  Decoding  Algorithm 


Figure  5.21b  Syndrome  Computation 
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5.21b  is  the  part  which  stores  the  syndromes  in  memory  and  com- 
putes the  required  even  syndromes.  Each  block  is  labelled  with 

(R;  nATII,  n , n ,.).  Tabulating  these  numbers  results  in  a 
ALU’  mem  mult  ° 

count  of  33737  ALU  operations  and  11044  memory  operations  for  the 
syndrome  calculations. 

The  second  stage  of  the  decoding  procedure  is  the  calculation 
of  the  error  locating  polynomial  ct(X).  The  roots  of  a(X)  give 
the  locations  of  the  errors.  a(X)  is  found  by  an  iterative  pro- 
cedure using  the  syndromes.  A flow  chart  of  this  iterative  pro- 
cedure is  given  in  Fig.  5.22.  This  implementation  assumes  a 
specific  memory  layout  for  the  polynomials  a(X),  cr'(X),  and  t(X). 
For  ct(X)  this  is:  0 1 ^ ...  in  consecutive  memory 

locations.  For  a ' (X)  this  is  cr|  a2  CT3  ***  a43  in  consecuti-ve 
memory  locations.  For  t(X)  this  is  0 0 Tq  T2  •••  i-n 

consecutive  memory  locations.  These  layouts  are  necessary  to  keep 
the  loop  structures  simple  within  the  various  blocks  of  Fig.  5.22. 

The  boxes  in  Fig.  5.22  are  labelled  with  the  3-tuple  (n  , 

/\1_»  u 

n , n , where  each  is  the  total  number  of  ALU  operations, 
mem’  mult  r ’ 

memory  operations,  or  multiplication  operations,  respectively, 
that  is  required  in  the  box  for  the  entire  execution  of  the  flow 
chart.  This  has  been  done  because  the  amount  of  effort  in  many 
boxes  varies  with  the  index  k and  the  usual  notation  is  inapprop- 
riate. Summing  all  of  these  entries  results  in  19374  ALU  opera- 
tions and  8343  memory  operations  for  the  determination  of  the 
error  locating  polynomial. 

The  last  significant  amount  of  computation  in  the  error  cor- 
rection procedure  is  the  Chien  search.  The  Chien  search  com- 
putes the  value  of  a(X)  for  each  a1  corresponding  to  a bit  in  the 
received  code  word,  cj  (a1)  = 0 indicates  that  an  error  has  been 


Enter 


Figure  5.22  Iterative  Procedure  for  a(X) 


found.  The  flow  chart  for  this  operation  is  given  in  Figs.  5.23a 
and  5.23b.  The  loop  of  Fig.  5.23a  computes  the  log  of  each  coef- 
ficient of  a (X)  and  performs  a correction  due  to  the  fact  that,  the 
(249,101 ;45)  code  is  a shortened  (255,107 ; 45 ) code.  The  correc- 
tion simply  serves  to  advance  the  Chien  search  past  the  missing 
bits.  Figure  5.23b  performs  the  actual  search  by  advancing  each 
term  of  the  polynomial  by  the  correct  power  of  a and  summing  the 
non-zero  terms.  If  the  term  £ = a (a1)  is  zero  then  a 1 is  entered 
in  the  error  array  and  the  error  count  is  incremented,  otherwise 
a 0 is  entered  in  the  error  array.  The  blocks  in  the  flow  chart 

are  labelled  with  (R:  n.TII,  n , n - ),  tabulating  these 

ALU  mem  mult 

results  in  51726  ALU  operations  and  16749  memory  operations. 

The  total  of  the  operation  counts  for  the  three  major  por- 
tions of  the  decoding  is  104837  ALU  operations  and  36136  memory 
operations . 

5 . 5 Demodulator  Processor  Time  Summary 

Table  5-4  presents  a summary  of  the  study  of  the  implemen- 
tation of  the  demodulator  on  a programmable  processor.  The 
tasks  listed  are  the  algorithms  described  in  the  subsections  of 
Section  5.4.  The  formulas  given  for  signal  presence  detection 
and  for  Doppler  estimation  include  all  operations.  To  find  the 
processor  time  required  per  sample  for  frame  sync  it  is  necessary 
to  add  the  times  for  Doppler  correction  and  for  the  frame  sync 
mathched  filter.  The  time  for  complete  data  demodulation  is 
found  by  summing  the  total  Doppler  correction  time,  the  64-point 
FFT  time,  the  demodulation  time,  the  Doppler  tracking  time  and 
the  slot  sync  time.  The  total  Doppler  correction  time  is  the 
per  sample  Doppler  correction  time  times  the  number  of  samples 
in  a frame.  If  two  Golay  code  words  are  included  then  twice 
the  Golay  soft-decode  time  must  be  added. 
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k—  22 
•l—  5 


(1;4,0,0) 


X^-  log  o23~^\  (22 ; 1 , 2 , 0) 

^^\(22;  1,0,0) 

= log  o\i“ . 


t (22 ; 3, 0 , 0) 


(X-l)  mod 


°23-k  X 
1—  1 + 5 
k—  k - 1 


(22 ; 2 , 1 , 0) 


^ \(22;1,0,0) 

k = 0 Ye.s 


k-*-  249 
c—  0 


(1 ; 3 , 0 , 0) 


Figure  5.23a  Chien  Search 
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The  memory  requirements  for  the  modem  have  been  examined 
briefly.  A data  memory  of  .5  K of  ROM  and  1 K of  RAM  is  suf- 
ficient. For  the  processor  model  of  Section  5.3  a 2 K micro- 
program memory  will  be  adequate. 

Table  5-5  lists  the  assumptions  made  and  presents  an  example 
of  the  use  of  the  formulas  in  Table  5-4.  The  assumptions  relat- 
ing to  the  processor  time  are  believed  to  represent  the  TRW 
Quintrell  processor.  This  table  indicates  that  all  demodulator 
functions  can  be  accomplished  within  the  speed  constraints  of  the 
TRW  processor. 


Ki 
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TABLE  5-5 


PROCESSOR  UTILIZATION  EXAMPLE 


Assume : 


• 3600  Hz  Sampling 

• 39  Data  Tones 

• ‘aLU  - 350  ”s 

• t = 350  ns 

memory 

• t A.=  700  ns 

mult 

• M = P = 16 


Signal  Presence  Detection: 


8.3  x 10  Sec/San.ple 
307,  Processor  Utilization 


Doppler  Estimation: 


4.6  x 10  Sec/Sample 
167,  Processor  Utilization 


Frame  Sync:  (Doppler 
Correction  and  Frame  Sync 
Matched  Filter) 


9.8  x 10  Sec/Sample 
357,  Processor  Utilization 


Data  Demodulation:  (Doppler 
Correction  + 64-Point  FFT  + 
Demodulation  + Doppler 
Tracking  + Slot  Sync) 


4.0  x 10  Sec/Frame 
187,  Processor  Utilization 


Data  Demodulation  and  Decoding: 
(Data  Demodulation  + 2 Golay 
Soft  Decoders) 


1.5  x 10  Sec/Frame 
677,  Processor  Utilization 


BCH  Decoder  for  KG  Sequence 


5.0  x 10  ^ Sec 
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SECTION  6 


ADDITIONAL  MODEM  DESIGN  CONSIDERATIONS 

i 

6 . 1 Introduction 

In  this  section  several  topics  which  are  related  to  the  ANDVT 
modem  are  discussed. 

Section  6.2  outlines  design  considerations  for  optimizing 
the  overall  waveform,  with  particular  emphasis  on  minimizing  the 
effects  of  the  high  peak-to-average  power  ratio  present  in  mul- 
tiple tone  HF  modems. 

Considerations  for  providing  the  ANDVT  modem  with  a data 
transmission  mode  of  operation  are  discussed  in  Section  6.2. 

Finally,  in  Section  6.3  the  application  of  coding  for  perform 
ance  monitoring  is  discussed. 

6 . 2 Signal  Level  Optimization 

For  an  N-tone  HF  modem  with  phase  shift  keying  the  trans- 
mitted signal  is  given  by 

N 

x(t)  = L A cos  (co.t  + 0.)  (6.1) 

i=l 

where  0^  represents  the  information  and  ox  represents  the  fre- 
quency of  the  i'th  tone  which  is  assumed  orthogonal  over  an 
integration  period  of  T second.  The  peak  value  of  this  signal  is 
NA  and  the  rms  value  of  VN  A.  Thins  the  peak-to-rms  ratio  in  terms 
of  N is  given  by 
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HO* 


— - 


peak-to-rms  ratio  = 20  Log  dB 


(6.2) 


which  for  a 40  tone  modem  is  16  dB.  If  one  designs  the  transmitted 
signal  to  be  distortion  free,  regardless  of  the  particular  phase 
relationships  between  the  individual  tones,  the  transmitter  must 

I 

operate  at  a mean  value  which  is  16  dB  below  its  peak  value.  Clearly, 
this  is  an  undesirable  mode  of  operation. 

The  problem  is  further  compounded  by  the  preamble  design  (see 
Section  2)  which  allows  for  a peak-to-rms  ratio  of  less  than  5 dB. 
Note  that  during  signal  presence  detection  and  Doppler  estimation, 
the  four  tone  preamble  can  be  designed  with  a crest  factor  of  1.77 
(2.48  dB)  using  180°  phase  increments  and  using  optimum  phase 
adjustemtns  a crest  factor  of  1.52  (1.84  dB)  can  be  obtained.  Dur- 
ing frame  estimation,  using  raised  cosine  pulses,  a peak-to-rms 
value  of  2.67  (4.26  dB)  is  obtained.  If  cosine  shaped  pulses  are 
used,  a peak-to-rms  value  of  only  3 dB  results.  It  is  desirable  to 
introduce  some  intentional  clipping  in  the  transmitted  signal  during 
the  occurrence  of  the  multiple  tone  signal  to  maximize  the  trans- 
mitted power.  (Note  that  HF  radios  with  some  peak  power  limitations 
are  assumed.)  In  conventional  modems  the  received  signal  is  inten- 
tionally clipped  to  maximize  the  dynamic  range  and  thus,  the  clip- 
ping of  the  transmicted  signal  should  be  somewhat  less  than  that 
required  by  the  receiver. 

A heuristic  discussion  of  these  issues,  based  on  a combination 
of  theoretical  arguments  and  experimental  results,  is  given  in  the 
Codem  I final  report  [6.1].  Following  these  arguments  one  may  con- 
clude that  several  dB  of  intentional  clipping  at  the  transmitter 
is  desirable.  Since  the  preamble  can  operate  at  a peak-to-rms 
ratio  well  below  5 dB,  it  would  be  desirable  to  operate  the  40  tone 
modem  with  at  least  11  dB  of  intentional  clipping  at  the  trans- 
mitted signal.  The  actual  level  of  acceptable  clipping  (which  is 
expected  to  be  somewhat  less  than  11  dB)  can  be  determined 
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experimentally  during  the  expected  follow-on  to  this  present 
contract.  The  clipping  nonlinearity  at  the  transmitter,  and  the 

I 

effects  of  this  nonlinearity  when  digitized  speech  is  transmitted 
over  a HF  channel  model,  make  this  problem  quite  intractable  math- 
ematically. Fortunately  this  problem  is  tractable  by  detailed 
simulations  which  include  the  entire  HF  modem,  the  HF  channel 
model,  and  the  speech  digitizer. 


6 . 3 A Data  Transmission  Mode  of  Operation  for  ANDVT 

It  is  desirable  to  have  a data  transmission  mode  for  data 
rates  of  2400,  1200,  and  600  b/s,  in  addition  to  the  speech 
transmission  mode  of  operation.  While  a frame  rate  of  44.44  f/s 
(frame  of  22.5ms  duration)  is  desirable  when  speech  is  encoded  at 
this  frame  rate,  the  question  arrises  as  to  what  frame  rate  should 
be  used  for  data  transmission.  Conventional  Tadil-A  modems  have 
a 16  tone  four-phase  DPSK  format,  and  a frame  rate  of  75  f/s. 

Thus  a data  rate  of  2400  b/s  is  achieved  before  coding  is  attempted. 
If  a 40-tone  four-phase  DPSK  modem  format  and  a frame  rate  of 
44.44  fps  is  used,  a maximum  data  rate  of  3555.55  b/s  can  be 
achieved.  This  format  enables  us  to  consider  coding  even  at  a 
data  rate  of  2400  b/s.  The  possibility  exists  that  once  a code 
is  chosen  for  the  digitized  speech  signal  the  same  code  may  be 
modified  for  data  transmission  applications.  However,  at  this 
point  a final  decision  on  the  code  for  ANDVT  has  not  been  made. 

The  remainder  of  this  section  presents  preliminary  trade-off 
results  for  convolutional  codes  and  the  (24,12)  block  codes  used 
with  four-phase  DPSK  modems  at  frame  rates  of  75  f/s  and  44.44  f/s. 


• J 


J 
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Hard  and  soft  decoding  is  considered  for  both  the  block  codes 
and  convolutional  codes.  The  soft  decision  performance  is  obtained 
by  a somewhat  optimistic  estimate  obtained  by  assuming  a factor  of 
two  increase  in  the  number  of  errors  that  can  be  corrected  [6.2]. 
Since  decoding  delay  is  constrained  only  by  complexity,  code  bits 
can  be  made  independent  with  sufficient  interleaving,  which  will 
be  assumed.  Figure  6.1  illustrates  the  benefits  of  coding,  in  par- 
ticular soft  decoding,  for  a four-phase  DESK  modulation  format.  The 
Doppler  variations  are  assumed  negligible  on  this  figure  and  thus 
the  performance  is  the  same  for  75  f/s  and  44.44  f/s  modems.  Note 
that  the  bit  error  probability  is  illustrated  in  terms  of  the 
signal-to-noise  ratio  per  tone  (E  /Nq)  which  is  closely  related  to 
the  signal-to-noise  ratio  per  bit  ( E^/Nq) , (see  Appendix  C) . In 
terms  of  E^/Nq  the  44.44  f/s  modem  has  a slight  advantage  (1.67  dB), 
since  the  guard  time  ratio  is  reduced  and  no  extra  tone  is  used  for 
Doppler  tracking. 

Figures  6.2  through  6.6  illustrate  the  effects  of  Doppler 
variations  for  a two-sided  ms  Doppler  spread  B = 2a  of  1,  2,  3, 

4,  and  5 Hz,  when  a frame  rate  of  44.44  f/s  is  utilized.  While 
Doppler  spreads  are  typically  below  2 Hz,  the  importance  of 
coding  for  operating  in  the  presence  of  Doppler  is  evident  from 
these  figures. 

Figures  6.7  through  6.11  illustrate  the  same  curves  when  a 
modem  with  a frame  rate  of  75  f/s  is  assumed.  This  modem  is  less 
sensitive  to  Doppler  spreads.  However,  when  soft  decoding  is 
used  and  B < 2 Hz,  the  performance  difference  between  a 75  f/s  and 
44.44  f/s  modem  is  negligible. 
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Figure  6-1  Performance  of  Rate-1/2  Codes  Over  a Slowly  Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation 
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Figure  6.2  Performance  of  Rate-1/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 1 Hz,  Frame 
Rate  = 44.44  Frames/Sec.) 
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Figure  6.3  Performance  of  Rate-1/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 2 Hz,  Frame 
Rate  = 44.44  Frames/Sec.) 
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Figure  6.5  Performance  of  Rate-1/2  Codes  Over  a Time- Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 4 Hz,  Frame 
Rate  = 44.44  Frames/Sec.) 
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Figure  6.6  Performance  of  Rate-1/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 5 Hz,  Frame 
Rate  = 44.44  Frames/Sec.) 
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Figure  6.7  Performance  of  Rate-l/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 1 Hz,  Frame 
Rate  = 75  Frames/Sec.) 
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Figure  6.8  Performance  of  Rate-1/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 2 Hz,  Frame 
Rate  = 75  Frames/Sec.) 
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Figure  6.9  Performance  of  Rate-l/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 3 Hz,  Frame 
Rate  = 75  Frames/Sec.) 
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Figure  6.10  Performance  of  Rate-1/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 4 Hz,  Frame 
Rate  = 75  Frames /Sec.) 
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Figure  6.11  Performance  of  Rate-l/2  Codes  Over  a Time-Varying  Rayleigh 
Channel  With  4-Phase  DPSK  Modulation  (B  = 5 Hz,  Frame 
Rate  = 75  Frames/Sec 
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Since  the  major  disadvantage  of  using  a 44.44  f/s  modem  is 
its  increased  sensitivity  to  Doppler  spreads  and  since  this  dis- 
advantage is  overcome  with  soft  decoding,  it  is  anticipated  that 
the  data  mode  for  ANDVT  will  also  use  a frame  rate  of  44.44  f/s. 
Further  studies  in  this  direction  should  be  undertaken  after  the 
coding  for  the  speech  digitization  is  finalized.  However,  pre- 
liminary results  indicate  that  if  sufficient  interleaving  is  used 
the  K = 7 convolutional  code  is  superior  to  the  Golay  code.  For 
lower  data  rates  this  convolutional  code  can  readily  be  modified 
[6.2].  Comparisons  between  this  code  and  a multiple  rate  struc- 
ture is  of  interest,  especially  if  a multiple  rate  code  is  selected 
for  speech  transmissions. 

6.4  Performance  Monitoring 

When  error  correction  coding  is  not  used  in  a system,  per- 
formance monitoring  is  achieved  by  observation  of  certain  modu- 
lator features  of  the  received  signal.  As  an  example,  for  a four- 
phase  DPSK  format,  performance  monitoring  can  be  achieved  by  not- 
noting  that  the  receiver  phase  of  each  tone  is  likely  to  fall  in  a 
window  centered  about  the  four  possible  phases  transmitted.  An 
even  simpler  example  is  the  measurement  of  the  ratio  of  "in-band" 
received  power  to  "out-of-band"  noise  and  setting  a threshold  for 
a loss  of  signal  based  on  this  ratio. 

When  coding  is  utilized,  Fig.  6.12  illustrates  a quite  simple 
approach  for  performance  monitoring  which  is  accurate  enough  to 

even  predict  the  channel  (raw)  error  rate  when  this  number  is 

-2 

below  10  . Under  these  conditions,  the  decoder  data  is  highly 

reliable  and  thus  by  comparing  the  raw  data  to  the  re-encoded  data 


. 


Figure  6.12  Performance  Monitoring  Based  on  Error 
Correction  Coding 
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an  accurate  estimate  of  the  channel  error  rate  is  possible.  For 

_2 

poor  channel  error  rates,  i.e.,  above  10  , the  decoded  data  may 

be  unreliable  and  the  channel  error  rate  estimate  may  also  be 
unreliable  and,  in  fact,  indicate  a higher  error  rate  than  is 
actually  present.  This  feature  is  still  quite  acceptable  for  most 
applications . 

Other  more  sophisticated  performance  monitors  are  obtainable 
by  including  channel  measurement  (soft  decision)  information  in 
the  performance  monitoring  algorithm,  such  as  the  use  of  the 
analog  weight  of  the  estimated  error  pattern  [6.3].  Techniques 
of  this  nature  can  be  applied  if  the  concept  illustrated  on  Fig. 
6.12  proves  to  be  inadequate  for  a specific  application. 
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APPENDIX  A 


HF  CHANNEL  MODEL 


Ideally,  the  response  of  one  path  of  the  HF  channel  to  a 
sinusoidal  signal  of  the  form 

x(t)  = A cos(tot-ffl)  (1) 


is 


y(t)  = A r ^ cos  (tot  + fl  + (Pj  ) + n(t)  (2) 

where  r.  is  the  time-varying  amplitude  of  the  path  (Rayleigh 
distributed),  tp  . is  a phase  due  to  Doppler,  time  delay,  etc., 
and  n(t)  is  a complex  white  Gaussian  noise. 

Assuming  a transmitted  signal. 


x(t)  = J2F  cos(cot  + 0) 


(3) 


the  received  signal  over  a multipath  channel  is 

N 

y(t)  = JlV  r . cos  (wt  + G + <p  . ) + n(t) 

j = l J J 


(4) 


,,  r 2- 


where  the  r.  are  independent  Rayleighs,  and  E|  r.  =P.,  the  power 
t h ^ ~ . 1-  J 

of  tile  j path  and  P.  = 1;  and  n(t)  is  the  total  noise  on  all 

i=l  J 

paths.  J 


A-l 


quadrature;  matched  filters 


at  the  receiver 


The  in-phase  and 
form  the  quantities 

T 

r = — f Jl  y (t)  cos  cut  dt 


(5) 


T 


rg  - - — f Jl  y(t)  sin  cot  dt 


(6) 


Tims 


T 


= ~ V2P  Jl  £ r f 

j = i J0 


, N 

± 2j¥  'y  r.  — cos(0+(p.)  + n 

1 j^i  J 2 j < 


coscot  cos  (cut  + 0 + ep  . )dt  + n, 

J ® 


and 


JP 


cos  (0  +cp  . ) 
J 


+ n 


r = /L  */P  r . sin(0+cp.)  + n 
s jT[  J J s 


Note  that 


(7) 


(8) 


and 


P 


1£ 


N, 


21  _ 0 

-nsJ  2T 


A- 2 


(9) 


(10) 


Thus  y , the  signal-to-noise  ratio  per  tone,  is  given  by: 


/N  \ N N 

J_o)  0 0 

\ 2T  / 


For  an  array  of  tones, 


N N 
TONE  PATH 


y(t)  = ^ ]>_,  y2P  r.  cos(u.  t 

i=i  j = l J L 


+ 0 . + n.(t) 

1 ij  lv 


Note  that  for  this  model,  r.  is  assumed  constant  over  all 

J 

frequencies  of  interest  while  <p.  . is  frequency-dependent,  r^  and 

r have  the  same  form  for  each  tone  with  the  r.'s  the  same  in 
s . J 

each  case. 

The  simplicity  of  Eq.  (12)  suggests  that  the  HF  channel  model 

is  well-suited  for  computer  program  implementation.  Define 
N 

|A0.|  tQ  tjie  set  q£  information-bearing  phase  differences 

1 i=l 

for  the  tone  array.  The  quantity, 


0n(k)  " £ ^n(j)  (13) 

J = 1 

which  is  the  actual  transmitted  phase  angle  at  the  k iteration, 
is  formed. 

The  received  phase  angle  for  the  nth  tone  over  the  jth  path 
t h 

at  the  k iteration,  i p (k) , is 

n 


i,.n(U)  - 0n(k)  + tan  + 2»[(W-TJ)kJ  - (>' f sc|1  + f0>T j. 


where  X. (k)  and  Y.(k)  are  the  real  and  imaginary  channel  gains 
J th  ^ th 

over  the  j path  at  the  k iteration;  r.  is  the  time  delay  of 
th 

the  j path,  v.  is  the  Doppler  offset,  (n  f + f^)  is  the  dis- 
J ^ r p 

tance  of  the  nth  tone  from  the  carrier,  and  T is  the  integrate  time 
for  the  matched  filters.  Each  [X^(kT)}^  and  each  {Y^(kT)}^  is 
obtained  by  filtering  a sequence  of  independent,  unit  variance 
Gaussian  random  variables.  The  filter  used  is  a second-order 
recursive  digital  filter  equivalent  to  the  cascade  of  two  iden- 
tical single  pole  filters.  If  w(kT)  is  the  input  sequence  and 
z(kT)  is  the  output,  the  filter  is  defined  by 

z(kT)  =A1z(kT-T)  + A2z(kT-2T)  + A3w(kT)  (15) 

where 


= 2e 


-2nBT 


= IE  11 

N 2 „ 


-2ttBTx  3/2 


a+e-2nllT)l/2 


and  B is  the  Doppler  (double-sided)  bandwidth  of  the  path. 

This  filter  is  initialized  prior  to  forming  Z(T)  by  loading 
Z(0)  and  Z(-T)  with  two  Gaussian  random  variables  having  a cor- 
relation coefficient  given  by 


1 +e 


-ttBT 

-2ttBT 


The  received  power  is: 


rj(k)  = /xj(  k)  + yj(k) 


C 


r 


Thus,  overall, 


rcn(k) 


NW“  1.2 2 


^ JXj(k)+Y.(k)  cos4)jn(k)  + ncn(k) 


(21) 


and 


N 


r GO 
sir 


path  rr  y 

£ VXj(k)+Yj(k)  sil^jn(k)+,1s, 


00 


(22) 


The  decision  statistics  appropriate  to  this  processing  are: 


l = r (k)  r (k  - 1)  + r (k)  r (k  - 1)  (23) 

cn  cnv  ' cnv  ' sn  sn 


l = r (k)  r (k  - 1)  - r (k)  r (k  - 1) 
sn  sn'  ’ cn'  cn'  sn 


th 


(24) 


where  K>  (4  ) is  the  component  of  the  n tone  at  frame  number  k 

cn  sn 

in-phase  (quandrature)  with  the  n tone  at  frame  number  (k  - 1) • 


A-5 


; 


APPENDIX  B 


CODING  CONSIDERATION  FOR  VOICE  PROTECTION 


In  this  appendix,  some  insight  into  the  code  selection  for 
both  multiple-rate  codes  and  block  codes  is  given.  Unfortunately, 
mathematical  results  are  obtained  for  a simplified  HF  channel 
model  which  does  not  allow  for  correlation  between  code  digits. 
Thus,  even  though  firm  conclusions  cannot  be  drawn  from  these 
results,  trade-offs  in  performance  can  still  be  ascertained. 

Note  that  only  single-burst-correcting  and  double-burst-correcting 
multiple-rate  codes  are  considered.  The  technique  used  in  this 
appendix  can  be  applied  to  general  multiple-burst-correcting 
multiple-rate  codes;  however,  the  computation  becomes  increasingly 
tedious  and  lengthy  as  the  burst-correcting  capability  increases. 

B. 1 Performance  of  Single-Burst-Correcting  and  Double- Burst- 

Correcting  Multiple-Rate  Codes 

Performances  of  single-burst-correcting  and  double-burst- 
correcting  multiple  rate  codes  are  derived  under  the  assumption 
that  there  are  no  undetected  errors  and  that  blocks  are  indepen- 
dent. Both  the  high  rate  code  and  the  low  rate  codes  are  taken 
into  consideration.  The  performance  curves  for  some  multiple- 
rate  codes  suitable  for  41  tones  per  frame  format  are  computed. 

The  results  are  compared  to  the  Golay  (24,12;8)  code  and  the  low 
rate  code  employed  as  a block  code.  Only  the  hard  decoding  algor- 
ithm is  used  in  obtaining  the  results  for  multiple  rate  codes. 

B.1.1  General  Error  Probabilities 

Let  Pft  be  the  probability  that  a high  rate  code  block  con- 
tains uncorrectable  errors  and  P^  the  corresponding  probability 
for  the  low  rate  code.  When  bits  within  a code  block  are  inde- 
pendent, we  have 


•] 

C 


(1) 


ph = 


N 

E 

i=eh+1 


,Nv  sN-i 

(^P  (1  - P) 


and 


P.  = 


N 

E 

i=e 


(i>pid 


P) 


N-i 


l+l 


(2) 


where  e^  and  e^  are  the  respective  error-correcting  capabilities 
of  the  high  rate  code  and  the  low  rate  code,  and  p is  the  bit  error 
rate.  Let  a,  and  a be  the  conditional  bit  error  rate  for  the 

h t v 

high  rate  code  and  the  low  rate  code,  respectively,  given  the  code 
block  is  in  error.  For  the  case  where  bits  are  independent  within 
a code  block,  we  have 


and 


a,  = 


1 2 ,N.  . i,..  v N-i 

N . S+1  (i)ip  (1-p) 

1=eh+1 

N 

E 

i=e,  +1 


,N.  in  vN-i 
(t)p  (1  - P> 


a.  * 


1 2 ,N.  . iM  .N-i 

N . <1>1P  (1'P> 

1 = e^+1 

2 ,N.  i.-  -N-i 

S (,)P  (1  - P) 

i=e  +1  1 

'0 


(3) 


(4) 
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It  is  more  convenient  to  define 


sh  " “hph 


(5) 


and 


‘ a<.p«. 


(6) 


Then  /3^  and  /3^  are  given,  respectively,  by  the  numerators  of 
Eqs.  (3)  and  (4). 

B.1.2  Single-Burst-Correctina  Multiple-Rate  Code 

A single-burst-correcting  multiple-rate  code  can  be  char- 
acterized by  a two-state  Markov  chain  as  shown  in  Fig.  1.  The 
states  are  h and  l.  When  the  chain  is  in  state  h,  a code  block 
is  decoded  as  a high  rate  code;  when  it  is  in  state  l , a code 
block  is  decoded  as  a low  rate  code.  A state  is  determined  if  the 
previous  state  and  transition  probability  are  known.  The  transi- 
tion (probability)  matrix  for  the  single  Markov  chain  is 


(7) 


Let  w.  and  w. 
h l 

respectively. 


be  the  steady-state  distribution  for  states  h and  l. 
Then 


(8) 
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Solving  Eq.  (8),  we  have 


w,  = 

h 1 - P 


h + Pt 


(9) 


and 


1 - P, 


w£  1 - P 


h + Pf 


(10) 


Recall  that  the  multiple-rate  code  is  an  (NjK^+K^jD^)  -* 

(N, 2K^  + K2 ; D^)  code.  The  most  significant  bits  are  protected 
by  the  multiple-rate  code  structure;  and  the  second  most  sig- 
nificant bits  are  protected  by  either  the  low  rate  code  or  the 
high  rate  code,  depending  on  the  state.  For  the  most  signi- 
ficant bits,  an  erroneous  block  results  if  two  consecutive  blocks 
are  in  error;  thus,  the  bit  error  rate  for  the  most  significant 
bit  is : 


pbtt  * ”hphah  + ViVt 


■ whptA  + VA 


va  + (1  - ph>vt 


1 - ph  + pt 


an 


:l 


where  p^  and  P^  are  the  decoded  bit  error  rates,  respec- 
tively, for  the  high  rate  code  and  the  low  rate  code  when 
employed  as  block  codes.  At  high  signal-to-noise  ratios, 

— 0,  w ^ — 1;  thus, 


P 


bit 


Vt 


(13) 


and 


bit 


lb 


(14) 


B.1.3  Double- Burst-Correcting  Multiple-Rate  Code 

Whether  a block  is  decoded  as  a low  rate  code  or  as  a 
high  rate  code  depends  on  the  decoded  results  of  the  previous 
two  blocks  for  the  double-burst-correcting  multiple-rate  code. 
As  a result,  five  states  are  needed  to  describe  the  Markov 
chain:  Hh,  H-t,  Hh,  LI,  and  Lh.  The  second  letter  of  the  state 

shows  that  the  present  code  block  is  decoded  as  the  low  rate 
code  (l)  or  the  high  rate  code  (h) . The  first  letter  denotes 
the  conditions  of  the  previous  block: 
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Note  that  it  is  impossible  to  have  state  HI  and  Li  because,  when 
the  previous  block  is  known  to  contain  a burst,  the  present  block 
must  be  decoded  as  the  high  rate  code.  Also, it  is  impossible  to 
have  state  Lh  because  the  state  Lh  implies  that  the  previous  two 
blocks  are  decoded  correctly;  hence,  the  present  block  must  be 
decoded  as  a low  rate  code.  Therefore,  only  L is  allowed  to 
follow  L.  The  state  diagram  is  shown  in  Fig.  2.  Let  w^,  w^,  w^  , 
wt-,  , and  wT  be  the  corresponding  steady-state  distributions  for 


cfflfoo  T-TVi  Hi  T Vi  on  H T I Then 


Solving  Eq.  (15),  we  have 


”Hh  ■ V <1-ph>/ip(,  + a-V<1-ph+V1  (16> 

"nh  - Vl-(1-ph>2i/‘pt  + a- Vfl‘ph+VI  <17> 

"h t ‘ “Eh  ‘ p<,(i-V2/tp<,  + <1-V(1-ph+V1  (18) 
wu  = <1-V(1-ph>2/(pt  + <1-ph)<1-ph+V1  <19) 

Recall  that  the  most  significant  bits  of  blocks  containing  a burst 
can  be  correctly  recovered  if  they  are  followed  by,  or  adjacent  to, 
two  error-free  blocks.  In  Table  1,  all  possible  cases  resulting  in 
an  erroneous  block  being  present  are  listed  (a  "1"  indicates  the 
block  contains  a burst  and  a "0"  denotes  the  block  is  correct) . 

Since  there  is  a burst  at  time  0,  all  blocks  at  time  t must 
be  decoded  as  high  rate  codes.  Also,  due  to  the  error  patterns 
shown  in  the  table,  all  code  blocks  at  time  2t  and  3t 
are  also  decoded  as  high  rate  codes.  For  cases  1,  2,  and  3,  the 
joint  probability  time  t to  time  3t  is 

T = 2P2(1  - Ph)  + = 2P2  - P2  (20) 

For  cases  4 through  8,  the  joint  probability  from  time  r to  time 
3t  is 

S = Ph(l-Ph)2  + 3P2(1-Ph)+P2  = Ph(l-Ph)+T  (21) 
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The  bit  error  rate  for  the  most  significant  bits  is 


1 

I 

Pbit  ' (”HhPh“h  + wlUPt“t  + + <WHh  + “Lh>Ph“hS 

‘ [uHh^h  + (uH<,  + *UWt|T  + + “Ch^h5 

" [ (wHh  + "Sh  + "Lh)Sh  + ("Ht  + "H^t1  <2Ph  ' Ph) 

+ (wHh  + "Lh>f)hPh(1  ‘ Ph)  <23) 

The  bit  error  rate  for  the  K2  second  most  significant  bits  is 

Pbit  * (”Hh  + “iih  + "Lh)phb  + (wHt  + WLt)p{.b  <24> 

At  high  signal-to-noise  ratio, 


Comparing  Eqs.  (26)  and  (27)  to  Eqs.  (13)  and  (14),  we  observe 
that  is  about  the  same  for  both  single-  and  double-burst- 

correcting  multiple-rate  codes  at  high  signal-to-noise  ratios; 
but  for  the  double-burst-correcting  multiple-rate  code 

would  be  much  smaller  than  that  for  the  single-burst-correcting 
code  at  high  signal-to-noise  ratios. 

B.1.4  Performance  Curves  for  Multiple-Rate  Codes 

Performance  curves  for  single-  and  double-burst-correcting 
multiple-rate  codes  are  plotted  in  Figs.  3 through  8.  The 
relevant  information  concerning  each  of  the  figures  is  given  in 
Table  2. 

In  Table  3 the  performance  of  various  coding  schemes  are 
compared  at  a signal-to-noise  ratio  per  received  tone  (E^Ng) 
of  15  dB. 

B.1.5  Recommendations 

1)  When  a block  code  is  employed,  either  the  (64,36;12) 
code  or  its  shortened  version  may  provide  a better  performance 
than  the  two  Golay  (24, 12; 8)  codes  used  in  parallel  when  hard 
decoding  is  employed.  In  Section  B.2  further  comparison  of  these 
two  codes  will  be  presented. 

2)  The  (43,15 ; 12)  - (43,30;6),  (64,36;12)  - (64,51;6),  and 
(64,36;12)  -»  (64,45;8)  multiple  rate  codes  appear  to  be  attrac- 
tive for  voice  protection  when  hard  decoding  is  employed.  These 
codes  warrant  simulation  for  soft  decoding  performance. 
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TABLE  2 


is 


RELEVANT  INFORMATION  FOR  FIGURES  3 THROUGH  8 


Figure 

Single-  and  Double- 
Burst  -Correcting 
Multiple-Rate  Codes 

K1 

K2 

Block  Codes  Compared 

3 

(43,15; 12)  - (43 , 30 ; 6) 

15 

0 

(24,12;8)  & (43,15 ; 12) 

4 

(64,36;12)  - (64,51;6) 

15 

21 

(24,12;8)  & (64, 36; 12) 

5 

(64,36;12)  - (64,45;8) 

9 

27 

(24,12;8)  & (64, 36 ; 12) 

6 

(64,36; 12)  - (64,39;10) 

3 

33 

(24 , 12 ; 8)  & (64, 36 ; 12) 

7 

(82,54;9)  - (82,68;5) 

14 

40 

(24,12 ; 8)  & (82,54;9) 

8 

(82,54;9)  - (82,61;7) 

7 

47 

(24 , 12 ; 8)  & (82,54;9) 
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TABLE  3 

COMPARISON  OF  VARIOUS  CODING  APPROACHES  USING  HARD  DECODING  FOR  41  TONES/FRAME  FORMAT 
AT  A SIGNAL-TO-NOISE  RATIO  PER  RECEIVED  TONE  = 15  dB  FOR  A SLOWLY  FADING  HF  CHANNEL 
USING  4-PHASE  DPSK  MODULATION 


= Number  of  most  significant  bits,  K2  = Number  of  second  most  significant  bits, 
= Number  of  third  most  significant  bits 
MRC  = Multiple-Rate  Code 
SBC  = Single-Burst-Correcting  MRC 
DBC  = Double-Burst-Correcting  MRC 
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Figure  3 Performance  Comparisons  of  a Hard  Decoded  (43, 15; 12)  - 
(43,30;6)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones/Frame) 
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Figure  4 Performance  Comparisons  of  a Hard  Decoded  (64, 36; 12)  - 
(64,51;6)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones/Frame) 
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Figure  5 


Performance  Comparisons  of  a Hard  Decoded  (64,36;12)  - 
(64,45; 8)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones/Frame)  „ 

B-17 


a 


1 


p-p-pr 


1 I 1 I'TT 


= 3 Most  Significant  Bits 

K„  = 33  Second  Most  Significant 
Bits 


\ > . vv 
\ :• _ >, 

\\Y 

\V\ 

\ '1  v\ 


3L  o3 

£ c *-*• 

rt  rt 
rt  w » 
rt*  <■*  0 
^ O rti 

*t  ° 2 

73  * 2 

tu  ro 

^ ^ H 

o£  1 
o w 

O' 

TO 


\\  v\ 

i \ 

\ \ 

\ \ 

\ 

l\  Mi 


ft  r* 


ngw  N 

° S rt 

cu  ^ 
n>  * w 
cd  „ 
o o 

ft 

rj  tn 
oq  rt* 

K*  <0 

rt  • 
rt*  «) 

'a  £ 

W-* 

fl>  n 

i rt 

5» 

rt 

n> 


A* 


7*  V 

r-J  cr 


c*  ^ 

rt  01 

rt*  n 

xfi 

*-* 

y~>  *■*» 

o 

®»  <£. 

?r 

, 2s  S 

So  * 

\ ■,  o St1  l 

V?  " % 2 g 

YL  ^ » £ 

\\V>  o ? « 

\ V • o*  b ^ 


VS 

\\* 

\ \o 

V \° 

\ \o 
\ \ w 


n 

o 

rt 

rt 

rt 

o 

(0 

o 

O 

rt 

rt 

rt 

rt* 

ff> 

2 

o 

rt 

■ 3 

0Q 

|j_f j 

Ij.1.1 

Lj.J_i_Ll.Lj.  I.VlL.l  I.  l.'Pi.l  J.i.J lJ-.I_Ll.I_L  1 J J I 


TOT  PL  SNR  PER  RECEIVED  TONE  LET./ NO)  IN  DB 

Figure  6 Performance  Comparisons  of  a Hard  Decoded  (64,36; 12)  - 
(64, 39; 10)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones/Frame) 
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Figure  7 Performance  Comparisons  of  a Hard  Decoded  (82, 54; 9)  — 
(82,68;5)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones /Frame) 
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Figure  8 Performance  Comparisons  of  a Hard  Decoded  (82,54;9)  -* 
(82, 61; 7)  Multiple-Rate  Code  with  Block  Codes  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
(41  Tones /Frame) 


B-20 


B . 2 Performance  Comparisons  of  Shortened  and  Unshortened 
BCH  (64,36;12)  Codes  with  the  Golay  (24,12;8)  Code 
for  Voice  Protection 

In  view  of  the  closed  hard  decoding  performance  results 
obtained  in  Section  B.l  for  the  (64,36;12)  code  and  the  Golay 
code,  shortened  and  unshortened  BCH  (64,36;12)  codes  are 
evaluated  and  compared  to  the  performance  of  two  side-by-side 
Golay  ( 24 , 1 2 ; 8 ) codes.  The  BCH  codes  used  are  the  (52,24;12) 
code,  (55,27;12)  code,  (58,30;12)  code,  (61,33;12)  code,  and 
(64, 36; 12)  code. 

The  hard  decoded  results  are  shown  in  Fig.  9.  The  (58,30;12) 

code  performs  approximately  the  same  as  the  Golay  code  for  a 

-3 

decoded  bit  error  rate  of  10  ; but  the  BCH  code  provides  the 

protection  of  30  bits  vs.  24  bits  in  the  case  of  the  Golay 

code.  The  two  shorter  BCH  codes  outperform  the  Golay  code  for 

-3 

a decoded  bit  error  rate  of  10  , and  the  two  longer  BCH  codes 

protect  more  bits  (33  bits  and  36  bits  vs.  24  bits  for  Golay 
codes).  Therefore,  for  the  hard  decoding  case,  all  five  of  these 
codes  are  strong  competition  for  the  Golay  codes. 

For  the  soft  decoding  case,  none  of  these  five  BCH  codes 
are  expected  to  be  comparable  to  the  Golay  code,  as  shown  in 
Fig.  10.  Therefore,  when  soft  decoding  is  employed,  the  Golay 
code  would  probably  yield  the  best  performance.  Detailed  simu- 
lations for  realistic  HF  channels  are  still  required  before  firm 
conclusions  can  be  drawn  on  the  effectiveness  of  these  coding 
techniques  for  digital  speech  transmissions. 
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Figure  9 Performance  Comparison  of  Hard  Decoded  Shortened  and 

Unshortened  BCH  (64,36;12)  Codes  and  the  Golay  (24,12;8) 
Code  over  a Slowly  Fading  HF  Channel  with  4-Phase  DPSK 
Modulation 


Figure  10  Estimated  Performance  of  Soft  Decoded  Shortened  and  Unshortened 
BCH  (64,36;12)  Codes  and  the  Golay  (24,12;8)  Code  over  a 
Slowly  Fading  HF  Channel  with  4-Phase  DPSK  Modulation 
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APPENDIX  C 


SIGNAL-TO-NOISE  RELATIONSHIPS  FOR  HF  MODEMS 


It  is  quite  important  to  establish  proper  definitions  of 
signal-to-noise  since  this  report  compares  results  for  different 
modulation  formats  such  as  the  conventional  16-tone  modems  and  an 
ANDVT  modem  designed  with  39  tones.  The  modulation  format  for 
the  preamble  design  is  based  on  using  only  four  parallel  tones 
followed  by  a wideband  pulse  sequence.  Thus,  any  definition  of 
SNR  must  also  be  suitable  for  preamble  evaluation. 

In  all  cases  it  is  meaningful  to  obtain  results  in  terms  of 
the  ratio  of  the  total  received  power  to  noise  power  density. 

This  signal-to-noise  density  ratio  (SNDR)  is  given  by  the  ratio 
of  the  total  received  power,  P,  to  the  noise  power  per  Hz  of 
bandwidth,  denoted  by  Nq,  as 

SNDR  = ~ (1) 

0 

For  preamble  evaluation  in  Section  2 results  are  given  in  terms 
of  this  SNDR.  The  above  is  called  a signal-to-noise  density  ratio 
since  the  denominator  is  the  noise  density,  rather  than  the  more 
conventional  noise  power  in  the  bandwidth  of  the  signal.  To  avoid 
the  ambiguity  which  results  from  different  bandwidth  definitions 
one  may  define  the  signal-to-noise  ratio  per  information  bit 


N 


0 


(2) 
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where  R is  the  actual  data  rate  of  the  modem.  For  speech  trans- 
mission R = 2400  b/s  and  thus  Eqs.  (1)  and  (2)  can  be  related  by 


+ 33 . sjdB  - for  R = 2400  b/s  (3) 

For  digitized  speech  the  modulated  signal  is  approximately 
contained  in  a bandwidth  of  2400  Hz, and  thus,  E^/Nq  can  also  be 
viewed  as  the  ratio  of signal-to-noise  power  measured  in  a band- 
width of  2400  Hz. 

In  addition  to  obtaining  results  in  terms  of  P/Nq  and 
Eb/Nfl,  which  are  easily  related,  it  is  customary  for  parallel  tone 
evaluations  to  present  results  in  terms  of  the  average  signal-to- 
noise  ratio  per  tone  defined  by 


N 


0 


E[n2]  + E[n2] 
c s 


(4) 


where  r is  the  tone  amplitude  and  n and  n are  the  cosine  and 

c s 

sine  noise  components  at  the  output  of  the  integrate-and-dump 
circuits.  This  SNR  is  also  closely  related  to  the  signal-to- 
noise  power  measured  in  a bandwidth  equal  to  that  of  the  modulated 
signal . 

To  relate  E /N„  to  E,/Nn  we  first  note  that  a multiple  tone 
to  b U 

received  signal  can  be  written  as 

tl 

y(t)  = L r.  cos  (o;.t  + 0.  ) + n(t)  (5) 

1=1  1 


where  n(t)  is  assumed  to  be  a white  Gaussian  noise  process 
normalized  such  that  E[n(t)n(t  + t) ] = (Nq/2)6(t). 

The  cosine  and  sine  components  per  tone  are  given  by 

T 

y = 2 J y(t)  cos  co.  t dt  = r.  cos  0^  + nc 

ci  o 11  Ci 

(6) 

T 

y = 2 f y(t)  sin  to.  t dt  = r.  sin  6.  + n 

•'s,  i i i s. 

1 0 1 

Note,  the  multiple  tones  are  orthogonal  over  the  integration  time 
T,  which  is  less  than  a baud  interval  by  the  allowable  guard 
time . 

The  average  signal  energy  per  tone  is  given  by 

E.  = E[r?  cos^  0.  + r?  sin^  0.]  = E[r?]  (7) 

t l li  l i 


and  the  energy  per  noise  component  is  given  by 
E[n^  ] = E[n^  ] = NQ/2 
The  total  received  power  can  be  written  as 


(8) 


(9) 


where  N is  the  total  number  of  tones  and  T is  the  integration  time 
per  tone.  Thus,  we  can  relate  (1)  and  (2)  to  Ej-/Nq  by  the  following 


03 


ft 


E. 


N, 


N “t 
T N 


(10) 


0 


and 


N 


0 


JL  _t 

TR  N 


(11) 
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For  a conventional  16  tone  modem*  operating  at  2400  b/s,and 
including  a Doppler  correction  tone  of  twice  the  amplitude, 
^effective  = ^0)  » equation  (11)  gives  the  relationship  between 


E, /N^  and  E^/N^  as 
b (J  t U 


20 


N, 


1 

110 


2400  \N0; 


(12) 


or 


N 


0 


— _ 37 

N 

0 


dB  - 16  tone  modem 


(13) 


For  an  ANDVT  modem  operating  at  2400  b/s  with  two  Golay 
codes,  using  39  tones  and  an  integration  of  18  ms,  we  have 


39 


N0  18  x 10"3  x 2400  W 


(14) 


D.  Chase,  "A  Combined  Coding  and  Modulation  Approach  for  Com- 
munication Over  Dispersive  Channels,"  IEEE  Trans,  on 
Communications . Vol.  COM-21,  No.  3,  March  1973,  page  159. 
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39  tone  ANDVT  modem 


(15) 


Thus,  in  this  case,  a comparison  based  on  Et/N  is  essentially 
equivalent  (adjustment  is  only  .07  dB)  to  the  fiar  comparison 
based  on  E^/Nq. 


If  the  ANDVT  modem  is  used  without  coding  a data  rate  of 
2400  bps  can  be  achieved  with  27  tones  which  are  keyed  every 
22.5  ms  with  one  of  four  phases  (2  bits/tone).  In  this  case 
the  relationship  between  E^/Nq  and  Et/Ng  is  given  by 


27 

18  x 10'3  x 2400 


(16) 


'b 

0 


1.04 


dB 


27  tone  ANDVT  modem 


(17) 


Note,  by  comparing  Eqs.  (13)  and  (17),  without  coding  the  ANDVT 
modem  offers  a 1.67  dB  advantage  in  E^/Ng  when  compared  to  the 
16-tone  modem.  This  gain  is  due  to  the  fact  that  the  loss  due 
to  guard  time  is  smaller  for  the  ANDVT  modem  and  also  due  to  the 
Doppler  tracking  algorithm  for  the  ANDVT  modem  which  does  not 
require  an  external  tone. 

The  ANDVT  modem  results  are  all  presented  in  terms  of  E^/Ng 
and  the  preamble  results  are  presented  in  terms  of  P/Nq. 

Equation  (3)  results  relate  these  two  measures.  The  results  for 
the  KG  synchronizing  sequence  are  presented  in  terms  of  Et/Ng, 
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which  has  the  desirable  property  that  guard  time  and  code  redun- 
dancy does  affect  this  SNR  measure.  As  a consequence  of  the  low 
data  rate  used  for  the  KG  sync  sequence,  E^/Nq  is  considerably 
greater  than  E^/Nq  for  this  case,  and  thus,  it  is  more  meaningful 
to  make  comparisons  based  on  P/Nq.  Equation  (10)  relates  P/Nq 
to  E^/Nq,  which  for  a 39-tone  preamble  and  integration  of  18  ms 
yields 


= [E^/Nq  + 33.58]  dB  - 39-tone  KG  preamble  (18) 

In  summary,  the  preamble,  KG  sync  sequence,  and  multiple 
tone  modems  can  be  compared  in  terms  of  P/Nq,  but  by  convention 
E^/Nq  and  E^Nq  are  sometimes  used.  Equation  (3)  provides  the 
appropriate  conversion  for  the  ANDVT  modem  results  given  in  Sec- 
tion 4,  and  Eq.  (18)  provides  the  conversion  for  the  KG  sync 
results  given  in  Section  3. 


APPENDIX  D 


BIT  AND  BLOCK  ERROR  RATE 
COMPUTATION  FOR  LINEAR  BLOCK  CODES 


A detailed  analyses  for  obtaining  a lower  and  an  upper  bound 
on  decoded  bit  error  rates  for  linear  block  codes  are  presented. 
The  tightness  of  these  bounds  are  discussed.  An  optimistic 
approximation  to  the  bit  error  rate  of  soft-decoded  linear  block 
codes  is  then  given.  Finally,  the  block  error  rate  expression 
for  hard-decoded  linear  block  codes  and  an  approximation  to  soft- 
decoded  block  error  rate  is  obtained. 

D. 1 A Lower  Bound  on  Hard-Decoded  Bit  Error  Rate 
for  Linear  Block  Codes 

Let  denote  the  decoded  bit  error  rate  for  an  (N,K;d) 
block  code  and  let 

e = L(d-1)/2J  (1) 

Then  the  expected  number  of  errors  per  decoded  block  is  Np^- 
This  expected  number  of  decoded  errors  can  be  lower  bounded  by 


„ ^ d"1  ,N.  . i/T  .N-i  , , dJe  ,N,  i ...  v N-i 

Np,  > L (.)ip  (1-p)  +d  £ (. )p  (1-p) 

i=e+l  i«d 


N 

N i N-i 

+ £ (”)(i-e)p1(l-p)N  (2) 

i = d+e+1 
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where  p is  the  raw  bit  error  rate.  This  lower  bound  is  obtained 
by  first  noting  that  all  possible  patterns  of  e or  fewer  errors 
are  corrected  since  the  block  code  is  an  e-error  correcting  code. 
Error  patterns  with  e+1  to  d-1  errors  are  assumed  to  be  detected, 
thus  no  additional  errors  are  introduced  by  the  decoder.  When 
there  are  d to  d+e  errors,  a lower  bound  is  achieved  by  assuming 
that  these  patterns  are  all  decoded  into  d errors.  Finally,  if 
there  are  i > d+e  errors,  the  decoder  is  assumed,  optimistically, 
to  leave  i-e  decoded  errors.  Therefore,  Eq.  (1)  is  strictly  an 
inequality.  Dividing  both  sides  by  N,  Eq . (2)  becomes 


pb 


d-1 

E 

= e+1 


,Nv.  i ,,  vN-i 
(i)xp  (1-p) 


+ N 


d+e 


i=d 


,N.  i/n  v N-i 

(t)P  (1-P) 


N 

E 

d+e+1 


,NW.  . i,,  vN-i 
(i)(i-e)p  (1-p) 


(3) 


which  provides  the  lower  bound  on  the  decoded  bit  error  rate  for 
hard  decoding. 

D. 2 An  Upper  Bound  on  Hard-Decoded  Bit  Error  Rate 
for  Linear  Block  Codes 

In  this  section,  an  upper  bound  on  the  hard  decoded  bit 
error  rate  for  linear  block  codes  will  be  obtained.  This  bound 
is  independent  of  the  weight  distribution  of  a block  code.  The 
tightness  of  the  bound  is  to  be  discussed. 
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A 


For  the  case  that  d is  odd,  the  simple  bound  is: 


PK  < 77  £ (^)  (i+e)pL(l-p)N  1 for  d odd 

1 = e+1 


This  can  be  seen  as  follows:  whenever  there  are  e+1  or  more 
errors,  the  decoder  is  assumed  to  produce  extra  e errors.  How- 
ever, the  decoder,  in  reality,  produces  e or  fewer  errors.  Thus 
Eq.  (4)  is  indeed  an  upper  bound.  It  can  be  seen  from  Eq . (3) 
and  Eq.  (4)  that 

(p,  ).  , , > (p,  ) , , for  d odd  (5) 

rb  lower  bound  d b upper  bound 

In  other  words,  for  d odd, 


d ^b^upper  bound  ' ^b  < ^b^upper  bound 


where 


1 ^ d^e  = e+1  ^ 2 

2 d 2 e+1  - 3 


When  the  number  of  code  words  of  minimum  weight  d is  known  to 
be  w^,  the  upper  bound  given  in  Eq . (4)  can  be  tightened  to: 


Pk  < 


, N . e+1.,  vN-e-1 

Wp  (i-p) 


I 


where 


‘5-d'>0 


(9) 


Note  that  when  there  are  e+1  errors  in  a code  word,  there  is  a 
probability  A that  the  pattern  is  not  detected,  which  results 
in  d decoded  errors,  and  a probability  1 - A that  the  pattern 
is  detected  and  the  e+1  errors  remain  unchanged.  The  expression 
for  A is  obtained  by  noting  that,  for  each  code  word  of  minimum 
distance  d away  from  the  transmitted  code  word,  ( ^)  patterns 
of  e+1  errors  are  undetected.  [There  are  ( ^)  possible  error 
patterns  of  weight  e+1.]  If  Eq . (8)  is  employed,  the  upper  bound 
is  very  tight,  since  the  first  term,  which  is  the  dominating  term 
for  small  values  of  p,  is  the  exact  expression  for  the  decoded 
bit  error  rate.  Unfortunately,  for  most  linear  codes,  the  value 
of  w^  is  unknown.  Therefore,  Eq . (4)  is  employed  to  obtain  all 
of  the  performance  curves  in  this  report. 

For  the  case  that  d is  even,  the  upper  bound  on  the  decoded 
bit  error  rate  is  given  by 


e+1  , N . e+1 . 1 . N-e-1  1 

pb  - ~ WP  (1-p)  + N 


for  d even 


N 

£ 

i = e+2 


(J)(i+e)pi(l-p)N"i 


(10) 


The  first  term  of  Eq.  (10)  is  due  to  the  fact  that  all  possible 


e+1  = ^ errors  can  be  detected  for  even  d 


Equation  (10)  is  as 
tight  as  Eq . (8)  for  odd  d.  However,  there  is  no  need  to  know 
the  value  of  w^  in  evaluating  Eq . (10),  since  the  dominating  term 
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of  Eq . (10)  is  identical  to  the  dominating  term  of  Eq . (3),  the 
lower  bound.  Thus,  for  small  values  of  p,  the  upper  and  lower 
bounds  are  almost  indistinguishable.  [Also,  Eq . (6)  can  be 
easily  shown  to  be  valid  for  even  d when  Eq.  (10)  is  employed 
to  compute  the  upper  bound.] 

D. 3 An  Estimate  of  Soft  Decoding  Performance  for  Block  Codes 

Unlike  the  case  of  hard  decoding,  the  exact  performance  of 
soft-decoded  block  codes  cannot  be  bounded  accurately.  Note  that 
the  performance  depends  not  only  on  the  code  structure  and  the 
raw  bit  error  rate  but  also  on  the  more  detailed  channel  char- 
acteristics, such  as  the  confidence  level  of  each  received  bit. 
However,  a less  accurate  estimate  than  those  obtained  for  the 
hard  decoding  can  be  easily  computed  by  assuming  that  all  the 
possible  hard  error  patterns  of  weight  d-1  or  less  can  be  cor- 
rected. Since  a soft  decoder  can.iusually  only  correct  a fraction 
of  hard  error  patterns  with  weight  d-1  or  less,  this  estimate 
obtained  for  the  soft  decoding  may  be  somewhat  too  optimistic 
for  some  codes.  Nevertheless,  this  is  a simple  and  easy  way  to 
obtain  an  approximation  to  lower  bounds  for  the  code  performance. 
Furthermore,  asymptotically,  an  optimum  soft  decoder  can  correct 
(d-1)  errors  [1].  [Note  that  for  a more  conservative  estimate 
of  the  soft  decoding  performance  the  code  is  assumed  to  be  able 
to  correct  d-i  or  fewer  hard  errors  with 

1 < 1 < (ID 


A properly  chosen  value  of  i may  yield  a good  estimate  for  the 
soft  decoding  performance  for  a meaningful  range  of  signal-to- 
noise  ratios.  One  way  to  determine  i is  to  simulate  the 
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performance  by  computer  at  relatively  low  signal-to-noise  ratios. 
Once  the  value  i is  determined,  it  is  rather  simple  to  arrive  at 
the  estimate  of  soft  decoding  performance  at  higher  signal-to- 
noise  ratios.] 

D.4  Block  Error  Rate  for  Linear  Block  Codes 

For  hard  decoding,  a block  error  results  if,  and  only  if, 
there  are  e+1  or  more  errors  per  code  block.  Thus  the  block 
error  rate  is  given  by 


L 


N 


P,  = £ (^)pi(l-p)N'i 

i=e+l 


(12) 


For  small  values  of  raw  bit  error  rate  p. 


d ~ / N v e+1...  .N-e-1  N 

Pb  “ <e+l)p  (1-p>  " Zn  Pb 


(13) 


An  optimistic  estimate  for  soft-decoded  block  error  rate  is 
given  by 


N 


n /N.  i...  .N-i 

pb  * * <i)p  (1_p) 

i=d 


(14) 


which  is  obtained  by  assuming  that  all  d-1  or  fewer  hard  errors 
can  be  corrected.  At  high  signal-to-noise  ratio 


N.  d N-d  N 
pb  ~ (d)P  (i-p)  ~ d pb 


(15) 
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APPENDIX  E 
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DOPPLER  TRACKING  DURING  DATA  TRANSMISSION 


E.  1 Introduction 

In  order  to  track  out  the  Doppler  offset  while  data  is 
being  transmitted,  a Doppler  discriminator  is  needed  to  pro- 
vide an  appropriate  error  signal  for  driving  the  voltage- 
controlled  oscillator  that  tracks  the  Doppler.  Conventional 
approaches  use  an  unmodulated  tone  and  conventional  frequency 
discriminator  for  a Doppler  discriminator.  This  approach  wastes 
channel  capacity  and  power  and,  in  addition,  does  not  take 
advantage  of  the  signal  diversity  present  in  the  received  wave- 
form. In  this  appendix  we  shall  discuss  a Doppler  discrimin- 
ator that  removes  these  objections.  However,  this  discrimin- 
ator introduces  new  requirements;  the  frame  sync  is  assumed  to 
be  operating  and  the  average  error  rate  needs  to  be  sufficiently 
below  1/2. 

f 

E. 2 Formulation  of  Received  Signal  and  Decision  Variables 

The  transmitted  signal  consists  of  a set  of  frequency- 
division-multiplexed  subcarriers,  each  carrying  4-DPSK  data. 
Thus,  the  k'th  modulated  subcarrier  can  be  represented  in  the 
(complex  envelope)  form 


dk(t)  =eJ2’kFt£eJ,>WRact(^) 


where  F is  the  frequency  spacing  between  subcarriers,  T is  the  pulse 
duration. 


Rect 


(!)  = 


i ; W * f 


0 ;|t|>  j 


and,  in  accordance  with  the  4-DPSK  format, 

j<Pk,t+l  11 

6 6 ~ Jl  ~ 2J2 


The  received  signal  corresponding  to  the  transmission  of  d^(t) 
will  be  distorted  by  the  multipath  in  the  vicinity  of  the  transitions 
between  pulses.  However,  in  the  receiver,  a synchronized  time  gate 
passes  only  the  central,  nondistorted  portion  of  each  pulse.  Any 
fading  present  at  the  subcarrier  location  still  appears  on  the 
received  gated  pulse.  Thus,  if  T(f,t)  is  the  time-variant  transfer 
function  of  the  channel,  the  complex  fading  at  the  subcarrier  fre- 
quency kF  is  given  by  T(kF,t).  Thus,  after  the  time  gate,  the 
received  modulated  subcarrier  (in  the  absence  of  noise)  can  be 
represented  in  the  form 


wk(t)  = T(kF,t)  e 


j 2frkFt 


Rect 


(H?) 


where  T is  the  duration  of  the  time  gate, 
g 
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The  received  noise  after  the  time  gate  is  given  by 


i(t)  -ERect 


where  n(t)  is  the  received  noise  before  the  time  gate. 

The  integrate-and-dump  operation  for  the  p baud  of  the 
kth  subcarrier  is  given  by 


pT+T  / 2 


pT-T  /2 
g 


-j2i7kFt  lw(t)+fn(t)]  dt 


where 

w(t)  * ^ wk^t^  ^ 

where  K is  the  number  of  subcarriers. 

Because  the  subcarrier  separation  frequency  F is  set  equal  to 
the  reciprocal  gate  duration, 


v - ± (8> 

Tg 

the  integration  in  (6)  separates  out  a signal  contribution  due  to 
the  p^  band  and  k^"  subcarrier  above,  to  wit, 

pT+T  / 2 pT+T  /2 

C g r 8 -j2irkFt  ... 

S = j T(kF,t)  e KP  dt  + f e n(t)  dt 

RP  dT-T  /2  pT-T  / 2 

PT  V 8 (9) 


T 


This  may  be  expressed  in  the  form 

S.  -e^u. 

kp  kp 


(10) 


where 


U.  = T.  + N. 
kp  kp  kp 


(11) 


in  which 


and 


/T+T  / 2 

g T(kF,t)  dt 

pT-T  U 
g 


pT+T  /2 

m _ f 8 -j2irkFt  J<pkp  . s 
NRp  J e J e v n(t)  dt 

pT-T  /2 
g 


(12) 


(13) 


Assuming  white  complex  Gaussian  noise  and  complex  Gaussian 

fading  statistics,  T,  , N,  , and  U.  are  complex  Gaussian  variables 
b ’ kp  kp  kp 

with 


TLN,  - 0 


cp  rs 


(14) 


Ng  N 


p rs 


0 ; k / p,  r ^ s 


2N0Tg  ; k-p,  r-s 


(15) 


where  is  the  real  one-sided  noise  power  density. 
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Under  the  reasonable  assumption  that  T(kF,t)  will  change 
little  over  one  gate  duration,  we  see  that 


Tkp  * T(kF>PT> 


The  SNR  at  the  output  of  the  integrate-and-dump  is 


1 KN0Tg 


where  P is  the  average  received  signal  power  for  all  tones  and  K is 
the  number  of  tones. 

The  decision  variable  used  in  4-DPSK  modulation  is 


Dkp  Skp  Sk,p+1  \p  + ^ ^kp 


j (Vi.  i _<Pj.  n) 

= e k,pfl  k»P  u*  U.  .. 

kp  k,p+l 


where  R^  and  are  the  real  and  imaginary  part  of  , respectively. 

Let  the  input  quaternary  data  be  encoded  into  two  binary  variables 
associated  with  the  sign  of  the  real  and  imaginary  parts  of 

exp(j<Pk  .p+i'^.p'  lsee  <3)|-  Thus> 

3^ktpfl  ^k,p3  1 . . . no. 

e ’ <rkp  J qkp)  (l9) 

where 


r,  = + 1 
kp  - 

q,  = + l 

kp  - 


Two  pairs  of  binary  decisions,  r^p  and  q^,  are  made,  as  follows: 


rkP  ‘ SsntRkp] 


Ikp  ' S*nl%>> 


where 


Sgn(x)  = 


1 x > 0 


-1  x * 0 


Under  the  ideal  conditions  of  very  slow  fading  and  no  noise. 


E,  = U*  U.  . - 
kp  kp  k,p+l 


|T(kF,pT)|  ; slow  fading  and 


no  noise 


and  we  see  that 


rkp  rkp 


qkp  qkp 


i.e.,  no  errors.  More  generally,  Ekp  has  a nonzero  phase  angle  akp> 


E.  = E.  e 
kp  kp 


and 


-Sgn[rkp  cos“kp-\p  sln“kp] 

(29) 

= Sgn[rkp  sin  \p+  qkp  coSo,kp] 

(30) 

The  phase  a,  is  a random  variable  due  both  to  the  additive  noise  and 

Kp 

the  channel  fading.  In  addition,  it  has  a nonzero  average  value  due 
to  the  presence  of  nonzero  instantaneous  frequency  of  the  received 
subcarrier.  Only  if exceeds  45°  will  errors  occur. 

Using  (16)  and  noting  that  the  noise  components  of  U,  and 

kp 

U,  , are  statistically  independent  (15),  we  can  show  that  the 

k,  p+l 

mean  value  of  is  given  by 


ss  2n  v T 


(31) 


where  v is  the  uncompensated  Doppler  shift  of  the  received  subcarrier. 

A typical  value  of  residual  Doppler  error  may  be  1 Hz.  With  T = 22.5  ms, 
this  produces  a value  of  ot^  of  0.14  radian,  or  around  8°. 

E. 3 Doppler  Discriminator 

The  proposed  Doppler  discriminator  utilizes  the  decision  variables 

R.  and  Q,  . In  essence,  a quantity  proportional  to  the  Doppler  error 
kp  kp 

is  estimated  on  each  pair  of  bauds  for  all  subcarriers  and  averaged 
over  several  successive  pairs  of  bauds.  The  basic  estimate  for 

a single  pair  is  given  by 


kp 


- ('V  - |RkpO  se-tKkpJ  s8"[Qkpl 


kp 


(32) 
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Let  us  analyze  this  estimate.  Using  (28),  (18),  (19),  (22),  and 
(23), 


= -4.  Ek  (i 

72  kP\ 


kP  72 


rkP  sln  “kP+ qkP  cos  “kp1  - |rkP  cos  % - qkP  sln  \p^kP  qkP 

(33) 


Assuming  that  no  error  has  been  committed, \ |<45°,  we  see  that 


r,  sin  a,  + q,  cos  a. ■, 
kp  kp  kp  1] 


rkPslnakp+qkPcosakP  : qkP>0 


' ‘rkPsln “kp‘qkPcos“kp;  qkP<0 


(34) 


r.  cos  a,  - q,  sin  a, 
kp  kp  kp  Rp 


r.  cos  a,  - q.  sin  a,  ; r,  >0 
kp  kp  kp  kp  kp 


-t,  cos  a,  + q,  sin  a.  ; r.  <0 
kp  kp  kp  kp  kp 


(35) 


Using  (34)  and  (35)  in  (33)  and  assuming  no  error  in  binary 
decisions,  we  find  that 


e,  = 72  E.  sin  a, 
kp  kp  kp 


(36) 


In  the  absence  of  noise  and  assuming  a reasonable  Doppler  error, 


e,  « 

kP 


72  |T(kF,pT) | |T(kF, (p+1)t)|  2n T u 


(37) 
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Summing  over  all  frames  we  obtain 


Note  that  weak  subcarriers  have  a corresponding  small  contribution 
to  the  error  signal.  The  actual  error  signal  used  would  involve 
a time  average  of  over  the  fading  with  an  appropriate  discrete 
filter, 

((>p=Sh4  (39) 

where  [h  } is  the  discrete  filter  impulse  response. 

A detailed  analysis  of  the  effects  of  bit  errors  has  not  been 
conducted  due  to  time  limitations.  At  this  point  it  is  believed 
that,  with  reasonable  amounts  of  averaging  time,  high  error  rates, 
well  beyond  those  at  which  intelligibility  has  been  lost,  can  be 
tolerated. 


